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ABSTRACT 



This experimental study was conducted at the Naval Post- 
graduate School to investigate the heat transfer character- 
istics of flat vertical plates in a rectangular gas - f luidi zed 
bed. The primary objective was to determine what effect 
variations in the bed width- to-height ratio had on heat trans- 
fer to the vertical flat plates forming the container walls. 

The experiment was conducted using a specially heated and in- 
strumented fluidized bed equipped with a movable side wall 
which permitted modification of the bed geometry. As the 
width of the bed was adjusted, the settled bed height was main- 
tained at a constant level by the addition or removal of bed 
material . 

A secondary objective of the study was to determine the 
effect of variations in the fluidization gas flow rate on heat 
transfer to the bed walls. Flow rates ranging from fixed bed 
fluidization to pneumatic conveying were studied. 

Pressure drop measurements as well as visual observations 
were used to determine minimum fluidization flow rates for 
each configuration. Heat transfer coefficients were calculated 
for each flow rate and bed geometry using temperature data ob- 
tained from a computer controlled thermocouple network. In 
addition, a study was made of variations in wall temperature 
with changes in vertical position. 
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I. INTRODUCTION 



A fluidization bed in its simplest form consists of a con- 
tainer filled with small particles which are subjected to an 
upward flow of fluid emerging from the bottom of the container. 
When the upward viscous force (drag) exerted on the particles 
by the moving fluid balances or exceeds the downward force 
of gravity, the particles become suspended and begin to move 
about the container in a fluid-like manner, hence the term 
"fluidization." This churning mixture of particles and gas 
is often referred to as the "emulsion phase." When the 
velocity of the fluid is increased to the point that the 
particles are carried out of the container, a transition is 
made to a different form of fluidization known variously as 
pneumatic conveying, entrained flow fluidization, or carryover. 
When the velocity of the fluid is too low to counteract the 
gravitational force and cause fluidization, the bed is referred 
to as being "packed" or "fixed" and the process is known as 
flow through a porous medium or fixed bed fluidization. 

Distinctly different modes of behavior are observed depend- 
ing on whether the fluid passing through the container is a 
liquid or a gas. Liquid fluidized beds are often referred to 
as homogeneous or particulately fluidized beds. This type of 
bed tends to display a very smooth transition through the 
various stages of fluidization from fixed to flui'dized to 
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carryover. Concurrently a very smooth transition in heat 
transfer characteristics and hydrodynamic behavior is observed. 
On the other hand, gas fluidized beds tend to display a more 
erratic performance. Most researchers attribute the difference 
in behavior to the differing particle to fluid density ratios. 
Gas fluidized beds are known as "aggregate" or "bubbling 
fluidized beds." In this type of bed, once fluidization has 
begun, the gas usually forms into bubbles which displace the 
bed particles as the bubbles rise. The appearance of a bubbling 
bed is remarkably similar to that of a boiling liquid, as shown 
in Figure 1. During conditions of very active fluidization 
the bubbles sometimes form together or coalesce to create a 
large bubble whose diameter approaches the width of the bed, 
as shown in Figure 2. This condition is known as "slugging" 
or "heaving" and can significantly effect bed performance. 
Another phenomenon of gas fluidized beds is shown in Figure 3. 

It is a condition known as "channelling." When this occurs, 
the gas travels through a fairly confined area forming a gas 
filled tube within the bed rather than passing evenly through- 
out the cross section. This behavior is believed to be 
caused by poor distributor design and, like slugging, can 
strongly effect the heat transfer characteristics of the bed. 

Fluidized beds offer a wide variety of beneficial char- 
acteristics such as high rates of heat transfer, high thermal 
inertia, large solid to fluid contact area, and isothermal 
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bed temperatures. In addition, once fluidized, the pressure 
drop across a bed remains relatively constant despite varia- 
tions in fluid flow rate. The pressure drop is a function 
of bed height and particle density and is independent of 
particle size. [Ref. 1] This permits increasing the surface 
area by using smaller particles without paying a penalty in 
increased pressure loss. 

Fluidized beds also suffer from certain limitations. The 
size of particles which can be fluidized is limited to a 
specific range extending from approximately 0.001 inch in 
diameter to about 0.25 inch diameter, depending on the density 
of the particles and the fluid. [Ref. 2] In addition, some 
materials cannot be fluidized regardless of size because of 
their tendency to clump together and form particles too large 
to fluidize, or because of their tendency to break down into 
particles which are so small that they become entrained in the 
fluidizing medium and are carried from the bed. Another dis- 
advantage is that fluidization requires the expenditure of 
energy to force the fluidizing medium through the bed at the 
required velocity. 

In general the beneficial aspects of fluidized beds far 
exceed the limitations and as a result, fluidized bed prin- 
ciples have been utilized in many diverse applications 
over the years. Early mining engineers used fluidization 
to separate solids of different densities by suspending 
them in rising streams of water. [Ref. 1] Prior to World 
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War II, chemical engineers discovered that fluidized beds 
made outstanding chemical reactors because of the extremely 
large contact area between the particles and the fluid. 

For example, a container measuring one foot by one foot by 
one foot, filled with 0.01 inch diameter spherical particles 
yields approximately 4,000 squares feet of surface area. 

This characteristic of fluidized beds has been utilized to 
great advantage in a wide number of chemical processes , most 
notable the catalytic cracking of petroleum. [Ref. 3] 

Because the particles within a fluidized bed are constantly 
churning about and mixing, the temperature at one point in 
the bed is usually in very close agreement with the tempera- 
ture at any other point. Industry has capitalized on this 
isothermal property by using fluidized beds to improve the 
heat treatment of metals in such processes as quenching, 
tempering, and carburizing. As a combined result of the 
large surface area and the thorough mixing of the particles, 
fluidized beds possess impressive heat transfer capabilities. 
Consequently, the power generation industry has begun utiliz- 
ing this technology to create a new generation of heat 
transfer devices such as steam generators which efficiently 
burn poor quality fuels while simultaneously yielding very 
low levels of harmful emissions. These devices are capable 
of burning high sulfur coal, pelletized garbage, and even 
raw crushed shale, just as it comes from the ground. Another 
application is seen in waste heat recovery units which are 
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fluidized by the exhaust gases from conventional boilers or 
gas turbines. The thermal energy which would otherwise 
have been wasted is transferred to a secondary fluid cir- 
culating through tubes emersed within the bed. This re- 
covered energy is then utilized to increase overall plant 
efficiency. [Ref. 4] 

The use of fluidized bed technology has expanded greatly 
in recent years as hundreds of new applications have been 
discovered, ranging from advanced ship propulsion plants to 
sophisticated solar energy conversion and storage devices. 

As the utilization of these devices increases so does the 
quest for a better understanding of the many complex phenom- 
enon associated with fluidization. 

Fluidized bed research has been conducted throughout the 
world and more has been written about fluidization than can 
be read in a lifetime, however much remains to be learned 
in order to unlock the full potential of this fascinating 
process. This thesis represents a modest attempt to add to 
the knowledge of fluidization in the area of heat transfer 
to vertical flat plates in gas fluidized beds. 
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II. OBJECTIVES AND METHOD OF APPROACH 



A. GENERAL STATEMENT OF THE PROBLEM 

The bulk of fluidized bed heat transfer research conducted 
to date has been related to the study of heat transfer to or 
from objects immersed within the bed, particularly tubes and 
tube bundles. Comparatively little research has been done 
on the study of heat transfer to the walls of the containing 
vessel itself. In addition, a large majority of the test 
apparatus utilized thus far to collect experimental data 
have been cylindrical in configuration. Very few studies 
have been conducted on heat transfer within rectangular 
fluidized beds; a geometry which is being used with increasing 
frequency in a wide variety of applications. 

The principle objective of this thesis was to gain a 
greater understanding of the heat transfer to the flat ver- 
tical plates which form the walls of a rectangular gas 
fluidized bed and to determine what effects, if any, varia- 
tions in bed geometry had on heat transfer performance. 

A secondary objective was to study the variations in heat 
transfer which occurred as the fluidization gas flow rate 
was adjusted over a wide range. 

B. METHOD OF APPROACH 

This investigation was primarily experimental in nature 
with limited utilization of analytical methods for comparison 
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purposes. To acquire the necessary data for this study, a 
computerized data acquisition system was linked to a network 
of up to 75 copper-constantan thermocouples installed in a 
specially instrumented and heated fluidized bed. This system 
was used to obtain a thermal "snapshot" of the test apparatus 
which showed the temperature at numerous key locations at 
a specific point in time. This data was then combined with 
measured values of input power and incoming air flow rates 
to produce a calculated heat transfer coefficient. Experi- 
mental data was obtained for a number of different air flow 
rates and bed conditions beginning with fixed bed fluidiza- 
tion and progressing through to the onset of entrained flow 
fluidization. Once data had been collected for the entire 
range of flow rates, the bed geometry was altered by moving 
an adjustable side wall within the test apparatus so that 
the width of the rectangular bed was increased. Additional 
particles were then added to the bed so that the bed depth 
remained at the same level as in the previous configuration. 

The data collection process was then repeated in its 
entirety for the new configuration. A total of four differ- 
ent bed geometries were studied in the course of this work. 
Over 5,700 temperature readings were recorded for eighty dif- 
ferent test runs. All experimental data, including narrative 
comments regarding observed bed behavior were recorded on 
magnetic disks. The recorded data was automatically pro- 
cessed by computer and then printed out in such a manner 
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that bed performance for a particular test run was readily 
apparent. The calculated results for the individual runs 
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on to the heat transfer investigation, an ex- 
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III. EXPERIMENTAL APPARATUS AND PROCEDURES 



A. DESCRIPTION OF FLUIDIZATION APPARATUS 

The fluidization apparatus consisted of a rectangular 
plexiglas box filled with 0.012 inch diameter silica glass 
beads as shown in Figure 4. The box was equipped with a 
movable side wall which permitted the inside dimensions of 
the bed to be varied from 6.25 inches wide by 6 inches deep 
by 18 inches high to 24 inches wide by 6 inches deep by 18 
inches high. The movable wall was equipped with felt gaskets 
which prevented the escape of beads through the gap between 
the movable side wall and the fixed front and rear walls. 

0.5 inch thick clear plexiglas (grade GM) was used through- 
out the construction and all components were assembled using 
plastic laminating cement with wood screw reinforcements. 

The floor of the apparatus also served as the distributor 
plate whose function was to evenly distribute the flow of 
incoming air over the entire cross section of the bed. The 
distributor plate was constructed from a plexiglas sheet 
with 0.125 inch diameter holes drilled 0.25 inches on center. 
Secured over this plate was a 140 mesh stainless steel wire 
screen to prevent the beads from falling down into the holes 
and to provide an even finer distribution of air. 

Located beneath the distributor plate was an air chamber 
or plenum which served to ensure that the air pressure on 
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Bed Thermocouple Probe Taps Air Outlet 
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Air Distributor Tube 



the underside of the distributor plate was the same at any 
point. This was accomplished by first directing the incoming 
air through a 2.5 inch inside diameter plexiglas distribution 
tube centered lengthwise within the 24 inch by 6 inch by 6 
inch plenum. This tube had 0.5 inch diameter holes drilled 
on 1 inch centers on its underside and served to evenly direct 
the incoming air downward over the entire length of the plenum 
In order to prevent the escape of air when the bed was con- 
figured to a width less than the 24 inch width of the distri- 
butor plate, the movable side wall was equipped with a 
horizontal baffle plate which laid flat over the top of the 
unused portion of the distributor plate. Once the desired 
position of the movable wall was set, the baffle plate was 
securely clamped to the distributor plate and the edges were 
sealed with silicon caulking. After passing through the 
distributor plate and finally the bed itself, the fluidizing 
air escaped from the apparatus via one or more of three 
plexiglas exhaust tubes positioned on the top of the apparatus 
These tubes were quipped with wire mesh covers to prevent the 
loss of glass beads due to carryover. 

In order to measure the air pressure at various depths 
the apparatus was equipped with static pressure taps posi- 
tioned along the rear wall with 2 inch vertical spacing 
between centers. Additionally taps were installed just 
below the distributor plate and immediately above it. These 
taps were connected by flexible plastic tubing to a manifold 
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which permitted selective connection of the taps to a 
Meriam Model 33KA35 30 inch water manometer as shown in 
Figure 5. 

The fluidizing apparatus was supplied with air from a 
Spencer axial flow turbo compressor rated at 500 SCFM. 

The fluidizing air was not filtered or dehydrated prior 
to entry. The volumetric flow rate of air entering the 
apparatus was measured by a calibrated Fisher and Porter 
flowmeter rated at 34 CFM. 

B. DESCRIPTION OF HEAT TRANSFER APPARATUS 

The test apparatus was equipped with two instrumented 
flat plate heater assemblies each measuring 5 inches wide 
by 10 inches high. These assemblies were mounted flush with 
the inside surface of both the movable side wall and the 
fixed side wall opposite it. Each heater assembly was 
powered by a Watlow strip heater rated at 250 watts. The 
heater consisted of a grid of nickel alloy resistance wire 
embedded in a thin" sheet of rubberized supporting material 
and backed by a 0.5 inch thick sheet of silicone sponge 
rubber insulation. This insulation was in turn bonded to 
a 0.5 inch thick clear plexiglas mounting block. The face 
of the assembly consisted of a 0.25 inch thick copper plate 
mounted flush with the heating surface of the strip heater. 
In order to avoid the use of screws or other conventional 
attachment methods which may have altered or impeded the 
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heat flow from the strip heater to the copper plate, the 
two were held together using a strongback arrangement as 
shown in Figures 6 and 7. This mounting method ensured 
that the strip heater was held firmly against the copper 
plate, thus helping to minimize thermal contact resistance 
between the two. To further reduce this resistance, a ther- 
mally conductive paste was applied to the heating surface 
of the strip heater prior to positioning the copper plate 
over it. Power was supplied to the left hand heater by a 
Lambda regulated DC power supply model LK345A, while the 
right hand heater was supplied by a Hewlett Packard DC 
power supply Model 6296A. Both power supplies were rated 
at 60 Volts and 2 Amps and were each equipped with an in- 
dividual voltmeter and ammeter. 

All thermocouples used in this apparatus were 10 gauge, 
ANSI Type T, copper - cons tantan , exposed junction type. Each 
heater assembly was equipped with twenty thermocouples 
arranged in five rows of four thermocouples each as shown 
in Figures 8 and 9. These thermocouples were mounted to the 
copper plate in such a manner so as to be least disruptive 
of the heat distribution patterns within the plate. The 
thermocouple junctions were soldered into small shallow holes 
drilled into the back of the plates so that the top of the 
junction was flush with the plate surface. The thermocouple 
leads were then routed straight back through holes drilled 
into the strip heater and its mounting block. Figure 10 
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Figure 6. Heater Assembly Mounting Diagram; Side View 
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Figure 7. Heater Assembly Mounting Diagram; Rear View 
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Figure 8. Right Hand Heater Thermocouple Arrangement 
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Figure 9. Left Hand Heater Thermocouple Arrangement 
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Figure 10. 
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shows the thermocouple mounting procedure. X-ray photographs 
were used to determine the exact location of the resistance 
wires embedded within the heaters so as to avoid striking 
them while drilling the holes. This thermocouple arrangement 
made it possible to read the plate inside surface temperature 
immediately adjacent to the surface of the strip heater while 
causing minimum disruption of the heat flux patterns within 
the plate itself. A determination of the plate outside 
surface temperature as seen by the bed was then accomplished 
by means of simple calculations as shown in section IIIE. 

The temperatures within the bed itself were monitored by 
15 retractable, 24 inch, stainless steel sheathed, exposed 
junction, Omega Engineering Type T thermocouple probes. The 
position of these probes was held fixed at a bed depth of 
7 inches above the distributor plate for all data runs. The 
number of bed probes actually utilized depended on the geo- 
metric configuration being studied, as shown in Figure 11. 

Incoming air temperature was monitored by 3 retractable, 

4 inch, stainless steel, sheathed thermocouple probes posi- 
tioned within the air plenum as shown in Figure 12. Out- 
going air temperature was monitored by 3 similar probes 
mounted one each in the center of the 3 air exhaust tubes, 
as shown in Figure 13. Because of the positioning of these 
tubes with respect to the movable side wall only one thermo- 
couple was used to read the air outlet temperature for all 
data runs . 
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NOTE: Arrangement is as viewed from top front. 
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Figure 12. Air Inlet Thermocouple Probe Arrangement 
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Figure 13. Air Outlet Thermocouple Probe Arrangement 
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Figure 14 



Figure 1-S 




Right Hand Heater Mounting Block Thermocouple 
Arrangment 




Left Hand Heater Mounting Block Thermocouple 
Arrangment 
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Figure 16. Right Hand Heater Sidewall Thermocouple Arrangement; 
Rear View 




Figure 17. Left Hand Heater Sidewall Thermocouple Arrangement; 
Rear View 



42 



Three insulated thermocouples were mounted to the outside 
surface of each of the plexiglas heater mounting blocks, as 
shown in Figures 14 and 15, in order to determine the mounting 
block outer surface temperatures. These temperatures were 
used to estimate the amount of heat lost through the back of 
the heater assemblies as explained in section IIIE. Thermo- 
couples were also mounted on the inside and outside surfaces 
of the plexiglas sidewalls, as shown in Figures 16 and 17, 
in order to estimate the losses to the atmosphere and to the 
bed itself from the thin plexiglas wall strips adjacent to 
the edges of the heater assemblies. A final set of thermo- 
couples was mounted to the inside and outside surfaces of 
the container rear wall, as shown in Figure 18, in order to 
obtain an estimate of losses through the rear wall. Because 
of symmetry it was assumed that the losses through the front 
wall were identical to the rear wall losses. 

All thermocouples used in the test apparatus were sampled 
by an HP 3497A data acquisition system controlled by an HP 85 
desk top computer. This system was capable of sampling and 
storing the output of seventy five thermocouples within seven 
seconds. Because of the steady state conditions in existance 
during each of data runs, the seven second sampling time was 
considered brief enough to treat the thermocouple readings 
as through they had all been taken simultaneously. 

Once all of the thermocouples had been sampled, the HP 85 
retrieved the stored thermocouple voltages and converted them 



to temperature values using equations previously developed. 
These temperature values, along with other stored parameters, - 
were then recorded on magnetic disks using a Hewlett Packard 
Model HP 82901M dual flexible disk drive. This recorded 
data was later processed and permanently displayed using an 
HP 82905B printer and an HP 7225A plotter. 

Computer programs used to control data acquisition and 
data processing are shown in Appendixes J and K, respectively. 
The arrangement of the heat transfer data collection and pro- 
cessing system is shown in Figure 19 while Figures 20 and 21 
show photographic views of the entire experimental system. 

C. EXPERIMENTAL APPARATUS CALIBRATION 

The inlet air flowmeter was calibrated using an ASME 
Herschel - type venturi with an inlet diameter of 4.26 inches 
and a throat diameter of 2.13 inches. Experimental data 
(shown in Appendix B) was collected for nine different data 
points and from this data, volumetric air flow rates were 
calculated using equations shown in Appendix H. A plot of 
air flowrate vs flowmeter reading is shown in Figure 22. A 
least squares curve fit was applied to the data points to 
obtain a linear equation which converts flowmeter readings 
to air flow rates in CFM. This equation was then incorporated 
into the data reduction computer program shown in Appendix K. 
The rotameter calibration data shows that a 100% flowmeter 
reading equates to 34. .0 CFM. 
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Figure 18. Rear Wall Thermocouple Arrangement 
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Figure 19. 
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Photograph of Experimental Apparatus 
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Figure 21. Photograph of Experimental Apparatus 
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ROTOMETER READING (X) 
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Figure 23. Photograph of Calibration Apparatus 






Thermocouples number 1 through 67 were calibrated using 
a Rosemount Commutating Bridge Model 920A and a Rosemount 
Calibration Bath Model 913A using liquid nitrogen and water. 
The HP 3497A and the HP 85 were included in the system during 
calibration, as shown in Figure 23. Data was collected for 
twenty different bath temperatures and a least squares curve 
fit was applied to the data to obtain a fourth order poly- 
nomial function for each of the thermocouples. The constants 
from these functions were then incorporated into the data 
acquisition computer program shown in Appendix J to allow 
automatic conversion of thermocouple output voltage to tem- 
perature, in degrees Fahrenheit. A sample of the calibration 
data in shown in Appendix D. 

Thermocouples number 68 through 79 were not calibrated 
but were instead assigned constants from the standard power 
expansion for Type T thermocouples. [Ref. 5] They were then 
compared in a water bath of known temperature to confirm 
proper operation. 

The water manometer used to determine static pressure 
drop within the bed was regularly rechecked for zero readings 
and was readjusted as necessary. 

D. EXPERIMENTAL PROCEDURE 

Preparation for the taking of heat transfer data was 
accomplished by adjusting the movable side wall to the proper 
position for the geometry to be studied. The baffle plate 
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was then clamped into position and the joints were sealed 
with silicone rubber caulking. This caulking was allowed 
to harden for a minimum of twenty four hours to ensure a 
leak tight seal. Carefully weighed quantities of glass 
beads were then added to the bed through the right hand 
air exhaust tube until the depth of the bed reached 12.5 
inches above the distributor plate. The air supply turbo 
compressor was then lined up and energized. Following this, 
the air inlet control valve was opened and adjusted until 
the desired air flow rate registered on the flowmeter. The 
heater power supplies were then energized and adjusted to 
the desired level. 

The power supply controls were set so that the voltage 
going to each of the heaters was the same. Current into 
the heaters was determined by the resistance of the nickel- 
heating element which was in turn determined by the temper- 
ature of the heaters. Once the voltage was equilized, the 
bed was allowed to come to steady state conditions. Prior 
to' the taking of data, the system was checked for thermal 
equilibrium by monitoring the temperature at specific points 
using a special computer program which sampled thermocouples 
at five-second intervals. Once it was confirmed that the 
heater assemblies and the bed itself were at a steady tem- 
perature, the data acquisition program was initiated. At 
the beginning of the program the computer requested that 
the operator provide certain information regarding the data 
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run to be recorded. The first item of information requested 
was the run number. ThLs number was structured so as to be 
as informative as possible. The first two digits represented 
the bed width in inches. These two digits were followed by 
a letter indicating which of the heaters was energized; "L" 
designated left, "R" stood for right, and "B" designated both. 
This letter was followed by two digits indicating the air 
flowmeter reading. Finally, a single letter was used for 
special purpose designations such as a repeated run without 
a change in configuration or flow rate. The next item of 
information requested by the program was the date and time. 
This was followed by requests for bed width and height, air 
flow rate, ambient temperature (read on a mercury in glass 
thermometer positioned near the test apparatus), voltage 
and current readings from the two power supplies to each of 
the heaters, the number of bed thermocouple probes in use, 
and the amount of bed expansion, if any. Finally, general 
comments concerning observed bed behavior were entered. 

As the push of a button, all of the thermocouples were 
then sampled and the corresponding temperatures were then 
recorded on disc along with the other data entered from the 
keyboard. Once the recording process was completed, the 
system was ready to begin another data run. If the next 
run to be studied was for the same bed configuration, all 
that had to-be done was adjust the air flow rate to the new 
reading and wait for the temperatures to stabilize. In 
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some cases, the power supplied to the heaters was also ad- 
justed in order to maintain as large a temperature differ- 
ential as possible between the incoming air and the heaters. 

In order to avoid excessive softening of the plexiglas 
adjacent to the heaters, the maximum heater temperature was 
limited to approximately one hundred eighty degrees 
Fahrenheit . 

After data was collected over the entire range of flow 
rates for a given geometry, the caulking was removed and the 
movable wall was repositioned for the next configuration. 

The entire data collection process was then repeated. In 
all, four different geometric configurations were studied. 

Bed widths of 6.25 inches, 8 inches, 10 inches, and 12 inches 
were used. In addition, data runs were taken in the 6.25 
inch configuration with the bed empty. Data was also col- 
lected with the underside of the distributor to the left of 
the movable wall masked off as a check to see if the baffle 
plate was performing effectively. A total of eighty data runs 
were completed during this study. Experimental heat transfer 
data as well as calculated results are shown in Appendix C. 

Pressure drop data was collected for each bed configura- 
tion by aligning the valve manifold (shown in Figure 5) so 
that the static pressure tap immediately above the distri- 
butor plate was connected to the 30 inch water manometer. 

The fluidizing air flow rate was then adjusted to the desired 
level and the pressure reading at the base of the bed was read 
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on the manometer and recorded. The manifold was then re- 
aligned so that the static pressure tap near the top of 
the bed housing was connected to the manometer. This read- 
ing was then read and recorded. By substracting this read- 
ing from that obtained at the base of the bed, the bed 
pressure drop was determined. Pressure drop data was ob- 
tained in this fashion over the entire flow range for each 
of the bed configurations. During the pressure drop data 
runs, the power to each of the heaters was secured. 

E. EXPERIMENTAL DATA PROCESSING 

All raw data was processed using an HP 85 desk top 
computer. The final results consisted of six pages of 
printed data for each run. The first page was a listing of 
experimental conditions and comments recorded during the 
data acquisition phase of the study. 

The second page was a listing of all thermocouple out- 
puts and their corresponding temperatures. 

The third page showed a temperature profile of the right 5 
hand heater surface as seen from within the bed. In addition, 
an average value of the thermocouple readings in each hori- 
zontal row was calculated and displayed. Thermocouple number 
four was found to be malfunctioning, therefore it's output 
was disregarded. In a separate display, the average row 
temperatures were plotted against bed depth to show how the 
heater temperatures varied from top to bottom. 



The fourth page of the results printout showed the left 
hand heater temperature profile in a format identical to 
that of page three. 

The fifth page showed the temperature profile of the bed 
itself, as measured by the thermocouple probes immersed 
within it. The profile is arranged as seen from above the 
bed . 

The sixth and final page contained the calculated results. 
The first value displayed on this page was the average heater 
temperature (T^) . This value was a simple arithmetic average 
of the five horizontal row temperatures previously calculated. 

T h = { R ( 1 ) + R ( 2 ) + R(3) + R (4 ) + R(5)}/5 (3.1) 

The average bed temperature (T^) was similarly calculated 
using the temperatures displayed on page five of the printout. 

T b = {T (50) + T (51 ) + T ( 5 2 ) + T ( 5 3 ) + T(54) + T(56) 

+ + T(X') }/ (X-50) (3.2) 

X = Highest number bed probe 

The electrical energy into the heater (q ie ) was calculated 
from the voltage and current measurements read off. The heater 
power supplies as shown below: 
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(3.3) 



l ie 



= I x E x 0.0569 ( 



Btu 



Watt -Min 



) 



x 



60 



✓Min', 
1 Hr J 



An estimate of the heat less from the heater to the atmos- 
phere, through the heater mounting block (q was obtained 
as follows: First an average temperature (T ^) was calculated 

for the three thermocouples attached to the outer surface of 
the heater mounting block. 



T b = { T (4 1 ) + T (4 2 ) + T(43)}/3 



(3.4) 



Next it v/ as assumed that the thermal contact resistance 
between the copper plate and the strip heater was negligible, 
therefore the temperature of the strip heater surface was 
equal to . The contact resistance between the strip heater 
and the insulation was also assumed to be negligible. 

The heater insulation and mounting block was then con- 
sidered to be a composite wall as shown in Figure 24. Fourier's 
Law of Conduction was then utilized to calculate the ’heat flux 
through the mounting block. 



^mb 



T h ' T mb 



L. 



(— V + (- 



J mb 



Vi 



k ,A , 
mb mb 



(3.5) 
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Figure 24. Heater Assembly Temperature Profile 
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0 . 5 inches 



= 5 in. x 10 in. = 50 in 



2 



k. = 0.025 



l 






the mounting block to the atmosphere was calculated using 
Newton's Law of Cooling. 



Then an estimate was made of the heat lost to the atmos- 
phere from the hot plexiglas side wall sections adjacent to 
the side edges of the heater (q ga ) • This estimate was based 
on the side wall outer surface temperature and used Newtons 
Law of Cooling and the convection heat transfer coefficient 
previously calculated for the heater mounting block. Since 
the side wall sections and the mounting block were in such 
close proximity it was assumed that the coefficients would 
be in close agreement with one another. 



h 




(3.6) 




(3.7) 



A s = 0.5 in. x 10 in. x x 2 



(one strip on each side of heater) 
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A similar loss estimate was then calculated for the heat 

flux to the atmosphere from the side wall sections adjacent to 

the top edge of the heater (q„ ). 

n t a 



Ota “ h mba x A t x < T (68) ’ T a } 



= 0.5 in. x 6 in. x 



1 Ft' 



144 in‘ 



(3.8) 



Yet another estimate was then calculated for the loss to 
the atmosphere from the side wall section adjacent to the 
bottom edge of the heater (q. ) . 



%a * h mba x A b x «(?!) ' V 



A^.= 0.S in. x 6 in. x 



1 Ft ‘ 



144 in 






(3.9) 



All atmospheric losses were then totaled. 



Oa ■ %,b ♦ q sa * q ta * q ba (3.10) 

The difference between the electrical power input, q g ^, 
and losses from the heater to the atmosphere, q , is the heat 
flux transferred to the bed. In order to identify that portion 
transferred directly from the heated copper plate to the bed, 
an estimate of the heat flux from the heated plastic strip 
surrounding the copper plate to the bed must be made. This 
is done by assuming that same heat transfer coefficient applies 



60 



to both the copper plate and to the plastic strip. The method 
selected was an iterative; wherein it was. first assumed that 
the heat flux passed through the copper plate, from which the 
heat transfer coefficient was calculated. Using this coeffi- 
cient, the heat flux from the heated plastic strip was cal- 
culated which in turn yielded a new heat flux from the copper 
plate. The process was repeated until the heat flux values 
did not change significantly. Thus following this procedure 
the heat flux into the bed from the heater q c ^ was assumed 
to be equal to the incoming electrical power q g( j minus the 
losses from the heater to the a-tmosphere q . 



q cb q ei q a 



(3.11) 



Once the heat flux through the copper plate q c ^ was known, 
the plate outside surface temperature T^ seen by the bed 
could be calculated. 




q cb 
E x 




(3.12) 



L = 0.25 in. x ( i-^ ) 
c v 12 in' 

k c = 232 - (0.032 x (T h - 70)} Btu/Hr- Ft -°F 

1 Ft ^ 

A c = 5 in. x 10 in. x (prTTn 2 '* 
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Once the outside plate surface temperature was known, the 
convection heat transfer coefficient from the plate to the 
bed (h ^3 could be calculated. 



h 



cb 



^cb 



V T P ■ V 



( 3 . 13 ) 



Using the plate to bed heat transfer coefficient (h ^3 
just determined, an estimate was made of the heat flux which 
reached the bed from the hot plexiglas strips adjacent to the 
heater sides ( q s ^ 3 • 



‘•sb ■ h cb x A s x {T(73) ■ V 



( 3 . 14 ) 



Similar estimates were then prepared for the flux entering 
the bed from the hot sidewall section adjacent to the top edge 
of the heater (q as well as the section adjacent to the 
bottom edge of the heater (q^) * 

q tb = h cb x A t x {T ( 69) • T b } (3 - 15) 

% b = h cb X A b X {T(70) - V (3 * 16) 

The heat flux into the bed from all sources other than the 
heater plate itself were then added to obtain a total value 

(C W- 
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^ob ^sb + ^tb + ^bb 



(3.17) 



The original copper plate heat flux value was then cor- 
rected to allow for that amount (q which did not pass 
through the plates. 

Icb ■ %i • % • ‘U (3 - 18) 

It was then necessary to recalculate all of the parameters 
dependent on this value. This process was repeated in an iter- 
ative fashion until the calculated flux through the plate did 
not change s ignf iciant ly . Once the iterative process had been 
completed, the final values were printed out and the entire 
calculation process was repeated for the left hand heater. 

The left hand heater calculations were similar to those used 
for the right hand with the exception of variable names and 
thermocouple numbers. Because of a limitation on the total 
number. of thermocouples which could be controlled by the 
HP 3497A, there were no thermocouples available to measure 
the inside and outside surface temperature of the plexiglas 
sidewall section adjacent to the bottom edge of the left 
hand heater. This prevented the calculation of a loss 
estimate based on a direct measurement. To compensate for 
this, it was assumed that the losses from the sidewall 
sections adjacent to the top and bottom edges of the heater 
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were identical to one another. Once the left hand heater 



computations were completed and the results displayed, cal- 
culations were begun on an energy balance. The objective 
was to determine the total amount of energy entering the bed 
and then compare that with the amount calculated to be leaving 
in order to validate the data and the calculations used. The 
first parameter to be determined was the density of the air 

within the bed p„ . 

a 



Since the air entering a fluidized bed is known to assume 
the bed temperature within a very short distance of entry, 
the air outlet temperature was used as the bed air temperature. 

Using this value of density, the amount of energy carried 
out of t-he bed by the fluidizing air stream (q aQ ) was 
calculated . 



1.325 x P a 



(T +459.69) 
a 

P = 29.92 in. Hg 
a ° 

T = T ( 5 5 ) 
a v 



(3.19) 




(3.20) 



(3.21) 



T • = { T ( 4 7 ) + T (4 8 ) + T (49 ) }/ 3 

3 . 1 

T = T ( 5 5 ) 
ao 
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The heat loss through the front and rear walls (q^ ) was 

then estimated. Only the rear wall was instrumented, however,, 
because of symmetry, it was assumed that the loss through 
the front wall would be identical. 



^f ra 



k fr x A fr x {T(75) ' T ( ?4 ) } 



L 



fr 



(3.22) 



k f = 0.112 Btu/Hr - Ft - °F 

1 Ft 

A ~ = 12 in. x W x (y^ — t — ) x 2 (to account for 

f r 12m. 

both front and rear walls) 

W = Bed Width 

L c = 0 . 5 in . 

t r 



The total heat flux. leaving the bed q Q was then calculated. 
It was assumed that the only flux paths out of the bed were 
via the front and rear walls and via the outgoing air stream. 
The side walls were considered to be at a sufficiently elevated 
temperature as to block the flow of heat outward from the bed. 
The floor of the bed was considered to be free of losses -since 
any heat passing through it would be returned to the bed by 
the fluidizing air stream as it passed through the distributor. 



q = q + q r (3.2 

H o H oa M fra 

The total energy into the bed from. all right hand sources 
was then calculated (q^ r ) . This included both the flux 
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through the plate and the flux through the sidewalls 
adjacent to the heater. 



^br - ^ei ' ^a 



(5.24) 



A similar calculation was performed for the left hand 
sources (q-^) and then the total energy into the bed from 
all sources was calculated (q^) . 



^b ^br + ^bl 



The superficial gas velocity (U) (in Ft/Sec) was calculated 
using the following equation: 



U 




Min*. 

Sec 



A 



d 



(3.26) 



= Distributor Area 



6 in. x W x 



2 

, 1 Ft 
^l44“Tn 



• 2 ) 



The superficial gas mass velocity (G) was then calculated. 



G = 



Lbm m 



Ft 



Hr-Ft' 



= ^ - U (-^-) x Pfl (^) x ( 



Lbnu ,3600 Sec 



) 



Sec 



Ft' 



Hr 



G = U x p x 3600 

3 . 



(3.27) 
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The final calculation of the program was the determination 



of the particle Reynolds Number (Re) 



p V D 

Re = £ 

R S. 



G x D 



3600 x y x g 



y = 3.96 x 10 ^ Lb f- Sec /Ft” 



g c = 32.174 Lbm- Ft / Lbf - Sec ' 

D = 0.012 in. x (— l— 
p 12 m. 



(3.28) 



Following the completion of this calculation the energy balance 
section of page six was printed out and the program was 
terminated. 
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IV. PRESENTATION AND DISCUSSION OF RESULTS 



A. FLOW PATTERNS 

The particle flow patterns observed in the course of this 
study were in general agreement with those described in the 
literature for a bubbling bed. The particles moved upward 
near the centerline of the bed as did most of the bubbles. 

Once the particles reached the surface of the bed they were 
thrown outward from the center and then proceeded to travel 
downward along the bed walls as shown in Figure 25. The 
area of greatest activity was roughly cylindrical in shape 
with the corner areas exhibiting considerably less particle 
motion. Observed particle stagnation regions are shown in 
Figure 26. No indications of channeling were observed visually 
at any time. At high fluidization flow rates in the narrow 
bed configuration, slugging was occasionally observed. This 
was characterized by the appearance of large bubbles and the 
sudden lifting of the bed surface to a level two to three 
inches above its previous position. Once the bubble reached 
the top of the bed, the surface would suddenly collapse and 
resume its former level until the next bubble appeared. 

The particles within the bed were observed to behave in 
a manner much like a viscous fluid. The layer of particles 
directly adjacent to the walls showed little, if any, move- 
ment. Because of their transparent composition it was possible 
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Figure 25. Observed Particle Circulation Flow Paths 




Figure 26. Observed Particle Stagnation Regions 
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to observe that the next layer of particles was moving at a 
slightly higher velocity and that the third layer inboard 
was moving even faster. This movement was not steady but 
tended to occur in surges. It is possible that a static 
charge contributed to the "No slip" condition observed at 
the boundary. After the fluidizing air flow was secured upon 
completion of a test run it was noted that numerous particles 
were clinging to the container walls in the upper portion of 
the bed. In an attempt to eliminate this condition, an anti- 
static film was sprayed on the upper portions of the bed walls. 
For a short period of time, the walls remained free of cling- 
ing particles, however, the effect did not last for more than 
two or three data runs. Following this, no further attempts 
were made to alter this phenomenon. Research by Miller and 
Logwinuk [Ref. 2] suggested that the presence of such a static 
charge could contribute to erratic heat transfer results, how- 
ever for the purpose of this study, any affect was considered 
to be negligible . - 

Closer examination of the bubble flow path showed that 
many bubbles tended to originate from the left side of the 
bed near the bottom of the movable side wall. As a result, 
there was a shift in overall particle activity to the left 
of center and a corresponding alteration of the quiescent 
zone near the bottom of the bed such that the majority of 
stagnant particles were located in the right hand corner. 
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Figure 27 . 



Photograph 



of Off-Center Bubble 
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This phenomenon was particularly noticeable in the wide bed 
configurations. An example of this off-center activity is 
shown in Figure 27. 

Examination of the apparatus revealed that a stream of 
air was being injected into the bed from the area between 
the baffle plate and the distributor. Air from the inlet 
plenum was passing through the unused portion of the dis- 
tributor to the left of the movable wall. Once the air was 
through the distributor, further upward progress was blocked 
by the baffle plate and the sealed joints. There was however, 
a very thin horizontal passage between the baffle and the top 
of the distributor as shown in Figure 28. The air flowed 
along this path until it emptied into the bed itself. Because 
of this "peripheral injection" there was an uneven distribution 
of fluidizing air and a resulting shift in the center of activ- 
ity toward the left hand side. In order to determine the effect 
of this phenomenon an entire set of data was taken with the 
unused portion of the distributor masked-off from the under- 
side so that no air could enter. During these ma'sked data 
runs it was observed that the bubble path was very well 
centered as was the quiescent zone near the bottom of the bed. 
All other data runs were taken with the distributor unmasked 
in order to study the difference in heat transfer performance 
between the two sides. A detailed study of particle cell 
flow patterns for this apparatus was presented by Morgan 
[Ref. 2] in 1981. • 
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Figure 28. Air Baffle Leakage Flow Path 



B. 



PRESSURE DROP DATA 



Measurements were taken of the pressure drop across the 
bed for all flow rates and bed configurations. This data 
can be found in Appendix A. Plots were constructed of LOG 
AP vs. LOG G for each bed configuration in order to study 
the variation in pressure drop with changes in the super- 
ficial mass velocity. These plots are shown in Figures 29 
through 32. An ideal pressure drop curve is shown in Figure 
35. Each of the plots display the same basic characteristics 
of the ideal curve in that the pressure drop increases 
steadily until the point of minimum fluidization (G f) is 
reached. From that point onward, the pressure tends to 
remain relatively constant despite further increases in air 
flow rate. This is because of the fact that once the bed is 
fluidized, any further increase in the flow rate imparts an 
increased lift to the particles because of the viscous drag 
force exerted on them by the moving air stream. As the lift 
increases, expansion of the bed occurs. and voidage (e) in- 
creases. As the voidage increases, the interstitial spacing 
between the particles also increases. This results in a 
reduction of the interstitial gas velocity, which in turn, 
lowers the drag exerted on the particles until the gravita- 
tional force and the lift are once again in balance. As a 
result of the increased voidage, the pressure drop across 
the bed remains constant or may even decrease. 
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Examination of the pressure drop curve for the 6.25 inch 
configuration (Fig. 29) shows a deviation from ideal behavior 
in that the pressure drops sharply after reaching the point 
of minimum fluidization. A pressure drop of this sort is 
normally associated with channeling. As indicated previously, 
no visual indications of channeling were detected, however 
it is possible that "intermittent channeling” was taking 
place. This type of channeling would be difficult to differ- 
entiate from rapid bubbling because the channel does not 
extend through the entire depth of the bed as it does in 
"through channeling." An alternate explanation for the drop 
in pressure is the effect of particle interlocking which some- 
times manifests itself during the initial fluidization of a 
bed. It is noted that subsequent data runs did not demonstrate 
such behavior. 

Based on the curve shown in Figure 29 the point of minimum 

fluidization for the 6.25 inch configuration occurred at G=159.6 
2 

Lbm/Hr-Ft (26% flow). It was- was noted the visual observation 
of fluidization onset occurred at the same air flow rate as 
was predicted by this curve. 

The pressure drop curve for the 8 inch configuration 
(Fig. 30) is in very close agreement with the ideal curve 
with no indications of channeling or slugging. Based on this 
curve the point of minimum fluidization occurred at G=132.6 
Lbm/Hr-Ft (28% flow). Fluidization was visually observed 

t 

to commence at 137 Lbm/Hr-Ft^ (29% flow). 
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LOG DELTA P ve LOG G 
6. 25-INCH CONFIGURATION 
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Figure 29. 



Pressure Drop vs Superficial Mass Velocity; 
6.25 in . Bed 
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LOG DELTA P ve LOG G 
0-INCH CONFIGURATION 
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Figure 30. Pressure Drop vs Superficial Mass Velocity; 
8.0 in . Bed 
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LOG DELTA P vs LOG G 
10-INCH CONFIGURATION 
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Figure 31. 



Pressure Drop vs Superficial Mas's Velocity; 
10.0 in . Bed 
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LOG DELTA P vs LOG G 
12- INCH CONFIGURATION 
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Figure 32. Pressure Drop vs Superficial Mass Velocity; 
12.0 in . Bed 
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Figure 33. Pressure Drop vs Superficial Mass Velocity;” 
Ideal 
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The curve for the 10 inch configuration (Fig. 31) was 
smoother than the previous configurations and matched the 
ideal curve closely except for a gradual decrease in pres- 
sure following fluidization onset. The reason for this 
decrease was not readily apparent. The curve shows that 
fluidization commenced at G=148 Lbm/Hr-Ft^ (401 flow). 

The observed starting point of fluidization was also at 
401 flow. 

The 12 inch configuration curve (Fig. 32) was also very 
smooth and similar to the ideal. Like the 10 inch curve, 
a slight decrease in differential pressure was noted. The 
starting point of fluidization (G £) indicated by the curve 
was G=146.5 Lbm/Hr-Ft" (48% flow). Fluidization was visually 
observed to commence at the same point. 

A combined graph showing the pressure drop curves for all 
four configurations is shown in Figure 34. With the exception 
of the 8 inch configuration, they show a slight gradual de- 
crease in the minimum fluidization superficial mass velocity 

(G r ) as the bed width is increased. A similar trend was 
m± 

observed by Morgan [Ref. 2] who attributed this to a reduction 
in the effect of the wall frictional drag as a fraction of 
the overall drag experienced by the bed. 

The average bed pressure drop at minimum fluidization 
(P m £) for all configurations was approximately 16.0 in. H^O 
or 0.577 psi. • Research by Trivedi et al [Ref. 6] indicated 
that the pressure drop across a fluidized bed should equal 
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Figure 34. Pressure Drop vs Superficial Mass Velocity; 
All Configurations 
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the weight of the bed material per unit cross sectional area. 
The settled bed density for the glass beads used in this ex- 
periment was 96 Lbm/Ft 0 which resulted in a weight - to -area 

2 

ratio of 0.681 Lbm/In . Thus, the pressure drop observed 
was slightly less than the theoretical value. Recalculation 
of the weight-to-area ratio based on the expanded bed density 
yielded essentially the same results. 

C. TEMPERATURE DISTRIBUTION 

The temperature distribution across each of the heater 
surfaces was examined for all data runs and two distinct 
trends were found to be evident. First; the four thermo- 
couple readings in any horizontal row were found to be in 
very close agreement with one another in all cases as shown 
in Appendix E. This even distribution of temperature hori- 
zontally across the plate suggested that the local rate of 
heat transfer from the plate was the same for any given point 
along the plate surface at a given bed depth. This conclusion 
was based on the assumption that the heat flux into the plate 
was also horizontally uniform. The second trend observed was 
that there was a variation in heater temperature from the top 
of the plate to the bottom and that this "temperature profile" 
displayed a distinct and predictable change in shape as the 
air fluidization flow rate was increased. Figures 35 through 
44 illustrate these trends. At low air flow rates the point 
of greatest temperature was always at the top of the plate. 
Moving downward, the temperature then gradually dropped in 
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a non-linear but smooth fashion to a level approximately three 
to five degrees less than that observed at the top of the 
plate. Examination of the plots of temperature vs bed depth 
for the low flow rate runs showed all of the curves to have 
the same shape with each curve displaying a positive slope. 

It was noted that as the flow rates increased, the slope of 
these curves became more positive until they assumed a ver- 
tical orientation and then reversed their slope to become 
negative. It was further noted that this slope reversal 
occurred at or near the point of fluidization onset. The 
change in slope of these curves meant that the peak tempera- 
ture within the plate was moving downward from the top of 
the plate toward the bottom as the fluidization air flow rate 
was increased. A possible explanation of such behavior lies 
in the fact that prior to fluidization onset, when the par- 
ticles were not in motion, the heater was primarily cooled 
by an upward flow of air from the distributor. Because of 
this upward flow path, the temperature of the heater tended 
to increase with an upward movement along the plate. Once 
the bed was fluidized, however, a significant amount of 
cooling was then provided by the moving particles sweeping 
along the plate surface. Since the particles were moving 
downward along the wall a thermal gradient was introduced 
which was opposite to that caused by the upward flow of air. 
The peak temperature within the plate would consequently be 
observed at the point where these two gradients met. If the 
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cooling due to the particles was greater than that provided 
by the upward air flow, the peak temperature would be found 
toward the bottom of the plate. Conversely, if the air flow 
cooling predominated, the peak temperature would be seen in 
the upper section of the plate. An increase in the fluidiza- 
tion air flow rate would provide more upward flowing cooling 
air while simultaneously increasing the downward moving 
particle activity. 

Experimental results showed that an increase in air flow 
rate caused the peak temperature . to move downward. This sug- 
gested that, as the flow rate was increased, the increase in 
particle cooling predominated over the increase in air cooling. 
The reason for this predominance may have been that the par- 
ticles had a much higher specific heat capacity than the air 
and a lower mean velocity. Consequently, they were able to 
remove a larger quantity of heat. Another possible cause for 
the particle cooling predominance as flow rates were increased 
was that the motion of the cooling air may have been altered 
once fluidization commenced so that more of the upward moving 
gas was concentrated in the center of the bed. Although there 
was no objective evidence that this was occurring, it may have 
been that any increase in flow rate was merely diverted to 
the center of the bed and consequently did not increase the 
amount of cooling air available near the wall. 

Another trend which was observed while studying the heater 
temperature distributions was that as the bed width increased. 
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the shift in the orientation of the temperature curve became 
less pronounced. In the 6 inch configuration the peak 
temperature point moved from the top of the bed downward to 
the very bottom as fluidization activity increased. The 
maximum temperature differential observed between the top 
and the bottom of the plate was approximately 17 degrees. 

In the 8 inch configuration, the peak temperature point 
also transited down the entire length of the plate, however 
it did so at a slower rate than in the 6 inch configuration. 
The maximum top- to-bottom differential was essentially the 
same as that seen for the 6 inch bed. With the bed con- 
figured to the 10 inch width, the peak temperature did not 
move down the entire length of the plate, but instead stopped 
at a level approximately 3 inches above the bottom edge. The 
shift in temperature orientation was even more gradual than 
in the 8 inch configuration. The maximum observed tempera- 
ture differential was approximately 11 degrees. While in 
the 12 inch configuration, very little change was observed 
in the temperature distribution as the flow rate was in- 
creased. The peak temperature shifted from the top of the 
plate to just above center, however the maximum temperature 
differential of 6 degrees was significantly less than that 
observed in the other configurations. 

In general, the right hand heater temperature distribu- 
tion proved to be more sensitive to changes in flow rate 
than did the left hand heater. 
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A characteristic unique to the left hand heater was that 
immediately prior to commencement of the downward transit of 
the peak temperature point, the temperature differential 
between the top and the bottom of the plate suddenly increased 
to approximately three times its previous value. As the flow 
rate was further increased, this differential abruptly dropped 
back down to its former level. This trend is shown in Figures 
40 through 43. It is believed that this fluctuation could 
be related to a possible reorientation of air flow patterns 
as fluidization commences. This phenomenon was not evident 
in the right hand heater temperatures. 

Examination of the data in Appendix C shows that for each 
flow rate and bed configuration, the average temperature of 
the left hand heater was lower than for the right hand heater. 

In many cases, the difference in cooling of the plates 
was quite significant. During the data runs with the distri- 
butor masked-off, it was observed that the left and right 
hand heaters were being cooled to approximately the same 
degree, since the average heater temperatures were in close 
agreement with one another. 

The temperature profile of the bed itself, as measured 
by the immersed probes, showed a non-uniform distribution 
prior to the commencement of fluidization for all configura- 
tions. The temperatures were highest near the heater sur- 
faces where particle heating was taking place by conduction. 
Because of the relatively low rate of heat transfer the 
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Figure 35. 



Right Hand Heater Temperature Profiles;' 6.25 in. 
Configuration 
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Figure 36. 



Right Hand Heater Temperature Profiles; 6.25 in. 
Masked Configuration 
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Figure 37. 



Right Rand Heater Temperature Profiles; 8.0 in. 
Configuration 
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Figure 38. . Right Hand Heater Temperature Profiles 
Configuration 



10.0 in . 



91 



Vertical Position 




Temp. 







represent flowmeter readings as a 
percentage of maximum rotometer capacity. 



Figure 39. 



Right Hand Heater Temperature Profiles; 12.0 in. 
Configuration 
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NOTE: Numerical values represent flowmeter readings as a 

percentage of maximum rotometer capacity. 



Figure 40. 



Left Hand Heater Temperature Profiles; 6.25 in. 
Configuration 
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NOTE: Numerical values represent flowmeter readings as a 

percentage of maximum rotometer capacity. 
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Left Hand Heater Temperature Profiles; 6.25 in. 
Masked Configuration 
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NOTE: Numerical values represent flowmeter readings as a 
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Figure 42. 



Left Hand Heater Temperature Profiles; 8.0 in. 
Conf igurat ion 
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Figure 43. 



Left Hand Heater Temperature Profiles; 10.0 in. 
Configuration 
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Figure 44. 



Left Hand Heater Temperature Profiles; 12.0 in. 
Configuration 
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temperature tended to decrease toward the center. Immediately 
upon commencement of fluidization, however, the temperature 
distribution within the bed became extremely uniform and 
remained so as fluidization flow rates increased. This was 
directly attributable to the thorough mixing induced by the 
fluidized particle motion. In most cases, the temperature 
of the outgoing airstream was within two or three of the 
temperature of the bed itself. 

D. HEAT TRANSFER COEFFICIENTS 

In order to establish a performance baseline, data was 
collected for a bed from which all of the particles had been 
removed and an average heat transfer coefficient of 3.8 Btu/ 

2 o 

Hr-Ft - F was calculated. All remaining runs were accom- 
plished with the apparatus filled to a settled bed depth of 
approximately 12 inches. The resulting heat transfer data 
for all runs is tabulated in Appendix C and plots of heat 
transfer coefficient vs. superficial mass velocity are 
shown in Figures 45 through S3. Each of these plots demon- 
strates the same basic trend in that initially, during low 

gas flow rate conditions, the heat transfer coefficient was 

2 o 

correspondingly low (approximately 4.5 Btu/Hr-Ft - F) and 
remained low until a flow rate approaching the minimum 
fluidization superficial mass velocity ( G £ ) was attained. 

From this point onward, the heat transfer coefficient in- 
creased in a roughly linear fashion as the gas* flow rate was 
increased . 
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8. 25- INCH CONFIGURATION 



Figure 4 
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6. 25- INCH CONFIGURATION-MASKED DISTRIBUTOR 
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Figure 46. Heat Transfer Coefficient vs Superficial Mass 
Velocity; 6. 25 in. Masked Configuration 
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Configuration 
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Figure 48. Heat Transfer Coefficient vs Superficial Mass 
Velocity; Left Hand Heater, 6. 25 in. 
Configuration 
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Figure 49. Heat Transfer Coefficient vs Superficial Mass 
Velocity; 8.0 in. Configuration 
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Figure 50. Heat Transfer Coefficient vs Superficial Mass 
Velocity; 10.0 in. Configuration 
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Figure 51. Heat Transfer Coefficient vs Superficial Mass 
Velocity; 12.0 in. Configuration 
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Figure 52. Heat Transfer Coefficient vs Superficial Mass 

Velocity; Right Hand Heater, All Configurations 
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A second significant trend which was observed was the 
fact that the heat transfer coefficient for the left hand 
heater was significantly higher than for the right hand 
heater, particularly as the bed width was increased. In 
the wide bed configurations the heat transfer coefficient 
for the right hand heater showed very little improvement 
in response to an increase in flow rate. With a high 
height- to-width ratio (narrow bed) , the heat transfer 
coefficient for the left hand heater was roughly twice 
the value for the right hand heater at high gas flow rates. 
With a low height-to-width ratio (wide bed) , the left hand 
heat transfer coefficient exceeded that for the right by 
a factor of 15 to 1 . 

The exception to this trend occurred during the data 
runs accomplished with the distributor masked-off. Under 
this condition, the heat transfer coefficients for the left 
hand and right hand heaters were in much closer agreement 
with one another than in all other cases. The effect of 
masking the distributor was to increase the heat transfer 
coefficient for the right hand heater by approximately 10 % 
while simultaneously reducing the heat transfer coefficient 
for the left hand heater by approximately 251. It would be 
expected that, with the distributor masked-off, the heat 
transfer coefficients would be the same for both the left 
hand and right hand heaters. Although there is a much 
closer match than in the unmasked configuration, it is 
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observed that the left hand heater still achieved higher rates 
of heat transfer. There are several possible explanations 
for this difference. The first of these is that the unused 
portion of the distributor may not have been completely sealed 
by the masking and as a result, some peripheral injection 
could have still have been taking place. Secondly, it was 
noted that there was a small amount of gas leakage past the 
movable side wall gaskets which may have altered the gas and 
particle flow patterns along the left wall. This, in turn, 
may have altered the heat transfer coefficients. Thirdly, 
it is possible that the movable side wall may not have been 
perfectly vertical. Although care was taken to accurately 
position and secure the side wall, deviations from a true 
vertical orientation of as much as 1° are considered possible. 
Research conducted in the Soviet Union by Filippovskii and 
Baskakov [Ref. 7] has indicated that plate angle has a con- 
siderable influence on heat transfer coefficients. Because 
of the small angle involved in this case, the contribution 
of this effect is not considered to be the principle cause 
of the difference in performance between the two plates. 

A final effect which may have made a minor contribution to 
this difference is the fact that the losses from the heater 
backing to the atmosphere were lower for the left hand heater 
than for the right. As a result, more of the electrical 
energy going into the heater was directed through the copper 
plate and on to the bed, thus resulting in a minor improvement 
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in the heat transfer coefficient. Because of the configura- 
tion of the apparatus, the left hand heater mounting block 
was enclosed within the unused portions of the apparatus 
walls. Because of this, there was a stagnation of the air 
in the vicinity of the mounting block as well as an elevation 
of the local ambient temperature. These effects combined to 
reduce the amount of heat lost through the mounting block. 
Because of the specialized conditions which were found to 
exist near the left hand wall, the performance of the right 
hand heater shall be considered as being representative of 
a typical flat plate. 

The magnitude of the heat transfer coefficients calculated 
as a result of this study were within the range of values 
found in the literature, however it should be noted that 
there is a wide variation in the coefficients found by 
other researchers due to the effect of differing bed sizes 
and geometries, particle size and composition, fluidizing 
gas properties, superficial mass velocities, and bed operat- 
ing temperatures. 

There is also a wide variety of empirical correlations 
available for the prediction of heat transfer coefficients, 
however each of these serves, for the most part, only to 
describe the particular test apparatus and conditions from 
which they were derived. For this reason, no attempt was 
made to fit the data to the correlations currently available 
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nor was an attempt made to develop a new correlation to 
describe the performance of the apparatus used in this 
experiment . 

A review of the literature shows that numerous researchers 
have observed heat transfer curves similar to that shown in 
Figure 54. This typical curve displays a peak heat transfer 
coefficient (h max ) at some optimum superficial mass velocity 
(G 0 p t )* Flow rates above this optimum value resulted in a 
reduction in the heat transfer coefficient. This has been 
attributed to the fact that as the voidage is increased, 
there is a point at which the beneficial contribution of 
increased particle activity is offset by the negative effect 
of reducing the solids concentration within the emulsion. 

In this study such an optimum flow rate was not achieved. 

It is believed that, due to the design of the fluidization 
apparatus, the range of gas flow rates was restricted by 
the relatively shallow "freeboard" region above the level 
of the settled bed. When the bed was undergoing active 
fluidization, particles within the "splash zone" reached 
the top of the container and blocked the air exit ports. 

This was believed to be taking place at gas flow rates 
below the optimum level. It is possible that a reduction 
in the settled bed depth and a consequent increase in the 
height of the freeboard region would permit higher super- 
ficial mass velocities and result in higher levels of heat 
transfer . 
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Figure 54. 



Heat Transfer Coefficient vs Superficial 
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The difference in the level of heat transfer between 
the left and right hand walls is clearly related to the 
high level of cooling observed at the left hand wall. 

This, in turn, is related to the large amount of particle 
activity observed to the left of center. The exact mech- 
anism of heat transfer which takes place at the wall cannot 
be stated with certainty. Although there are many theories 
as to the mechanism of heat transfer within a fluidized 
bed, the bulk of the models postulated to date fall into 
one of two broad categories. The first of these models 
proposes that the heat transfer is primarily through a gas 
film adjacent to the heat transfer surface. The heat is 
transferred by steady state conduction through the film to 
the downward moving particles. This film is kept thin by 
the scrubbing action of the solid particles which come in 
contact with the surface. This reduction in the boundary 
layer thickness is thought to be the key to achieving high 
rates of heat transfer. The second model theorizes that 
heat transfer is primarily by unsteady conduction to "packets" 
which are made up of several solid particles grouped together 
for a short period of time. These packets are thought to 
be in contact with the heat transfer surface for only brief 
periods of time before they move away from the surface and 
are replaced or "replenished" by fresh packets. As a result, 
there is always a high temperature differential between the 
surface and the packets and this, in turn, yields high levels 
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of heat transfer. A significant factor shared by both 
the packet model and the film conduction model is that 
they both rely on the free motion of particles to effect 
high rates of heat transfer. It is believed that in this 
apparatus a high level of particle activity near the wall 
promotes the frequent exchange of particles from within the 
bed core where it is comparatively cool. As the frequency 
of particle or packet replenishment increases, the tempera- 
ture differential between the wall and the particles also 
increases. This, in turn, increases the overall heat transfer 
Conversely, when there is comparatively little particle activ- 
ity near the wall, the replenishment rate is low and the same 
particles remain in contact with the hot wall for longer 
periods of time as they move downward. Because of the long 
residence time, the temperature differential is diminished 
and the heat transfer coefficient is reduced. In this study, 
the level of particle activity near the right hand wall was 
comparatively low and so was the heat transfer coefficient. 

A' study of the temperature profiles for the right hand plate 
shows a steady increase in plate temperature with a downward 
change in vertical position. This profile is consistent 
with that expected for long particle residence times since 
the downward moving particles would experience a steady in- 
crease in temperature as they approached the bottom of the 
plate. It is noted that when the configuration of the bed 
is tall and narrow, the bubble flow path is very close to 
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the bed walls. This causes a high degree of particle activity 
near the wall and a corresponding increase in the particle 
replenishment rate. As the bed height - to -width ratio is 
decreased by increasing the bed width, the bubble flow path 
moves further away from the wall. This causes a reduction 
in the amount of activity near the wall and a corresponding 
decrease in the replenishment rate. This scenario is con- 
sistent with the observed reduction in heat transfer coef- 
ficient as the bed width was increased. Conversely, the 
amount of activity seen by the left hand plate was consis- 
tently high as a result of the peripheral injection effect 
described previously. The degree of activity and consequently 
the heat transfer coefficient, remained unaffected by an in- 
crease in bed width. In addition, the heater temperature 
profile showed no sharp transition from top to bottom, thus 
indicating that the particle residence time was short. 

Based on the findings of this study, it may be concluded 
that peripheral injection is an effective means of promoting 
particle activity or turbulence along the wall of the con- 
tainer. The idea of turbulence promoters in fluidized beds 
is not a new one. Early research by the Babcock and Wilcox 
Company [Ref. 8] included the use of spirally shaped strips 
to create a rotational motion which forced the solids toward 
the walls by centrifugal force. The result was a significant 
increase in heat transfer. Recently Sokolov et al [Ref. 9] 
in the Soviet Union studied the use of angled louvers 
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positioned within the bed in order to redirect the motion of 
the particles outward toward the walls, thus destroying the 
stagnant layer which was insulating the surface. The result 
was a marked improvement in heat transfer performance. 
Although more study is obviously warranted, the use of pe- 
ripheral injection would appear to offer benefits similar 
to other enhancement techniques while eliminating the need 
for in-bed structures which occupy valuable space and which 
require periodic replacement as a consequence of particle 
errosion. In addition, the amount of turbulence created at 
the wall by the injection technique could easily be adjusted 
to optimum levels by controlling the flow rate of the in- 
jection gas. Such a system could be used to offset the low 
rates of heat transfer at the walls of wide fluidized beds. 

As discussed previously, an energy balance was conducted 
in order to confirm the validity of the numerical results. 
The results of this energy balance were disappointingly 
poor. In most cases, the measured heat flux leaving the 
bed with the outgoing airstream was found to be only half 
of that calculated to be entering the bed via the electrical 
heaters. It is believed that the disagreement was primarily 
due to incorrect measurement of the fluidizing gas flow 
rate. Despite the poor energy balance, the conclusions 
reached as a result of this study are considered to be valid 
since they were b ased on observed trends rather than on the 
magnitude of specific numerical results. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

Based on the experimental results of this study the 
following conclusions may be drawn: 

1. The plate-to-bed heat transfer coefficient for a flat 
vertical plate serving as a containing wall decreases with 
a decrease in the bed height - to-width ratio. 

2. The plate-to-bed heat transfer coefficient increases 
with an increase in superficial mass velocity above the mini- 
mum required to initiate fluidization. 

3. The plate-to-bed heat transfer coefficient remains small 
and relatively constant at superficial mass velocities below 
the minimum required to initiate fluidization. 

4. The plate-to-bed heat transfer coefficient increases 
with an increase in the particle activity or turbulence in 
the vicinity of the wall. 

■» 

5. The particle turbulence near the vertical plate, and 
consequently the heat transfer coefficient, may be enhanced 
by the injection of air from beneath the plate in a direction 
perpendicular to the primary fluidizing air flow path. 

6. The temperature profile of a vertical flat plate may be 
used to provide insights as to the nature of particle activity 
in the vicinity of the plate. 
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7. The particle flow pattern within a rectangular fluidized 
bed is in very close agreement with that described in the 
literature for cylindrical fluidized beds with upward motion 
of the particles taking place along the centerline of the 
container and downward return flow of the particles taking 
place along the periphery. 

8. Slugging or heaving of the bed was more prevalent in the 
narrow bed configuration as opposed to the wide bed configura- 
tion where this phenomenon was not observed. 

9. Fluidized bed particle motion along a non-moving boundary 
matches the behavior of a viscous fluid in that a "no-slip" 
condition is observed. 

10. The bed pressure drop characteristics for a rectangular 
fluidized bed match those found in the literature for a cy- 
lindrical fluidized bed in that the pressure drop across the 
bed increases linearly with an increase in air flow rate 
until the point of fluidization onset. Once fluidization 
has commenced, the pressure drop remains constant. 

11. The temperature of the heater surface at a given depth 
remains constant for any point in a horizontal direction. 

12. The bed temperature will assume a uniform temperature 
distribution once fluidized. 

B. RECOMMENDATIONS FOR FURTHER STUDIES 

In order to better pursue the ojbectives of this study at 
sometime in the future, the following recommendations and/or 
suggestions are made: 
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1. Replace the three circular air exhaust tubes on the top 
of the fluidization apparatus with a rectangular exhaust 
port which runs the entire twenty-four inch length of the 
bed . This would provide a non- rest rictive exit path for 
the fluidizing air which would minimize the amount of swirl- 
ing and mixing in the region above the bed. 

2. Repair severed wire on thermocouple #15. 

3. Replace malfunctioning thermocouples #4 and #40. 

4. Install a minimum of four new thermocouples to monitor 
air exit temperature. At present only one thermocouple can 

f 

be used for this purpose (T(55)) for bed configurations of 
12 inches or less. 

5. Provide a more positive means of securing the baffle 
plate to the distributor plate so as to prevent leak- by into 
the bed. 

6. Install a bypass valve around the existing air inlet con- 
trol valve to permit fine adjustment of fluidization flow rates. 

7. Install two additional thermocouples in the left hand 
movable wall i/4 inch below the bottom of the heater in 
order to obtain a more accurate estimate of wall losses. 

8. Install two additional thermocouples in the front wall 
of the bed so as to obtain a more accurate estimate of wall 
losses . 

9. Substitute the existing power supplies with units, which 
are capable of providing additional current and which are 
equipped with controls to vary power as well as voltage. 
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This would permit equalization of input power to the two 



heaters and allow the desired temperature differential 
between the plate and the bed to be maintained. 

10. Install a movable wall in the air plenum beneath the 
distributor plate so that the flow path followed by the in- 
coming air is identical for the left and right sides. 

11. Prior to the commencement of future data runs, apply 
caulking material to the movable side wall so as to prevent 
leakage of fluidizing air past the gaskets. 

12. Modify the fluidization apparatus so that a controlled 
amount of injection air may be applied to each of the heater 
plates. Such a modification would permit additional study 
of the effects of the peripheral injection phenomenon. 



C. CLOSING REMARKS 

Gas fluidization is a fascinating process which promises 
to provide a multitude of beneficial applications. Full 
enjoyment of the potential which exists depends on a more 
complete understanding of the many complex phenomenon 
associated with fluidization. It is hoped that this thesis 
has, in some small way, contributed to that understanding 
or helped to stimulate additional thought and research in 
this exciting field of study. 
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Table I. Bed Pressure/Gas Mass Flow Data (6.25 in. Bed) 
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Table IV. Bed Pressure/Gas Mass Flow Data (12.0 in. Bed) 
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APPENDIX B. 



EXPERIMENTAL ROTOMETER CALIBRATION DATA 



♦♦ROTOMETER CALIBRATION** 
ASME HERSCHEL-TYPE VENTURI 
4.26/2. 13 IN DIA 

20 PERCENT READING 



♦♦EXPERIMENTAL DATA** 



ATMOSPHERIC PRESSURE = 29.76910 IN Hg 
VENTURI INLET TEMPERATURE = 72.0 DEG F 
VENTURI INLET STATIC GAGE PRESSURE = 10.70 IN H20 
VENTURI DIFFERENTIAL PRESSURE = .010 IN RED GIL 
SPECIFIC GRAVITY OF RED OIL = .334 



♦♦CALCULATED VALUES** 



ATMOSPHERIC PRESSURE = 

VENTURI INLET GAGE PRESSURE = 
VENTURI INLET ABSOLUTE PRESSURE = 
VENTURI DIFFERENTIAL PRESSURE = 
VENTURI INLET AIR DENSITY = 
VENTURI PRESSURE DROP RATIO (X) = 
EXPANSION FACTOR (Y) = 

ASSUMED REYNOLDS NUMBER = 
DISCHARGE COEFFICIENT <C> = 
VELOCITY OF APPROACH FACTOR (E) = 
FLOW COEFFICIENT <K> = 

AIR MASS FLOW RATE = 

AIR VOLUMETRIC FLOW RATE = 

VENTURI INLET REYNOLDS NUMBER = 



14.6259 PSI 
.33654 PSI 
14.23932 PSI 
.00030 PSI 
.07230 LBM/FT3 
. 00002 
1.00 
3000 
.889 
1 . 03280 
.91316 
.01021 LBM/SEC 
3.469 CFM 
3002 
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**ROTOMETER CALIBRATION** 
ASME HERSCHEL-TYPE VENTURI 
4.26/2. 13 IN DIA 

25 PERCENT READING 



** EXPERIMENTAL DATA** 



ATMOSPHERIC PRESSURE = 29.76910 IN Hg 
VENTURI INLET TEMPERATURE = 72.0 DEG F 
VENTURI INLET STATIC GAGE PRESSURE = 12.35 IN H20 
VENTURI DIFFERENTIAL PRESSURE = .011 IN RED OIL 
SPECIFIC GRAVITY OF RED OIL = .834 



** CALCULATED VALUES** 



ATMOSPHERIC PRESSURE = 

VENTURI INLET GAGE PRESSURE = 
VENTURI INLET ABSOLUTE PRESSURE 
VENTURI DIFFERENTIAL PRESSURE = 
VENTURI INLET AIR DENSITY = 
VENTURI PRESSURE DROP RATIO (X) 
EXPANSION FACTOR (Y) = 

ASSUMED REYNOLDS NUMBER - 
DISCHARGE COEFFICIENT (C) = 
VELOCITY OF APPROACH FACTOR (E> 
FLOW COEFFICIENT <K> = 

AIR MASS FLOW RATE = 

AIR VOLUMETRIC FLOW RATE = 
VENTURI INLET REYNOLDS NUMBER = 



14.6259 PSI 
. 446 1 4 PSI 
14. 17972 PSI 
.00033 PSI 
.07200 LBM/FT3 
. 00002 
1.00 
3150 
. 891 



1 . 03280 
.91970 
.01070 LBM/SEC 
8.916 CFM 
3147 
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**ROTOMETER CALIBRATION** 
ASME HERSCHEL-TYPE VENTURI 
4.26/2. 13 IN DIA 

30 PERCENT READING 



**EXPER I MENTAL DATA** 



ATMOSPHERIC PRESSURE = 29.76910 IN Hg 
VENTURI INLET TEMPERATURE = 72.0 DEG F 
VENTURI INLET STATIC GAGE PRESSURE = 14.75 IN H20 
VENTURI DIFFERENTIAL PRESSURE = .015 IN RED OIL 
SPECIFIC GRAVITY OF RED OIL = .334 



* *CALCULATED VALUES** 



ATMOSPHERIC PRESSURE = 

VENTURI INLET GAGE PRESSURE = 
VENTURI INLET ABSOLUTE PRESSURE 
VENTURI DIFFERENTIAL PRESSURE = 
VENTURI INLET AIR' DENSITY = 
VENTURI PRESSURE DROP RATIO (X) 
EXPANSION FACTOR (Y) = 

ASSUMED REYNOLDS NUMBER = 
DISCHARGE COEFFICIENT (C) = 
VELOCITY OF APPROACH FACTOR (E) 
FLOW COEFFICIENT (K) = 

AIR MASS FLOW RATE = 

AIR VOLUMETRIC FLOW RATE = 
VENTURI INLET REYNOLDS NUMBER = 



14.6259 PSI 
.53284 PSI 
14.09302 PSI 
.00045 PSI 
.07156 LBM/FT3 
. 00003 
1 . 00 
3700 
.895 



1 . 03280 
. 92435 
.01252 LBM/SEC 
10.497 CFM 
3682 



127 



**ROTOMETER CALIBRATION** 
ASME HERSCHEL-TYPE VENTURI 
4 . 26 / 2 . 13 IN DIA 

35 PERCENT READING 



** EX PER I MENTAL DATA** 



ATMOSPHERIC PRESSURE = 29.76910 IN Hg 
VENTURI INLET TEMPERATURE = 72.0 DEG F 
VENTURI INLET STATIC GAGE PRESSURE = 17.50 IN H20 
VENTURI DIFFERENTIAL PRESSURE = .020 IN RED OIL 
SPECIFIC GRAVITY OF RED OIL = .834 



* CALCULATED VALUES** 



ATMOSPHERIC PRESSURE = 

VENTURI INLET GAGE PRESSURE = 
VENTURI INLET ABSOLUTE PRESSURE 
VENTURI DIFFERENTIAL PRESSURE = 
VENTURI INLET AIR DENSITY = 
VENTURI PRESSURE DROP RATIO (X) 
EXPANSION FACTOR (Y) = 

ASSUMED REYNOLDS NUMBER = 
DISCHARGE COEFFICIENT (C) = 
VELOCITY OF APPROACH FACTOR (E) 
FLOW COEFFICIENT <K) = 

AIR MASS FLOW RATE = 

AIR VOLUMETRIC FLOW RATE = 
VENTURI INLET REYNOLDS NUMBER. = 



14.6259 PSI 
.63219 PSI 
13.99367 PSI 
.00060 PSI 
.07105 LBM/FT3 
. 00004 
1 . 00 
4250 
. 899 
1 . 03280 
. 92797 
.01446 LBM/SEC 
12.211 CFM 
4253 
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**ROTOMETER CALIBRATION** 
ASME HERSCHEL-TYPE VENTURI 
4.26/2. 13 IN DIA 

40 PERCENT READING 



** EXPERIMENTAL DATA** 



ATMOSPHERIC PRESSURE = 29.76910 IN Hg 
VENTURI INLET TEMPERATURE = 72.0 DEG F 
VENTURI INLET STATIC GAGE PRESSURE = 20.30 IN H20 
VENTURI DIFFERENTIAL PRESSURE = .025 IN RED OIL 
SPECIFIC GRAVITY OF RED OIL = .834 



**CALCULATED VALUES** 



ATMOSPHERIC PRESSURE = 

VENTURI INLET GAGE PRESSURE = 
VENTURI INLET ABSOLUTE PRESSURE 
VENTURI DIFFERENTIAL PRESSURE = 
VENTURI INLET AIR DENSITY = 
VENTURI PRESSURE DROP RATIO (X) 
EXPANSION FACTOR <Y> = 

ASSUMED REYNOLDS NUMBER = 
DISCHARGE COEFFICIENT <C> = 
VELOCITY OF APPROACH FACTOR (E) 
FLOW COEFFICIENT <K> = 

AIR MASS FLOW RATE = 

AIR VOLUMETRIC FLOW RATE = 
VENTURI INLET REYNOLDS NUMBER = 



14.6259 PS I 
.73334 PS I 
13.89252 PSI 
.00075 PSI 
.07054 LBM/FT3 
. 00005 
1 . 00 
4750 
. 901 



1 . 03280 
. 93086 
.01616 LBM/SEC 
13.744 CFM 
4753 
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**ROTDMETER CALIBRATION** 
ASME HERSCHEL-TYPE VENTURI 
4.26/2. 13 IN DIA 

45 PERCENT READING 



**EXPER I MENTAL DATA** 



ATMOSPHERIC PRESSURE = 29.76910 IN Hg 
VENTURI INLET TEMPERATURE = 72.0 DEG F 
VENTURI INLET STATIC GAGE PRESSURE = 24.30 IN H20 
VENTURI DIFFERENTIAL PRESSURE = .031 IN RED OIL 
SPECIFIC GRAVITY OF RED OIL = .834 



* *CALCULATED VALUES** 



ATMOSPHERIC PRESSURE = 

VENTURI INLET GAGE PRESSURE = 
VENTURI INLET ABSOLUTE PRESSURE 
VENTURI DIFFERENTIAL PRESSURE = 
VENTURI INLET AIR DENSITY = 
VENTURI PRESSURE DROP RATIO (X) 
EXPANSION FACTOR (Y) = 

ASSUMED REYNOLDS NUMBER = 
DISCHARGE COEFFICIENT (C) = 
VELOCITY OF APPROACH FACTOR <E) 
FLOW COEFFICIENT <K) = 

AIR MASS FLOW RATE = 

AIR VOLUMETRIC FLOW RATE = 
VENTURI INLET REYNOLDS ‘NUMBER = 



14.6259 PSI 
.87784 PSI 
13.74802 PSI 
.00093 PSI 
.06981 LBM/FT3 
. 00007 
1 . 00 
5300 
. 904 
1.03280 
. 93365 
.01795 LBM/SEC 
15.431 CFM 
5281 
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**ROTOMETER CALIBRATION** 
ASME HERSCHEL-TYPE VENTURI 
4.26/2. 13 IN DIA 

50 PERCENT READING 



♦♦EXPERIMENTAL DATA** 



ATMOSPHERIC PRESSURE = 29.76910 IN Hg 
VENTURI INLET TEMPERATURE « 72.0 DEG F 
VENTURI INLET STATIC GAGE PRESSURE = 27.30 IN H20 
VENTURI DIFFERENTIAL PRESSURE = .037 IN RED OIL 
SPECIFIC GRAVITY OF RED OIL = .834 



♦♦CALCULATED VALUES** 



ATMOSPHERIC PRESSURE = 

VENTURI INLET GAGE PRESSURE = 
VENTURI INLET ABSOLUTE PRESSURE 
VENTURI DIFFERENTIAL PRESSURE = 
VENTURI INLET AIR DENSITY = 
VENTURI PRESSURE DROP RATIO (X) 
EXPANSION FACTOR (Y) = 

ASSUMED REYNOLDS NUMBER - 
DISCHARGE COEFFICIENT <C) = 
VELOCITY OF APPROACH FACTOR <E) 
FLOW COEFFICIENT <K) = 

AIR MASS FLOW RATE = 

AIR VOLUMETRIC FLOW RATE = 
VENTURI INLET REYNOLDS NUMBER - 



14.6259 PSI 
.98621 PSI 
13.63965 PSI 
.00111 PSI 
.06926 LBM/FT3 
. 00008 
1 . 00 
5750 
.906 
1 . 03280 
. 93592 
.01958 LBM/SEC 
16.967 CFM 
5760 
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**ROTOMETER CALIBRATION** 
ASME HERSCHEL-TYPE VENTURI 
4.26/2. 13 IN DIA 

55 PERCENT READING 



** EXPERIMENTAL DATA** 



ATMOSPHERIC PRESSURE = 29.76910 IN Hg 
VENTURI INLET TEMPERATURE = 72.0 DEG F 
VENTURI INLET STATIC GAGE PRESSURE = 31.85 IN H20 
VENTURI DIFFERENTIAL PRESSURE = .044 IN RED OIL 
SPECIFIC GRAVITY OF RED OIL = .834 



* *CALCULATED VALUES** 



ATMOSPHERIC PRESSURE = 

VENTURI INLET GAGE PRESSURE = 
VENTURI INLET ABSOLUTE PRESSURE 
VENTURI DIFFERENTIAL PRESSURE = 
VENTURI INLET AIR DENSITY = 
VENTURI PRESSURE DROP RATIO (X) 
EXPANSION FACTOR (Y) = 

ASSUMED REYNOLDS NUMBER = 
DISCHARGE COEFFICIENT <C> = 
VELOCITY OF APPROACH FACTDR (E) 
FLOW COEFFICIENT <K> = 

AIR MASS FLOW RATE = 

AIR VOLUMETRIC FLOW RATE = 
VENTURI INLET REYNOLDS NUMBER = 



14.6259 PSI 
1. 15058 PSI 
13.47528 PSI 
. 00 1 33 PS I 
.06842 LBM/FT3 
. 000 1 0 
1 . 00 
6250 
.909 



1 . 03280 
. 93830 
.02128 LBM/SEC 
18.662 CFM 
6259 
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**ROTOMETER CAL I BRAT I ON* * 
ASME HERSCHEL-TYPE VENTURI 
4.26/2. 13 IN DIA 

60 PERCENT READING 



^EXPERIMENTAL DATA** 



ATMOSPHERIC PRESSURE = 29.76910 IN Hg 
VENTURI INLET TEMPERATURE = 72.0 DEG F 
VENTURI INLET STATIC GAGE PRESSURE = 37.00 IN H20 
VENTURI DIFFERENTIAL PRESSURE = .051 IN RED OIL 
SPECIFIC GRAVITY OF RED OIL = .834 



* CALCULATED VALUES** 



ATMOSPHERIC PRESSURE = 

VENTURI INLET GAGE PRESSURE = 
VENTURI INLET ABSOLUTE PRESSURE 
VENTURI DIFFERENTIAL PRESSURE = 
VENTURI INLET AIR DENSITY = 
VENTURI PRESSURE DROP RATIO <X) 
EXPANSION FACTOR (Y) = 

ASSUMED REYNOLDS NUMBER = 
DISCHARGE COEFFICIENT (C) = 
VELOCITY OF APPROACH FACTOR (E) 
FLOW COEFFICIENT (K) = 

AIR MASS FLOW RATE = 

AIR VOLUMETRIC FLOW RATE = 
VENTURI INLET REYNOLDS NUMBER = 



14.6259 PSI 
1.33663 PSI 
13.28923 PSI 
.00154 PSI 
.06748 LBM/FT3 
. 00012 
1 . 00 
6700 
.911 
1 . 03280 
. 94036 
.02280 LBM/SEC 
20.277 CFM 
6707 
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APPENDIX C. EXPERIMENTAL HEAT TRANSFER DATA 



<D 



l n 
om 



cd 

cd 



CD 

<-M 

<f) 

Cd 

5-. 



> 

CD 

i — ! 

-C 

cd 



cm m 


U- 


<-M CD 


o 


0) *M 


i 


o a 


CM 


U CD 




IE 


\ u- 


• 


3 ' 


4M r: 


•M *-« 


in —j 


oq m 



VM 

4-t d) 
O 

2 u o 

m m • 



u- 

o 

I 

CM 

^ u- 



+-> m 4 -j m 



eT 

. CQ 

lE— • • ccc 

CO 

> 

< 



£• u 
E 0) 

CD 

— H s ° 

IE— * cO v — ' 

< 3 



Cl CD 

O > nJ 



*1 

3 



5 E 

o u 

r-H CO 
Cl- W 



J- 

0) 

-M CO 
CD C 

E 

O 03 

+J Cl / — N 
O Qj ©\° 
CZ 06 < ' 



v£> O vO I/O O 
r^rNsjocN^rN^mco 

CM fs 

• * • N H O CN >J o CO 

<r<r<r— 'CMco<r<ruouo 



^ O CO (O o\ 
«— 'OMOvON^vOCNON 
«— • on o\ ^ on • • • • • 

o <r r* a> 

cococovoctn^cmcmcmcm- 











a\ 






<JS 


a> 


<r 


4r 


00 


o 


CO 




CO 


a> 


<T 


00 


UO 


o 


v£> 


p- 


uo 


• 


<T\ 


• 


• 


• 


• 


• 


• 


• 


• 




• 


o 




o 


o 


oo 


a> 


uo 


CO 


O 


00 


o 


o 


o 


o 


r- 


CO 


00 






CT\ 






1— • 





CNMOCNONONCNnCN 

cnoocsicNvoocncNr^o 

^4<rr-uocMcoaooouo<* 



cd 




. 






















CD 




CD ^ 


























E a 


























CD <M 


UO 


o 


CO 


CO 


<r 


<r 


U0 


uO 


UO 








H 3 / — \ 


vO 


CO 


o 


00 


00 


ON 


--4 


V«0 


<r 




rH 


DC 


<D IM 






















cd 


. 2 


• mo 


CM 


CO 


CT> 


MT 


ON 


a> 


i—4 


r-* 


<r 


i—4 


+j 


IH 


CO ^ ’ 


r-* 


r*. 


v «o 


r^. 


UO 


CO 


CO 


CM 


CM 


CM 


c 

<d 




< s 






»—^ 




i—4 




«—4 






^■4 


E 




/ — \ 






















•H 




CM 


























+-> 






















CD 




CO 


CM 


CM 


UO 


CO 


uO 


i—4 




00 


o 


r- 




tr> +j i 


00 




<r 


ON 


vO 


r- 


uO 


CO 


UO 


CM 


CD 




rt -M ^h 






















X 




^ o m 


CM 


«— i 


CO 


—4 


o 


CM 


CM 


CM 


CO 


<r 


CU 


CJ 




i— 4 


CM 


CO 


<r 


UO 


vO 


r^. 


00 


ON 


o 






c/l rH £ 


•—4 


--4 


«— » 


—4 






i-H 






CM 



n co n oo 
<ra\<rcocoooin*-<co<r 
vO cni cr> uo cm oo • • • 

O ' CM 

OrNfNCO(jNl^-Hi-( — r-t 



ooocM<r^oooocN<rvo 

-HCNcNCNcNcNcncncnen 



134 



Table VI. Experimental Heat Transfer Data (6.25 in. Masked Bed) 
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Rotometer Air 
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Table VIII. Experimental Heat Transfer Data (10. 
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APPENDIX D 



SAMPLE THERMOCOUPLE CALIBRATION DATA 



THERMOCOUPLE # 14 



EMP 


PLH T 


TC 


BHTH 


c, C T 


PES 


1ST 


OUTPUT 


TEMP 


* * * 


44 4:444 


444444 


4 444 


l 


• j y 


4555 


“7 “7 “7 


1 8 . 


59 


2 


^ ,• 


sspif: 


. 774 


19 . 


61 




u o 


1879 


1 . 024 


25 . 


x 2 


4 


2 y 


8 83 8 


1.221 


7. M 


cr -7 

i . . 


5 


2 9 


4 589 


1 549 


7* O 


7 9 


6 


v Ci 


9557 


1 . 759 


43 . 


35 


r 


30 . 


f i - ’ C* 


2 . 107 


51 . 


4 9 




31 . 


3107 


2 . 34 0 


56 . 


o 5 


u 


u! 


1 88 7 


2 . i <. — f 


65 . 


8 3 


10 


*1 


7447 


2 . 974 


71 . 


20 


1 1 


7 7 


6 8 6 4 


3.401 


80 . 


85 


12 


"7 "7 


9723 


3.531 


i~i ~ 


cr -7 


1 3 


7 "7 


3729 


3 . 257 


i r 


50 


1 4 


j 


4592 


2 . 34 4 


. 8 8 . 


T 7 


15 


3 1 * 


9132 


2 . 805 


O %L 


Ci 2 


1 6 


31 


0491 


2 .225 


54 . 


24 


17 


3 O . 


5237 


2 002 


49 . 


05 


18 


2 9 


7 0 7 0 


1 . 652 


40 


ij 


1 9 


2 9 


2117 


1 . 444 


T cr 

!• ■ J . 


96 


20 


2 y . 


4275 


1.121 


2 y . 


19 



u 

cp 7 5 
0 ) 

D 



a 

e 

cu 



5 9 



5 1 

4J 

03 

cq 4 v 
3 4 
2 6 
1 8 




r- 



o-j »*-j '*-j «>j r-:> f\> 



Thermocouple Output (mv) 
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f: * X P 0 L V N 0 MI AL F U H C T I 0 N * * X 
< ORDER =4) 

V = - . 3 6 3 2 1 7 5 5 4 1 + 26 . 87 5 5 1 3 37 1 3 
X X - 1 . 5 3 5257363 8 * X 2 + 

. 306 27483 .£X A 3 - A .3 04 83 8 788 8 
tX* 4 



X V ( R C T > V ( L'.HLC ) 

ft. XX XXX XXXXXXXX 

733 18.53 13.5366 





. 7 i* 4 


13 
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c* ”7 

i . i 


7 0 


. 57 81 


1 


5 43 


*7 •!; 


7 s 


“7 jZj 


. 4253 


1 


. 7 53 


43 


“7 cr 


43 


. 3807 




.187 


51 


43 


51 


. 4334 




3 40 


5 6 


CT 

ii . 1 


56 


354 3 


2 


. » ’ 2 i'' 


65 


6 7 


65 


. 6616 


2 


. 374 


71 


. 20 


71 


. 2135 


T 


401 


S O 


. 65 


3 0 


. 67 17 


"7 


. 531 


t 


“7 


ft "7 


.511 5 


T 


257 


r r 


. 5 0 


i i'* 


. 5034 


2 


. 844 


6 S 


"7 "7 


6 3 


. 2384 




. 605 


4! 2 


Q 2 


r . 2 


3 3 3 3 


, * 


. 225 


54 


24 


54 


2127 


2 


. 0 0 2 


42 


05 


42 


. 0553 


1 


.652 


4 0 


'I 1 'J 


4 0 


8624 


1 


. 444 


“7 c" 

._i _i 


36 


~7 cr 

.Ji . j 


. 3230 


1 


121 


28 


. 13 


23 


. 1 554 


M 


OF 3 


QUA RE 


ERRORS = 






0 


12263 


527703 
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APPENDIX E. 



SAMPLE HEAT TRANSFER DATA REDUCTION COMPUTER PRINTOUTS 



1. Non- fluidized 6.25 in. configuration. 



^EXPERIMENTAL DATA** 
RUN # 06B20 
06/08/82 - 0425 



RUN NUMBER = 

BED WIDTH = 

STATIC BED HEIGHT = 

AIR FLOW RATE = 

AMBIENT TEMPERATURE = 
RIGHT HEATER VOLTAGE = 
RIGHT HEATER CURRENT = 
LEFT HEATER VOLTAGE = 
LEFT HEATER CURRENT = 
BED EXPANSION = 



06B20 
6.25 INCHES 
12.5 INCHES 
7.29 CFM 
74.50 DEG F 

44.0 VOLTS 
.82 AMPS 

45.0 VOLTS 
1.00 AMPS 

0.0 INCHES 



COMMENTS : NON-FLU I D I ZED/NON-EXPANDED 
BED WITH SMALL LOCAL SPOUT AT 10-10 
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1 

2 

•3 

4 

5 

6 

7 

8 

9 

10 

1 1 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

'—i 

jL 

23 

24 

'ne* 

xLU 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

33 

39 

40 



** EX PER I MENTAL DATA** 
RUN # 06 B 20 
06 / 08/82 - 0425 

THERMOCOUPLE READINGS 



OUTPUT 


TEMP 


TO# 


OUTPUT 


TEMP 


(mV) 


(DEG F) 




(mV) 


(DEG F) 



174.65 41 2.110 124.98 



sl ’> ■ o o 1 


174.55 


42 


2.076 


123. 54 


3 . 336 


174.76 


43 


1 . 958 


118. 62 


.977 


79.23 


44 


2.087 


123. 37 


3 . 330 


174. 55 


45 


2 . 112 


124.92 


3 . 337 


174.81 


46 


2 . 002 


120. 34 


3.343 


175. 05 


47 


.911 


73 . 51 


O . O O kJ 


174. 71 


48 


.911 


73 . 49 


3.312 


173.78 


49 


.910 


73.46 


3 . 322 


174. 17 


50 


1 . 430 


96 . 22 


3.319 


174. 07 


51 


1 . 408 


95 . 30 


3.316 


173. 94 


52 


1 . 385 


94.33 


3.278 


172. 36 


53 


. 953 


75 . 34 


3 . 282 


172. 50 


54 


. 948 


75 . 1 1 


3.287 


172. 72 


55 


1 . 202 


86 . 16 


3 . 286 


172. 60 


56 


1 . 422 


95 . 85 


3 • 230 


170. 36 


57 


1.417 


95 . 62 


3 . 235 


170. 55 


60 


2 . 035 


121. 70 


3.238 


170. 72 


69 


2 . 483 


1 40. 40 


O • 410^. 


170.58 


70 


1 . 952 


118. 20 


3 . 129 


166.36 


71 


2.111 


124. 90 


3 . 132 


166.51 


72 


1 . 890 


115.58 


3.127 


166.26 


73 


2 . 072 


123. 26 


3 . 129 


1 66 . 35 


74 


. 936 


74.29 


3 . 130 


166.41 


75 


. 940 


74.47 


3 . 131 


166.47 


76 


2 . 119 


125. 24 


3 . 128 


166.34 


77 


2.331 


136. 13 


3 . 123 


1 66 .31 


78 


1 . 872 


114.81 


3.110 


165. 59 


79 


2.023 


121 .41 


3.113 


165.70 








3.111 


1 65 . 63 








3 . 108 


165. 43 








3.077 


164.22 








3 . 075 


164.14 






' 


3.075 


164.11 








3.072 


163.96 








3.016 


161.74 








3.021 


161 . 93 








3 . 0 1 9 


1 6 1 .83 








3 . 020 


162.13 
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♦♦EXPERIMENTAL DATA* * 
RUN# 06B20 
06/08/32 - 0425 



RIGHT-HAND HEATER 
TEMPERATURE PROFILE 







COL 1 


COL 2 


COL 3 


COL 4 


AVE 


ROW 


1 


174.65 


174.55 


174. 76 


****** 


174. 


66 


ROW 


2 


174. 55 


174.81 


175. 05 


174. 71 


174. 


78 


ROW 


o 


173.78 


174. 17 


174.07 


173. 94 


173. 


99 


ROW 


4 


172. 36 


172. 50 


172. 72 


172. 60 


172. 


54 


ROW 


5 


1 70 . 36 


1 70 . 55 


170.72 


170. 58 


170. 


55 



NOTE: TEMPERATURES SHOWN ARE IN DEGREES FAHRENHEIT 
AND ARE ARRANGED AS SEEN FROM WITHIN THE BED 



RIGHT-HAND HEATER 
TEMP vs DEPTH 




TEMP (DEG F) 



NOTE: PLOTTED TEMPERATURES ARE AVERAGES OF EACH 
HORIZONTAL ROW 
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♦♦EXPERIMENTAL DATA** 
RUN# 06B20 
06/08/32 - 0425 



LEFT-HAND HEATER 
TEMPERATURE PROFILE 







COL 1 


COL 2 


COL 3 


COL 4 


AVE 




ROW 


1 


166. 36 


166.51 


166.26 


166.35 


166. 


~T — T 

•_’» / 


ROW 




166.41 


166.47 


166.34 


166.31 


166. 


38 


ROW 




165. 59 


165. 70 


1 65 . 63 


165.43 


165. 


59 


ROW 


4 


164.22 


164.14 


164.11 


163. 96 


1 64 . 


10 


ROW 


5 


161.74 


161 . 93 


161.83 


162. 13 


161 . 


91 



NOTE: TEMPERATURES SHOWN ARE IN DEGREES FAHRENHEIT 
AND ARE ARRANGED AS SEEN FROM WITHIN THE BED 



LEFT-HAND heater 
TEriF vs DEPTH 

X 

X 

X 

X 

X 

S iO .Si in S 

U"J ^ s*0 

TEMP COEG F) 



NOTE: PLOTTED TEMPERATURES ARE AVERAGES OF EACH 
HORIZONTAL ROW 
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NOTE: 



**EXPER I MENTAL DATA** 
RUN# 06B20 
06/08/82 - 0425 



FLUIDIZED BED 
TEMPERATURE PROFILE 



95.85 



95. 62 



75. 34 



75.11 



96.22 
95. 30 
94.33 



TEMPERATURES SHOWN ARE IN DEGREES FAHRENHEIT 
AND ARE ARRANGED AS SEEN FROM ABOVE THE BED 
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* ^CALCULATED RESULTS** 
RUN # 06B20 
06/08/82 - 0425 



RIGHT-HAND HEATER 



AVERAGE HEATER TEMPERATURE = 
AVERAGE PLATE SURFACE TEMPERATURE 
AVERAGE BED TEMPERATURE = 
ELECTRICAL ENERGY INTO HEATER = 
LOSS FROM HTR BACKING TO ATM = 
h ( HTR BACKING-ATM) = .52 

LOSS FROM HTR PERIMETER TO BED = 
HEAT FLUX THRU PLATE = 
h (PLATE-BED) = 3.39 



173.30 DEG F 
173.28 DEG F 
89.68 DEG F 
123.11 BTU/HR 
11.22 BTU/HR 
BTU/HR-FT-'-2-DEG F 
13.50 BTU/HR 
98.39 BTU/HR 
BTU/HR-FT-"'2-DEG F 



LEFT-HAND HEATER 



AVERAGE HEATER TEMPERATURE = 
AVERAGE PLATE SURFACE TEMPERATURE = 
AVERAGE BED TEMPERATURE = 
ELECTRICAL ENERGY INTO HEATER = 
LOSS FROM HTR BACKING TO ATM = 
h (HTR BACKING-ATM) = .42 

LOSS FROM HTR PERIMETER TO BED = 
HEAT FLUX THRU PLATE = 
h (PLATE-BED) = 4.77 



164.87 DEG F 
164.84 DEG F 
89.68 DEG G 
153.54 BTU/HR 
9.20 BTU/HR 
BTU/HR-FT" S 2-DEG F 
19.77 BTU/HR 
124.58 BTU/HR 
BTU/HR-FT" S 2-DEG F 



ENERGY BALANCE 



AVERAGE AIR INLET TEMPERATURE 
AIR OUTLET TEMPERATURE = 
AVERAGE BED TEMPERATURE = 
q ( AIR OUT) = 
q (LOSS F/R WALL) = 
q (TOTAL OUT OF BED) = 
q (TOTAL INTO BED) = 
SUPERFICIAL VELOCITY = 
SUPERFICIAL MASS VELOCITY = 
PARTICLE REYNOLDS NUMBER = 



73.49 DEG F 
86. 16 DEG F 
39.68 DEG F 
96.97 BTU/HR 
.50 BTU/HR 
97.47 BTU/HR 
256.23 BTU/HR 
.47 FT/SEC 
121.92 LBM/HR-FT--2 
2 . 44 
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2 . 



Fluidized 6.24 in. configuration. 



** EX PERT MENTAL DATA** 

RUN # 06B26 
06/08/82 - 0510 

RUN NUMBER = 06B26 



BED WIDTH = 


6.25 INCHES 


STATIC BED HEIGHT = 


12.5 INCHES 


AIR FLOW RATE = 


9.23 CFM 


AMBIENT TEMPERATURE = 


74.50 DEG F 


RIGHT HEATER VOLTAGE = 


60.0 VOLTS 


RIGHT HEATER CURRENT = 


1.12 AMPS 


LEFT HEATER VOLTAGE = 


60.0 VOLTS 


LEFT HEATER CURRENT = 


1.28 AMPS 


BED EXPANSION = 


. 1 INCHES 



COMMENTS: FLUIDIZATION WITH 2-3 INCH 
DIA BUBBLES ERUPTING AT LOCATION 5-5 
AND VICINITY/ACTIVE BOT LEFT 



147 



♦♦EXPERIMENTAL DATA** 
RUN # 06B26 
06/08/92 - 0510 

THERMOCOUPLE READINGS 



TO# 


OUTPUT 


TEMP 


TO# 


OUTPUT 


TEMP 




(mV) 


(DEG F) 




(mV) 


(DEG F) 


1 


2. 729 


150. 44 


41 


2. 146 


126. 48 


X. 


2.728 


150.36 


42 


2. 1 16 


125.21 


o 


2.724 


150. 13 


43 


1 . 996 


120.22 


4 


. 961 


78. 19 


44 


2. 160 


126.91 


5 


2 . 855 


1 55 . 53 


45 


2. 154 


126.67 


6 


2.859 


155. 67 


46 


2. 041 


121.97 


7 


2.858 


155. 63 


47 


. 922 


74. 00 


8 


2.848 


155.21 


48 


. 922 


73.98 


9 


2. 991 


1 60 . 95 


49 


. 922 


73. 99 


10 


3 . 000 


161.31 


50 


1 . 564 


101 . 98 


1 1 


2. 996 


161 . 15 


5 1 


1 . 555 


101 . 63 


12 


2. 988 


1 60 . 85 


52 


1.551 


101 . 49 


13 


3. 107 


165. 53 


53 


1 . 556 


101.66 


14 


3. 106 


165.48 


54 


1 . 588 


103. 00 


15 


3. 106 


165.49 


55 


1 . 526 


1 00 . 1 2 


16 


3 . 1 06 


165. 43 


56 


1 . 551 


101 . 40 


17 


3. 141 


166. 80 


57 


1 . 550 


101.35 


18 


3.144 


166.92 


68 


1 . 927 


117.14 


19 


3. 140 


166.81 


69 


1 . 958 


118.46 


20 


3. 137 


166. 79 


70 


1 . 930 


117.27 


21 


2. 259 


1 30 . 99 


71 


2.368 


135.64 




2. 290 


132. 28 


72 


1 .916 


116.68 


'TT 


2.291 


132.23 


73 


2. 039 


121 . 37 


24 


2. 297 


132.51 


74 


1 . 045 


79 . 1 0 


'"'i C=- 
wl 


2. 340 


134. 30 


75 


1 . 260 


88. 51 


26 


2. 336 


1 34 . 1 6 


76 


1 . 839 


113.41 


27 


2 . 335 


134. 10 


77 


1 . 356 


114.14 


23 


2. 343 


1 34 . 63 


73 


1 . 673 


1 06 . 34 


29 


2.346 


134.56 


79 


1 . 925 


117.06 


30 


2 . 335 


1 34 . 06 








31 


2. 339 


134. 26 








32 


2.354 


134.79 








• j* O 


2.293 


132. 51 








34 


2. 285 


131 . 94 








35 


2. 291 


132.13 








36 


2. 310 


132.97 








;t"7 




129.73 








33 


2.213 


129. 16 








39 




129.44 








40 


ncn 
^ i ^ J7 


131.11 
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# ♦EXPERIMENTAL DATA* * 
RUN# 06B26 
06/08/32 - 0510 



RIGHT-HAND HEATER 
TEMPERATURE PROFILE 







COL 1 


COL 2 


COL 3 


COL 4 


AVE 


ROW 


1 


150.44 


1 50 . 36 


150.18 


****** 


1 50 . 


33 


ROW 


2 


155. 53 


155.67 


155. 63 


155. 21 


155. 


51 


ROW 


3 


1 60 . 95 


161 . 31 


161.15 


160. 85 


161 . 


07 


ROW 


4 


165.53 


165. 48 


165. 49 


165. 43 


165. 


48 


ROW 


5 


166. 80 


166.92 


166.81 


166. 79 


166. 


83 



NOTE: TEMPERATURES SHOWN ARE IN DEGREES FAHRENHEIT 
AND ARE ARRANGED AS SEEN FROM WITHIN THE BED 



RIGHT-KRND HER TER 
TEMP vs D E F T H 

% 

X 
% 

% 

X 

sms m s m 

m m vD ^ n- n. 

TEMP (DEG F ) 

NOTE: PLOTTED TEMPERATURES ARE AVERAGES OF EACH 
HORIZONTAL ROW 



z I. 0 
C 9 

8 

l 7 

— 6 

" 5 

4 
3 
2 
t 



co 

o 

a. 



a 
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^•EXPERIMENTAL DATA** 
RUN# 06B26 
06/08/82 - 0510 



LEFT-HAND HEATER 
TEMPERATURE PROFILE 







COL 1 


COL 2 


COL 3 


COL 4 


AVE 


ROW 


1 


1 30 . 99 


132.23 


132. 23 


132. 51 


132. 02 


ROW 




1 34 . 30 


134.16 


134.1 0 


1 34 . 63 


134.30 


ROW 


o 


1 34 . 56 


1 34 . 06 


1 34 . 26 


1 34 . 79 


134.42 


ROW 


4 


132. 51 


131 . 94 


132. 18 


132. 97 


132. 40 


ROW 


5 


129.73 


129. 16 


129.44 


131.11 


129.37 


JOTE: 


TEMPERATURES SHOWN 


ARE IN 


DEGREES FAHRENHEIT 



AND ARE ARRANGED AS SEEN FROM WITHIN THE BED 



LEF T-HAND HEATER 
TEMP vs DEPTH 

* 

* 

* . 

* 

* 

iS -T <3. LO <S D 

oj cm ro ro ^ 

TFtlF (DEG F ) 

NOTE: PLOTTED TEMPERATURES ARE AVERAGES OF EACH 
HORIZONTAL ROW 
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NOTE: 



* *EXPER I MENTAL DATA** 
RUN# 06B26 
06/08/82 - 0510 



FLUIDIZED BED 
TEMPERATURE PROFILE 



101 . 40 



101 . 66 



101.98 
101.63 
101 . 49 



101.35 103.00 



TEMPERATURES SHOWN ARE IN DEGREES FAHRENHEIT 
AND ARE ARRANGED AS SEEN FROM ABOVE THE BED 
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* ^CALCULATED RESULTS** 
RUN # 06B26 
06/08/82 - 0510 



RIGHT-HAND HEATER 



AVERAGE HEATER TEMPERATURE = 
AVERAGE PLATE SURFACE TEMPERATURE 
AVERAGE BED TEMPERATURE = 
ELECTRICAL ENERGY INTO HEATER = 
LOSS FROM HTR BACKING TO ATM = 
h (HTR BACKING-ATM) = .36 

LOSS FROM HTR PERIMETER TO BED = 
HEAT FLUX THRU PLATE = 
h (PLATE-BED) = 9.97 



159.34 DEG F 

159.79 DEG F 

101.79 DEG F 
229.29 BTU/HR 

7.92 BTU/HR 
BTU/HR-FT"'-2-DEG F 
20.53 BTU/HR 
200.79 BTU/HR 
BTU/HR-FT--2-DEG F 



LEFT-HAND HEATER 



AVERAGE HEATER TEMPERATURE = 

AVERAGE PLATE SURFACE TEMPERATURE = 

AVERAGE BED TEMPERATURE = 

ELECTRICAL ENERGY INTO HEATER = 

LOSS FROM HTR BACKING TO ATM = 
h (HTR BACKING-ATM) = .07 BTU/HR-FT'"2-DEG F 

LOSS FROM HTR PERIMETER TO BED = 33.49 BTU/HR 

HEAT FLUX THRU PLATE = 227.02 BTU/HR 

h (FLATE-BED) = 21.26 BTU/HR-FT""2-DEG F 



132.60 DEG F 
132.54 DEG F 
101.79 DEG G 
62.05 BTU/HR 
1.54 BTU/HR 



ENERGY BALANCE 



AVERAGE AIR INLET TEMPERATURE 
AIR OUTLET TEMPERATURE = 
AVERAGE BED TEMPERATURE = 
q (AIR OUT) = 
q ( LOSS F/R WALL) = 
q (TOTAL OUT OF BED) = 
q (TOTAL INTO BED) = 
SUPERFICIAL VELOCITY = 
SUPERFICIAL MASS VELOCITY = 
PARTICLE REYNOLDS NUMBER = 



73.99 DEG F 
100.12 DEG F 
101.79 DEG F 
247.04 BTU/HR 
26.35 BTU/HR 
273.39 BTU/HR 
481.33 BTU/HR 
.59 FT/SEC 
150.65 LBM/HR-FT "2 
3 . 0 1 
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3. 



Non- fluidized 12.0 in. configuration. 



^^EXPERIMENTAL DATA** 
RUN # 12B30E 
06/17/82 - 0055 



RUN NUMBER = 

BED WIDTH = 

STATIC BED HEIGHT = 

AIR FLOW RATE = 

AMBIENT TEMPERATURE = 
RIGHT HEATER VOLTAGE = 
RIGHT HEATER CURRENT = 
LEFT HEATER VOLTAGE = 
LEFT HEATER CURRENT = 
BED EXPANSION = 



12B30E 
12.00 INCHES 
12.5 INCHES 
10.53 CFM 
78.60 DEG F 

46.0 VOLTS 
.85 AMPS 

56.0 VOLTS 
1.02 AMPS 

0.0 INCHES 



COMMENTS: NON-FLUIDI ZED/NON-EXPANDED 
BED/NO BUBBLES BOTTOM LEFT 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 1 

12 

1 3 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 



** EX PER I MENTAL DATA** 
RUN # 12B30E 
06/17/82 - 0055 

THERMOCOUPLE READINGS 



OUTPUT 


TEMP 


TC# 


OUTPUT 


TEMP 


(mV) 


(DEG F) 




(mV) 


(DEG F) 


O . t_J U 6 


183. 47 


41 


2. 123 


125.52 


3.550 


183. 20 


42 


2. 145 


126. 42 


w* . J JO 


183. 34 


43 


2.041 


122. 10 


1 . 140 


89. 84 


44 


1.785 


111. 20 


3.554 


183. 39 


45 


1 . 862 


114.45 


o . j j9 


183. 55 


46 


1.795 


111.65 


o . j6 j 


1 83 . 83 


47 


.967 


75.99 


3.556 


183. 41 


48 


. 966 


75. 92 


3. 541 


182.84 


49 


. 966 


75. 94 


3. 551 


183.22 


50 


1 . 438 


96.56 


3.545 


183.03 


5 1 


1.429 


96.21 


3.541 


182. 80 


52 


1 . 403 


95. 1 1 


3.512 


181 . 62 


53 


1 . 273 


89.41 


3. 515 


181.69 


54 


1.277 


39.56 


3.516 


181.78 


55 


1.257 


88. 55 


3.514 


181 . 61 


56 


1.261 


88 . 86 


3. 468 


179.84 


57 


1 . 260 


88 . 80 


3.471 


179.90 


58 


1.265 


89.01 


3.468 


179.81 


59 


1.274 


89.40 


3.458 


179.52 


60 


1.255 


88. 87 


2.886 


156.60 


61 


1.257 


88. 95 


2.886 


1 i_j6 . 6 o 


62 


1.262 


89.15 


2. 883 


156.46 


63 


1 . 239 


88. 13 


2. 886 


156.58 


68 


1.954 


118. 29 


2. 893 


156.89 


69 


2. 477 


140. 15 


2.894 


- 156.94 


70 


2. 173 


127. 50 


2. 891 


156.82 


71 


2.527 


142.22 


2.890 


156.75 


72 


1.974 


119.13 


2.886 


156.59 


73 


2 . 200 


128.63 


2. 888 


156.65 


74 


1 . 095 


81.29 


2. 887 


156.62 


75 


1.165 


84 . 36 


2. 886 


156.51 


76 


1.977 


119.26 


2.869 


155. 85 


77 


2. 092 


124. 10 


2. 867 


155.76 


78 


1.773 


110.61 


2.869 


155.81 


79 


1 . 740 


109. 20 


2. 866 


155.67 








2.828 


154. 16 








2.836 


154.47 








jL . O— 'Jj 


1 54 . 30 








2. 831 


154.51 . 









154 



**EXPER I MENTAL DAT A* * 
RUN# 12B30E 
06/17/82 - 0055 



RIGHT-HAND HEATER 
TEMPERATURE PROFILE 







COL 1 


COL 2 


COL 3 


COL 4 


AVE 




ROW 


1 


183. 47 


183. 20 


1 83 . 34 


****** 


183. 


34 


ROW 




183. 39 


1 83 . 55 


183. 83 


183.41 


183. 


55 


ROW 


3 


182. 84 


183.22 


183.03 


182.80 


182. 


97 


ROW 


4 


131.62 


181.69 


181 . 78 


181.61 


181 . 


68 


ROW 


5 


179.84 


179.90 


179. 81 


179. 52 


179. 


77 



NOTE: TEMPERATURES SHOWN ARE IN DEGREES FAHRENHEIT 
AND ARE ARRANGED AS SEEN FROM WITHIN THE BED 



R I GHT-HRND HEfiTER 
TEMP vs DEPTH 




NOTE: PLOTTED TEMPERATURES ARE AVERAGES OF EACH 
HORIZONTAL ROW 
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**EXF'ER I MENTAL DATA* * 
RUN# 12B30E 
06/17/82 - 0055 



LEFT-HAND HEATER 
TEMPERATURE PROFILE 







COL 1 


COL 2 


COL 3 


COL 4 


AVE 




ROW 


1 


156. 6 V 


1 56 .63 


156.46 


156.58 


156. 


57 


ROW 




156. 39 


156.94 


156.32 


136. 75 


156. 


35 


ROW 


3 


156. 59 


156. 65 


156.62 


156. 51 


156. 


59 


ROW 


4 


155.35 


155. 76 


155. 81 


155. 67 


155. 


77 


ROW 


5 


154.16 


154. 47 


1 54 . 30 


154.51 


154. 


36 



NOTE: TEMPERATURES SHOWN ARE IN DEGREES FAHRENHEIT 
AND ARE ARRANGED AS SEEN FROM WITHIN THE BED 



LEFT-HAND HEATER 
TEMP vs DEPTH 




TEMP (DEG F > 



NOTE: PLOTTED TEMPERATURES ARE AVERAGES OF EACH 
HORIZONTAL ROW 
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^^EXPERIMENTAL DATA* * 
RUN# 12B30E 
06/17/82 - 0055 



FLUIDIZED BED 
TEMPERATURE PROFILE 



38. 95 


89.01 


88.36 


39.41 


96.56 


89. 15 


89 . 40 






96.21 


88. 13 


88.37 


88.80 


89.56 


95. 1 1 



NOTE: TEMPERATURES SHOWN ARE IN DEGREES FAHRENHEIT 
AND ARE ARRANGED AS SEEN FROM ABOVE THE BED 
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♦♦CALCULATED RESULTS* * 
RUN # 12B30E 
06/17/82 - 0055 



RIGHT-HAND HEATER 



AVERAGE HEATER TEMPERATURE = 
AVERAGE PLATE SURFACE TEMPERATURE 
AVERAGE BED TEMPERATURE = 
ELECTRICAL ENERGY INTO HEATER = 
LOSS FROM HTR BACKING TO ATM = 
h (HTR BACKING-ATM) = .61 

LOSS FROM HTR PERIMETER TO BED = 
HEAT FLUX THRU PLATE = 
h <PLATE-3ED) = 3.33 



182.26 DEG F 
182.23 DEG F 
90.62 DEG F 
133.41 BTU/HR 
12.85 BTU/HR 
BTU/HR-FT"‘2-DEG F 
14.77 BTU/HR 
105.80 BTU/HR 
BTU/HR-FT^C-DEG F 



LEFT-HAND HEATER 



AVERAGE HEATER TEMPERATURE = 
AVERAGE PLATE SURFACE TEMPERATURE 
AVERAGE BED TEMPERATURE » 
ELECTRICAL ENERGY INTO HEATER = 
LOSS FROM HTR BACKING TO ATM = 
h (HTR BACKING-ATM) = .63 

LOSS FROM HTR PERIMETER TO BED = 
HEAT FLUX THRU PLATE = 
h <PLATE-BED) = 7.29 



156.03 DEG F 
155.98 DEG F 
90.62 DEG G 
194.90 BTU/HR 
9.90 BTU/HR 
BTU/HR— FT' N 2— DEG F 
19.57 BTU/HR 
165.43 BTU/HR 
BTU/HR-FT''2-DEG F 



ENERGY BALANCE 



AVERAGE AIR INLET TEMPERATURE 
AIR OUTLET TEMPERATURE = 
AVERAGE BED TEMPERATURE = 
q (AIR OUT) = 
q ( LOSS F./R WALL) = 
q (TOTAL OUT OF BED) = 
q (TOTAL INTO BED) = 
SUPERFICIAL VELOCITY = 
SUPERFICIAL MASS VELOCITY = 
PARTICLE REYNOLDS NUMBER = 



75.95 DEG F 
88.55 DEG F 
90.62 DEG F 
138.71 BTU/HR 
16.49 BTU/HR 
155.20 BTU/HR 
305.56 BTU/HR 
.35 FT/SEC 
91.38 LBM/HR— FT"-2 
1 . 83 
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Fluidized 12.0 in. configuration. 



4 . 



** EX PER I MENTAL DATA** 

RUN # 12B48E 
06/17/82 - 0155 

RUN NUMBER = 12B48E 



BED WIDTH - 


12-00 INCHES 


STATIC BED HEIGHT = 


12.5 INCHES 


AIR FLOW RATE = 


16.37 CFM 


AMBIENT TEMPERATURE = 


81.50 DEG F 


RIGHT HEATER VOLTAGE = 


46.0 VOLTS 


RIGHT HEATER CURRENT = 


.85 AMPS 


LEFT HEATER VOLTAGE = 


46.0 VOLTS 


LEFT HEATER CURRENT = 


1.02 AMPS 


BED EXPANSION = 


0. 1 INCHES 


COMMENTS: FLUIDIZATION 


WITH 2-3 INCH 



BUBBLES ERUPTING AT LOCATION 4-4 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 1 

12 

1 3 

14 

15 

16 

17 

18 

19 

20 

21 

nn 

r—y T 

.c •-> 

24 

25 

26 

27 

28 

29 

30 

3 1 

32 

.jvO 

34 

35 

36 

37 

33 

39 

40 





^^EXPERIMENTAL DATA** 








RUN # 12B48E 








06/17/82 


- 0155 








THERMOCOUPLE 


READINGS 






OUTPUT 


TEMP 


TC# 


OUTPUT 


TEMP 


(mV) 


(DEG F) 




(mV) 


(DEG F) 


3.471 


180.14 


41 


2. 279 


132. 01 


3. 474 


180. 21 


42 


2. 244 


1 30 . 53 


3.478 


1 80 . 39 


43 


2. 097 


124. 44 


1 . 083 


86.14 


44 


2. 121 


125. 29 


3. 472 


180. 16 


45 


2.136 


125. 92 


3.479 


180.41 


46 


2. 024 


121 . 26 


3. 484 


1 80 . 63 


47 


. 959 


75. 64 


3. 476 


180.27 


48 


. 959 


75.62 


3. 445 


179.05 


49 


. 960 


75. 67 


3. 455 


179.43 


50 


1 . 073 


80. 65 


3.450 


179.27 


51 


1.122 


82. 84 


3. 445 


1 79 . 04 


52 


1 . 044 


79.41 


3. 394 


176. 96 


53 


. 989 


76. 94 


3. 399 


177.13 


54 


1 . 089 


81.33 


3. 398 


177. 12 


55 


1.128 


82. 94 


3 . 393 


176. 84 


56 


1 . 097 


81 . 67 


o . 


174.67 


57 


1 . 1 00 


81. 79 


3 . 343 


174.84 


58 


1 . 1 00 


81.79 


3 . 339 


174.72 


59 


1 . 1 05 


82. 01 


3. 327 


174.35 


60 


1 . 096 


81 . 90 


2. 831 


154.38 


61 


1.167 


85.01 


2. 833 


154.49 


62 


1 . 096 


81 . 37 


2. 833 


154.44 


63 


1 . 095 


81.82 


2. 839 


154.68 


68 


2. 133 


125. 82 


2.866 


155. 80 


69 


2. 600 


145.23 


2. 866 


155.81 


70 


2. 088 


123. 93 


2. 864 


155. 73 


71 


2. 571 


1 44 . 03 


2.864 


1 55 . 70 


72 


1.946 


117.95 


2.869 


155. 90 


73 


2.063 


122.88 


2. 867 


1 55 . 80 


74 


1 . 031 


78.48 


2.863 


155.65 


75 


1.021 


73 . 04 


2 . 863 


155.58 


76 


2. 024 


121 . 24 


2. 833 


154.39 


77 


2.246 


1 30. 55 


2.827 


154.14 


78 


1 . 806 


112.01 


2.821 


153.87 


79 


2.033 


121 . 62 


2.817 


153.69 








2.7 66 


151.65 








2. 764 


151.55 








2 . 756 


151.22 








2. 755 


151.43 
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O'N'0'rMMn<t'rt<tK)«tMHhrt o »-• n r-i cn n r-o n m ro to ** <*■ in 
r-i a- n -o -o <i o oo ^ O' n a- o n n o O' o ro co ai n o~ o o- ro o r i in 



♦♦EXPERIMENTAL DATA* * 
RUN# 1 2B43E 
06/17/82 - 0155 



RIGHT-HAND HEATER 
TEMPERATURE PROFILE 







COL 1 


COL 2 


COL 3 


COL 4 


AVE 


ROW 


1 


1 80 . 1 4 


180. 21 


180. 39 


****** 


180. 


24 


ROW 




180. 16 


180.41 


180. 63 


180. 27 


180. 


37 


ROW 


j 


179.05 


179.43 


179.27 


179.04 


179. 


20 


ROW 


4 


176.96 


177. 13 


177. 12 


176. 84 


177. 


01 


ROW 


5 


174.67 


174.84 


174.72 


174.35 


174. 


65 



NOTE: TEMPERATURES SHOWN ARE IN DEGREES FAHRENHEIT 
AND ARE ARRANGED AS SEEN FROM WITHIN THE BED 



RIGHT-HRND HEATER 
TEMP vs DEPTH 




TEMP (DEG F) 



NOTE: PLOTTED TEMPERATURES ARE AVERAGES OF EACH 
HORIZONTAL ROW 
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* >fcEXPER 1 MENTAL DATAy^ 
RUN# 1 2B43E 
06/17/32 O 1 5 5 



LEFT-HAND HEATER 
TEMPERATURE PROFILE 







COL 1 


COL 2 


COL 3 


COL 4 


MV E 




ROW 


1 


154.33 


154.49 


154.44 


154.63 


154. 


50 


ROW 




155. 30 


155. 31 


1 —i5 . 73 


1 55 .70 


1 =• c- 
1 JwJ. 


76 


ROW 


Li 


155. 90 


1 55 . SO 


155. 65 


155. 53 


155. 


73 


ROW 


4 


154.39 


154.14 


153.37 


153. 69 


154. 


02 


ROW 


5 


151.65 


151.55 


151 . 22 


151 . 43 


151 . 


46 



NOTE: TEMPERATURES SHOWN ARE IN DEGREES FAHRENHEIT 
AND ARE ARRANGED AS SEEN FROM WITHIN THE BED 



LEFT-HAND HEATER 
TEMP vs DEPTH 



2 1 



z 

o 



01 

o 

Ci- 



ct 

UJ 




TEMP (DEG F) 



NOTE: PLOTTED TEMPERATURES ARE AVERAGES OF EACH 
HORIZONTAL ROW 
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NOTE: 



* ♦EXPERIMENTAL DATA* * 
RUN# 12B48E 
06/17/32 - 0155 



FLUIDIZED BED 
TEMPERATURE PROFILE 



35.01 


31.79 


31.67 


76.94 


80.65 


31.37 


32. 01 






32.84 


31.32 


31.90 


31.79 


81.33 


79.41 



TEMPERATURES SHOWN ARE IN DEGREES FAHRENHEIT 
AND ARE ARRANGED AS SEEN FROM ABOVE THE BED 
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**CALCULATED RESULTS* * 
RUN # 12B48E 
06/17/82 - 0155- 



RIGHT-HAND HEATER- 



AVERAGE HEATER TEMPERATURE = 
AVERAGE PLATE SURFACE TEMPERATURE 
AVERAGE BED TEMPERATURE = 
ELECTRICAL ENERGY INTO HEATER = 
LOSS FROM HTR BACKING TO ATM = 
h (HTR BACKING-ATM) = .51 

LOSS FROM HTR PERIMETER TO BED = 
HEAT FLUX THRU PLATE = 
h < PLATE-BED) = 3.17 



178.29 DEG F 
178.27 DEG F 
81.46 DEG F 
133.41 BTU/HR 
10.82 BTU/HR 
BTU/HR-FT"'2-DEG F 
16. 12 BTU/HR 
106.48 BTU/HR 
BTU/HR -FT-''2-DEG F 



LEFT-HAND HEATER 



AVERAGE HEATER TEMPERATURE = 
AVERAGE PLATE SURFACE TEMPERATURE 
AVERAGE BED TEMPERATURE = 
ELECTRICAL ENERGY INTO HEATER = 
LOSS FROM HTR BACKING TO ATM = 
h (HTR BACKING-ATM) = .35 

LOSS FROM HTR PERIMETER TO BED = 
HEAT FLUX THRU PLATE = 
h (PLATE-BED) = 5.10 



154.29 DEG F 
154.26 DEG F 
81.46 DEG G 
160.10 BTU/HR 
6.44 BTU/HR 
BTU/HR-FT"'2-DEG F 
24.67 BTU/HR 
128.98 BTU/HR 
BTU/HR-FT'-2-DEG F 



ENERGY BALANCE 



AVERAGE AIR INLET TEMPERATURE 
AIR OUTLET TEMPERATURE = 
AVERAGE BED TEMPERATURE = 
q (AIR OUT) = 
q (LOSS F/R WALL) = 
q (TOTAL OUT OF BED) = 
q (TOTAL INTO BED) ‘ = 
SUPERFICIAL VELOCITY = 
SUPERFICIAL MASS VELOCITY = 
PARTICLE REYNOLDS NUMBER = 



75.64 DEG F 
82.94 DEG F 
31.46 DEG F 
126.26 BTU/HR 
0.00 BTU/HR 
126.26 BTU/HR 
276.25 BTU/HR 
.55 FT/SEC 
143.54 LBM/HR-FT" 2 
2. 37 
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5. A typical data run. (Bed thermocouple probe #57 reads 
significantly higher than other bed thermocouples. 
Reading exceeds highest heater thermocouple reading . ) 

NOTE: ALL VALUES WERE RECALCULATED AFTER OMITTING 

THERMOCOUPLE NR. 57 READING. REVISED RESULTS 
ARE PRESENTED IN TABLES AND CURVES PREVIOUSLY 
SHOWN . 

**EXPER I MENTAL DATA** 

RUN # 06B36M 
06/08/82 - 1630 

RUN NUMBER = 06B36M 



BED WIDTH = 


6.25 INCHES 


STATIC BED HEIGHT = 


12.5 INCHES 


AIR FLOW RATE = 


12.48 CFM 


AMBIENT TEMPERATURE = 


74.50 DEG F 


RIGHT HEATER VOLTAGE = 


60.0 VOLTS 


RIGHT HEATER CURRENT = 


1. 12 AMPS 


LEFT HEATER VOLTAGE = 


60.0 VOLTS 


LEFT HEATER CURRENT = 


1.28 AMPS 


BED EXPANSION = 


0.25 INCHES 



COMMENTS: VIOLENT FLUIDIZATION AT CARRYOVER 
LIMIT 
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1 

n 

4m 

o 

4 

5 

6 

7 

8 

9 

10 

11 
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13 

14 
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18 
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21 
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23 

24 
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* ^EXPERIMENTAL DATA** 
RUN # 06B36M 
06/08/82 - 1630 



THERMOCOUPLE READINGS 



OUTPUT 


TEMP 


TC# 


OUTPUT 


TEMP 


(mV) 


(DEG F) 




(mV) 


(DEG F) 


1.930 


117. 46 


41 


1.773 


110.81 


1 . 936 


117.68 


42 


1.774 


110.85 


1. 951 


118.31 


43 


1.726 


108.83 


.980 


79 . 43 


44 


1 . 966 


118.82 


2 . 004 


120.53 


45 


2.108 


124.75 


2.004 


120.50 


46 


2.007 


120.55 


2. 005 


120.57 


47 


.936 


74.62 


2 . 000 


120.34 


48 


. 935 


74.56 


2.096 


124.34 


. 49 


.935 


74. 57 


2.089 


124.06 


50 


1 . 631 


104. 85 


2.081 


123. 72 


51 


1.592 


103. 22 


2.088 


123. 96 


52 


1.621 


104.49 


2. 180 


127.75 


53 


1 . 608 


103. 89 


2. 185 


127. 92 


54 


1.609 


103.90 


2. 164 


127.07 


55 


1.579 


102.38 


2.166 


127. 11 


56 


1.629 


104.75 


2.209 


128.38 


57 


2.260 


131.38 


2.219 


129.27 


68 


1.582 


102. 44 


2.214 


1 29 . 09 


69 


1.658 


105.70 


2.208 


128.96 


70 


1.942 


1 17.73 


1.936 


117. 53 


71 


2.025 


121 . 28 


1.948 


118. 03 


72 


1 . 565 


101.71 


1 . 947 


117. 96 


73 


1.626 


104.32 


1.944 


117.82 


74 


1.313 


90.81 


2 . 005 


1 20 . 38 


75 


1.497 


98. 78 


1 . 987 


119.65 


76 


1.768 


110.39 


1.997 


120. 04 


77 


1.705 


107. 71 


2.017 


120.87 


73 


1.682 


106. 72 


2 . 061 


122.74 


79 


1.740 


109.20 


2 . 053 


122.35 








2.051 


122.30 








2 . 056 


122. 44 


i 






2.071 


123. 09 








2 . 067 


122. 88 








2. 075 


123.21 








2.079 


1 23 . 38 








2. 080 


123. 42 








2 . 058 


122.50 








2 . 055 


122. 37 








2.094 


124.25 
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* ^EXPERIMENTAL DATA** 
RUN# 06B36M 
06/08/82 - 1630 



RIGHT-HAND HEATER 
TEMPERATURE PROFILE 







COL 1 


COL 2 


COL 3 


COL 4 


AVE 


ROW 


1 


117.46 


117.68 


118.31 


****** 


117.82 


ROW 




120. 53 


120.50 


120.57 


120. 34 


120.48 


ROW 


3 


124. 34 


124.06 


123. 72 


123. 96 


124.02 


ROW 


4 


127.75 


127. 92 


127.07 


127. 11 


127. 46 


ROW 


5 


128. 88 


129.27 


129.09 


128.96 


129.05 



NOTE: TEMPERATURES SHOWN ARE IN DEGREES FAHRENHEIT 
AND ARE ARRANGED AS SEEN FROM WITHIN THE BED 



R I GHT-HAND HEATER 
TEMP vs DEPTH 




NOTE: PLOTTED TEMPERATURES ARE AVERAGES OF EACH 
HORIZONTAL ROW 
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^EXPERIMENTAL DATA* * 
RUN# 06B36M 
06/08/82 - 1630 



LEFT-HAND HEATER 
TEMPERATURE PROFILE 







COL 1 


COL 2 


COL 3 


COL 4 


AVE 


ROW 


1 


117.53 


1 18.03 


117.96 


117. 32 


117.34 


ROW 


jL 


1 20 . 38 


119.65 


120. 04 


120.87 


120.24 


ROW 


o 


122. 74 


122. 35 


122. 30 


122.44 


122.46 


ROW 


4 


123. 09 


122. 88 


123. 21 


123.38 


123.14 


ROW 


5 


123. 42 


122. 50 


122.37 


124.25 


123.1 3 


NOTE: 


TEMPERATURES SHOWN 


ARE IN 


DEGREES FAHRENHEIT 



AND ARE ARRANGED AS SEEN FROM WITHIN THE BED 



LEFT-HAND HERTER 
TEMP vs DEPTH 




NOTE: PLOTTED TEMFERATURES ARE AVERAGES OF EACH 
HORIZONTAL ROW 
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NOTE; 



**EXPER I MENTAL DATA** 
RUN# 06B36M 
06/08/82 - 1630 



FLUIDIZED BED 
TEMPERATURE PROFILE 



! 104.75 103.89 104.85 ! 

I I 

I I 

» I 

> I 

I 103.22 ! 

l I 

l l 

I I 

l I 

! 131.38 103.90 104.49 ! 

I I 

I I 

» I 

I I 

TEMPERATURES SHOWN ARE IN DEGREES FAHRENHEIT 
AND ARE ARRANGED AS SEEN FROM ABOVE THE BED 



169 



**CALCULATED RESULTS* * 
RUN # 06B36M 
06/08/82 - 1630 



RIGHT-HAND HEATER 



AVERAGE HEATER TEMPERATURE = 
AVERAGE PLATE SURFACE TEMPERATURE 
AVERAGE BED TEMPERATURE = 
ELECTRICAL ENERGY INTO HEATER = 
LOSS FROM HTR BACKING TO ATM = 
h (HTR BACKING-ATM) = .19 

LOSS FROM HTR PERIMETER TO BED = 
HEAT FLUX THRU PLATE = 
h (PLATE-BED) = 40.18 



123.77 DEG F 
123.71 DEG F 
108.07 DEG F 
229.29 BTU/HR 
2.96 BTU/HR 
BTU/HR-FT-~-2-DEG F 
8.13 BTU/HR 
218.20 BTU/HR 
BTU/HR-FT-'2-DEG F 



LEFT-HAND HEATER 



AVERAGE HEATER TEMPERATURE = 
AVERAGE PLATE SURFACE TEMPERATURE 
AVERAGE BED TEMPERATURE = 
ELECTRICAL ENERGY INTO HEATER = 
LOSS FROM HTR BACKING TO ATM = 
h (HTR BACKING-ATM) = .00 

LOSS FROM HTR PERIMETER TO BED = 
HEAT FLUX THRU PLATE = 
h (PLATE-BED) = 56.11 



121.36 DEG F 
121.29 DEG F 
108.07 DEG G 
262.05 BTU/HR 
.00 BTU/HR 
BTU/HR-FT A 2-DEG F 
4.41 BTU/HR 
257.64 BTU/HR 
BTU/HR-FT"'2-DEG F 



ENERGY BALANCE 



AVERAGE AIR INLET TEMPERATURE 
AIR OUTLET TEMPERATURE = • 
AVERAGE BED TEMPERATURE = 
q (AIR OUT) = 
q ( LOSS F/R WALL) = 
q (TOTAL OUT OF BED) = 
q (TOTAL INTO BED) = 
SUPERFICIAL VELOCITY = 
SUPERFICIAL MASS VELOCITY = 
PARTICLE REYNOLDS NUMBER = 



74.53 DEG F 
102.38 DEG F 
108.07 DEG F 
353.81 BTU/HR 
22.30 BTU/HR 
376.11 BTU/HR 
488.33 BTU/HR 
.30 FT /SEC 
202.78 LBM/HR-FT‘2 
4.05 
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APPENDIX F . 



EQUIPMENT LISTING 



DEVICE 


MANUFACTURER 


MODEL # 


SERIAL n 


Axial Flow 
Turbo Compressor 


Spencer Turbine 


550 SCFM 




Rotometer 


Fisher 5 Porter 


34 SCFM 




30- Tnch 

Water Manometer 


Mer iam 


33KA35 


A69513 


Computer 


Hewlett - Packard 
Company 


HP 85 


2151 A 50182 


Data Acquisition 
System 


Hewlett - Packard 
Company 


HP 3497A 




Dual Flexible Disk 
Drive 


Hewlett - Packard 
Company 


HP 82 90 1M 




Plotter 


Hewlett - Packard 
Company 


HP 7225A 




Printer 


Hewlett - Packard 
Company 


HP 82905B 


214 3J 1296 5 


Copper Const ant an 
Thermocouples 


Omega Engineering 


10 Gauge 




Regulated DC 
Power Supply 
60V, 2 A (Left) 


Lambda 


LK345A 


C1821 


Regulated DC 
Power Supply 
60V, 2 A (Right) 


Hewlett - Packard 
Company 


6296A 


1 5 52A026 18 


Computating 

Bridge 


Rosemount 


9 2 0A 


110 


Calibration 

Bath 


Rosemount 


913A 




Strip Heaters 
250 Watts 


Watlo 


250 Watts 


8044 
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APPENDIX G. 



EXPERIMENTAL UNCERTAINTY ANALYSIS 

In order to arrive at an estimate of the difference 
between measured or calculated values and the true values 
of various parameters, the following uncertainty analysis 
was accomplished: 

1 . TEMPERATURE MEASUREMENTS 

A. Thermocouple Readings - Since the thermocouple out- 
put voltage was measured automatically by the Data Acquisition 
System, human error in reading the output was avoided; however, 
the values arrived at by the System were still subject to 
random error as well as inaccuracies in the empirical corre- 
lation developed during the calibration of the thermocouples. 
Because of these effects, it is believed that the thermo- 
couple readings are accurate to within ±°F. 

B. Ambient Temperature Reading - This measurement was 
accomplished with a mercury- in- glas s thermometer having mini- 
mum graduations of 1°F and a range of -20°F to +140°F. It 

is believed that the thermometer could be read correctly to 
within ±0.5°F. Allowing for the error within the device, 
the uncertainty of the thermometer reading is considered 
to be ±1°F. 
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2. PRESSURE MEASUREMENTS 



A. Bed Pressure Drop - The minimum manometer scale read- 
ing was 0.1 inch with a total range of 30 inches. Readings 
of 0.05 inches were possible. For static measurements, the 
pressure readings are considered to be accurate within 0.2%; 
however, during conditions under which the water column was 
fluctuating, the uncertainty was increased to as much as 2% 
of the full scale reading. 

B. Venturi Inlet Static Pressure - The minimum scale 
reading was 0.1 inch with a total range of 60 inches. Readings 
of 0.05 inches were possible. The pressure readings are con- 
sidered to be accurate to within 0.1% of full scale. 

C. Venturi Differential Pressure - The minimum scale 
reading was 0.001 inches over a range of 2 inches. The 
reading was considered to be accurate to within 0.1% of 
full scale. 

D. Atmospheric Pressure - Minimum scale readings of 
0.001 in. Hg over a range of approximately 32 inches were 
observed. As a consequence, the uncertainty of the reading 
was considered to be 0.1% of full scale. 

3. DISTANCE/SIZE MEASUREMENTS 

A. Bed Width/Depth - Measurements were considered . 
accurate to within ±0.0625 in. or 0.5% of the apparatus 
depth . 
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B. Bed Probe Depth - Due to the bending motion imparted 
to the probes by the particle activity, the depth measure- 
ments were considered accurate to within the 0.25 in. or 
1.5% of the bed depth. 

C. Bed Expansion - Due to the rapidly fluctuating ex- 
panded bed surface level, the bed depth is considered to 
have an uncertainty of ±0.5 in. or 5% of bed depth. 

D. Venturi Inlet/Throat Diameter - These measurements 
are considered accurate to within ±0.001 in. or 0.2® of 
diameter . 

E. Particle Diameter - Determined by optical microscope 
examination conducted by Morgan, [Ref. 2], D^ = 0.0122 (510 
microns) ± 1%. 

F. Plate/Wall Thickness - Measurements were considered 
accurate to within 0.02 in. or 5% of the wall thickness. 

4 . VOLTAGE/ CURRENT MEASUREMENTS 

Minimum ammeter scale graduations were 0.1 Amp with a 
total range of 2 Amps. Because the power supplies were not 
calibrated or compared with a known standard an uncertainty 
of ±10% of the full scale deflection is assigned to the 
current reading. The minimum voltmeter scale reading was 
1 volt with a range of 60 volts. Once again, due to the 
lack of calibration, an uncertainty of 5% is assigned of 
full scale deflection. 
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5. THERMAL CONDUCTIVITIES 



Due to the variation in conductivities found in the 
literature, the following uncertainties are assigned: 

A. k c = 232 Btu/Hr - Ft - °F ±0.51 

B. k mb = 0.112 Btu/Hr-Ft -°F ±10% 

C. k i = 0.025 Btu/Hr - Ft - °F ±10% 

D. k fr = 0.112 But/Hr - Ft - °F ±10% 

6 . DENSITY MEASUREMENTS 

Particle Density - Established by Morgan to be 154.6 
Lbm/Ft^ ±3%. 
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APPENDIX H. ROTOMETER CALIBRATION COMPUTER PROGRAM LISTING 



ROTOMETER CALIBRATION” 



10 PR I N TER IS 701,1 00 
20 5«.S 34 
30 CLEAR* 0 DISP " 

40 DISP “ " 

50 DISP "THIS PROGRAM CONVERTS VENTURI PRESSURE DATA INTO ROTOMETER PLOW RATES’ 1 
60 PAUSE 

70 CLEAR 0 DISP "ENTER PatmtlN Ha ) " 

30 INPUT P 

PO CLEAR e DISP "ENTER Tin (DEG F> " 

100 INPUT T 

110 CLEAR 9 DISP "ENTER Y" 

120 INPUT Y 

130 ‘ CALCULATE Patm(PSI) 

140 Pl»Pt 14. 7/29. 92 

150 ! CALCULATE Tin (DEG R) 

160 Tl-T+459.69 

170 CLEAR 9 DISP "PRESS CONTINUE TO START PROCESSING" 

180 PAUSE 
190 CLEAR 
200 INPUT 
210 CLEAR 



DISP "ENTER FLOWMETER READING” 



DISP 



•ENTER Pig (IN K20) " 

220 INPUT 12 

230 CLEAR 9 DISP "ENTER DELTA P(IN RO) " 
240 INPUT 13 



250 
2 60 



70 P2-I2*14. 7/33. 91/12 



290 



290 P3-P1-P2 



300 



320 

330 

340 



CALCULATE FLOWRATES 
CALCULATE Pig (PSI) 



CALCULATE Pi (PSI) 



CALCULATE DELTA P(PSI) 



310 P4«I3*62. 4*. 834/12/144 



CALCULATE X 
X-P4/P3 

CALCULATE RHOi 



350 D-P3*144/53.34/Tl 

360 ! CALCULATE MASS FLOW (LBM/SEC) 

370 CLEAR 9 DISP "ENTER ASSUMED REYNOLDS NUMBER" 

380 INPUT R 

390 CLEAR 9 DISP "ENTER DISCHARGE COEFFICIENT ASSOCIATED WITH ASSUMED REYNOLDS N 
UMBER" 

400 INPUT C 
410 E« 1.032796 
420 K-C*E 

430 M«K*3. 56327/ 144*Y*S0R <2*32. 174*D*P4* 1 44) 

440 ! CALCULATE VOLUME FLOW RATE (FT3/SEC) 

430 Vl-M/D 

460 ! CALCULATE VOLUME FLOW RATE (FT3/MIN) 

470 V2«V1*60 

480 ! CALCULATE REYNOLDS NUMBER 

490 R2-V1 / 1 4 . 23309* 144*D*4.26/ 12/ . 000000379/32. 174 
500 CLEAR 9 DISP "ASSUMED REYNOLDS NUMBER «",R 
310 « 

320 DISP "CALCULATED REYNOLDS NUMBER -",R2 
330 ! 

540 DISP "DO YOU WISH TO REPEAT CALCULATION USING A DIFFERENT REASSUMED REYNOLDS 
NUMBER? (Y/N)" 

330 INPUT Wl* 

560 IF W1*«"Y" THEN 360 
370 PRINT CHR* (12) 
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580 PRINT 
590 IMAGE 
oOO PF I NT 
© i 0 I MAGE 
c20 PRINT 
©30 IMAGE 
©40 PRINT 
©50 IMAGE 
©60 PRINT 
o70 IMAGE 



USING 5 Q 0 

2/.2SX. " * *ROTOMETEF. CAL I BRAT I ON* * " 

USING ©10 

23X,”A5ME HERSCHEL— YPE VENTURI” 

USING ©30 

33X. "A. 26/2. 13 IN DIA"620 PRINT USING ©30 : F 
USING ©50 ; F 

/.32X.DD, X, "PERCENT READING*' 

USING 670 

32.X . ** ” 



©80 PRINT USING 690 



690 IMAGE 
700 PRINT 
710 IMAGE 
720 PRINT 
730 IMAGE 
740 PRINT 
750 IMAGE 
760 PRINT 
770 IMAGE 
780 PRINT 
790 IMAGE 
800 PRINT 
810 IMAGE 
820 PRINT 



2/.31X, " * ^EXPERIMENTAL DATA* * " 

USING 710 ; P 

2/, 15X, "ATMOSPHERIC PRESSURE , 14 X , DD. DDDDD, X , "IN Hg " 

USING 730 ; T 

15X "VENTURI INLET TEMPERATURE «" , 13X , DD. D. X , "DEG F" 

USING 730 ; 12 

15X, "VENTURI INLET STATIC GAGE PRESSURE - " , 2X , DD . DD, X . " IN H20 
USING 770 ; 13 

15X, "VENTURI DIFFERENTIAL PRESSURE , 3X , D . DDD, X , "IN RED OIL" 
USING 790 : S 

15X, "SPECIFIC GRAVITY OF RED OIL *",16X,D.DDD 
USING 810 

2/ , 31 X , " * ^CALCULATED VALUES**" 

USING 830 ; PI 



830 IMAGE 2/, 15X, "ATMOSPHERIC PRESSURE « " , 1 7X , DD, DDDD , X , "PSI " 

840 PRINT USING 850 ; P2 

850 IMAGE 15X, "VENTURI INLET GAGE PRESSURE * " , 1 OX , D, DDDDD, X , "PSI " 

860 PRINT USING 870 ; P3 

870 IMAGE 15X, "VENTURI INLET ABSOLUTE PRESSURE 5X , DD. DDDDD , X , "PSI " 

880 PRINT USING 890 ; P4 

890 IMAGE 13X, "VENTURI DIFFERENTIAL PRESSURE 8X , D . DDDDD, X , " PSI M 
900 PRINT USING 910 ; D 

910 IMAGE 15X, "VENTURI INLET AIR DENSITY 8X . D. DDDDD, X , "LBM/FT3" 

920 PRINT USING 930 ; X 

930 IMAGE 15X. "VENTURI PRESSURE DROP RATIO (X) 10X , D. DDDDD 

940 PRINT USING 950 ; Y 

950 IMAGE 15X, "EXPANSION FACTOR <Y> -",24X,D.DD 
960 PRINT USING 970 : R 

970 IMAGE 15X, "ASSUMED REYNOLDS NUMBER *",21X,DDDD 
980 PRINT USING 990 ; C 

990 IMAGE 15X, "DISCHARGE COEFFICIENT (C) «",18X,D.DDD 
1000 PRINT USING 1010 ; E 

1010 IMAGE 15X, "VELOCITY OF APPROACH FACTOR (E) 1 OX , D, DDDDD 

1020 PRINT USING 1030 : K 

1030 IMAGE 1 5X , "FLOW COEFFICIENT <K> 21 X, D. DDDDD 

1040 PRINT USING 1050 ; M 

1050 IMAGE 15X « "AIR MASS FLOW RATE «", 1 5X . D. DDDDD , X , "LBM/SEC" 

1060 PRINT USING 1070 ; V2 

1070 IMAGE 1 5X , "AIR VOLUMETRIC FLOW RATE , 14X , DD. DDD. X , "CFM" 

1080 PRINT USING 1090 ; R2 

1090 IMAGE 15X, "VENTURI INLET REYNOLDS NUMBER *"13X,DDDD 

1100 CLEAR e DISP "DO YOU WISH TO PROCESS ANOTHER SET OF READINGS? (Y/N) " 
1110 INPUT W2* 

1120 IF W2*« H Y" THEN 170 
1130 END 
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APPENDIX I. 



THERMOCOUPLE CALIBRATION COMPUTER PROGRAM LISTING 



1 0 
2 0 
3 0 
40 



50 

80 



70 

£0 



SO 
1 0 0 



1 1 0 
1 20 
1 3 0 
1 4 0 
1 5 © 

1 C* 0 

170 
1 £0 



1 90 

2 0 0 
210 



220 
2 3 0 

24 0 
250 

2 8 O 
270 
2 8 0 

2 9 8 

3 0 O 
310 
3 2 0 



340 

358 

380 



n r i p »* ,j 

Dll? " D ft '.V C H L " 

njcp ii H 

013? "THIS CALIBRATION PROGR 
AM IS FOR USE WITH THE HP-34 
97 DATA ACQUISITION SYSTEM" 
DISP " " 

OISP “UP TO 50 THERMOCOUPLES 
MAY BE CAL I BPATED SI MULT AH I 
OUSLY " 

DISP " " 

DISP "UP TO 28 DIFFERENT BAT 
H TEMPERATURE SETTINGS MAY B 
E USED" 

DISP " " 

DISP "THERMOCOUPLE NUMBERS S 
HOULD AGREE WITH CHANNEL HUM 
BERS" 

PAUSE 
DISP “ " 

OPTION BASE 1 

MASS STORAGE IS "=0701" 

S H 0 R T R < 2 0 , 5 6 > , M < 2 0 , 2 > 

DIM T<20,2> 

OISP " " 

DISP "HOW MANY THERMOCOUPLES 
ARE B E I H G CAL I B RATE J ? " 

INPUT K 
D I S P " " 

CUSP "HOW MANY BATH THEMPEPA 
TURE SETTINGS ARE BEING USED 
? " 

INPUT N 
Hi= " AF " 

DISP " " 

DISP "WHAT 
NUMBER?" 

DISP " (. T W 0 
INPUT Y$ 

HiC33=Yi 
C=VAL < Yi > 

HTC33="AL" 

DISP " " 

DISP "WHAT 
UMBER?" 

DISP " (TWO 
INPUT HiC73 
Hi CSG = " VT 1 VR5 " 

Ci = Hi 



IS FIRST CHANNEL 
DIGITS) " 

IS LAST CHANNEL N 
DIGITS) " 
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370 DISP " " 

380 DISP " WHAT DO YOU WISH TO CA 
LL YOUR DATA FILE?" 

390 DISP " < S I X CHARACTERS MAX','" 
40O INPUT DT 
410 FOR 1=1 TO H 
420 DISP " " 

430 DISP " TEMPERATURE SETTING 
# " ; I 

440 DISP " " 

450 DISP "WHAT IS THt PLATINUM 6 
ft T H S E N S 0 R R E S I S TAN C E ? " 

460 INPUT T < I , 1 > 

470 DISP " " 

480 DISP "WHAT IS THE CORRESPOND 
ING BATH TEMPERATURE <C>? U 
490 INPUT J (. 1 1 2':> 

500 DISP " " 

510 DISP "TO START SAMPLING PRES 
S CONTINUE" 

520 PAUSE 

530 clear 709 

540 OUTPUT 709 ; C* 

5 “• 0 F 0 R J = 1 T 0 K 
5 6 O 0 U T P U T 7 0 9 ; "AS" 

=•70 ENTER 709 ; B 

5 8 O R < I , J > = B X 1 0 0 0 
590 NEXT J 

£ l ? u BEEP 
610 } EXT I 
620 DISP 

6 3 O C LEAR C C I S P " S A M P L I N G 

COMPLETE" 

640 DISP " " 

650 DISP "RESULTS HILL BE PRINTE 
D MOMENTARILY" 

6 6 8 DISP " " 

670 DISP "THERMOCOUPLE OUTPUT IS 
IN! MILIUOLTS" 

680 DISP " " 

690 DISP "BATH TEMP. IS IN CELS I U 

11 

700 DISP " " 

710 PRINT "DATA IS STORED IN FIL 
E NAMED " 

720 PRINT USING 730 ; CHR* < 34 > , D 
f , C H R T < 3 4 > 

~30 IMAGE 11 X , A . 6 A .• A 

7 4 CREATE D $ .• 1 > 1 6 5 0 0 
750 ASSIGN* 1 TO Of- 
760 E=C- 1 

770 FOR 1=1 TO K 
7 8 0 FOR • J=1 TO N 
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730 
300 
310 
328 
330 
3 4 0 
350 



M < J , 1) =R < J j I ') 

IK J, 2>=T<' J, 2> 

PRINT# 1 ; fl < J , 1) , M < J , 2 ; 
NEXT J 



F = E+ I 
PRINT 
IMRGE 



USING 350 



X * * 4 : * * * 4 : * * t " 



360 PRINT USING 370 
370 I URGE "********************* 
*********" , /, X 
330 PRINT USING 350 ; F 
850 IMRGE /,6X, "THERMOCOUPLE # " , 
X ^ 2 D 

503 PRINT USING 510 

510 1 M R G E 6 X , " *■ * * *■ * * * * * ^4 * * * * * * * " 

yr y 

520 PRINT USING 538 
530 IMRGE "TEMP" , 4X, "FLRT" , 6X, "T 
C" , 5 X .• " BATH " 

540 PRINT USING 550 
550 IMRGE X, "SET" , 3X, "RESIST" , 3X 
, "OUTPUT" , 3X, "TEMP" 

560 PRINT USING 570 

570 I MRGE " * * * * " , 3X , " ****** " , 3X , . 

" * * * * * X " , 3 X , " * % * '4: " , / 

530 FOR L= 1 TO N 

350 PRINT USING 1000 L,T<L,1>, 

M c‘_ , 1 > , M < L , 2 > 

1 O O O I AGE X , 0 0 .. 4 X .. 0 D 0 0 D 0 , 2 X , D 0 
. DC j , 3X . DO . DO 
1010 NEXT L 

1020 CL ERR 5 DISP " THERMO CO 

UPLE #" , F 

1030 DISP " READY FOR LEAST SQUAR 
ES CURVE FITTING" 

104O DISP " " 

1050 DISP "ENTER O FOR R PONER F 
UNCTION <Y = RX-'B> " 

1 0 60 DISP " " 

107O DISP "ENTER 1 FOR. AN EXPONE 
* NTIRL FUNCTION < Y = fle- BX ) " 
10-80 DISP " " 

1090 DISP "ENTER 2 THRU 10 FOR 
POLYNOMIAL OF N TERMS ■' D 
PEE = N - 1 ) " 

HO© D I S P " < Y =C 1 -r C 2 X * C 3 X 2 + C 4 X ~ 3 
+ . . . . C < N > X y - N - 1 > ) " 

1110 D I M Y < 2 0 > .. X *• 2 O • R < 1 O .. 1 0 > , B ( 

10, 10) , C < 1 0 > , D v 1 0 > 

1120 MAT Y = Z E R C N > 

1.130 MAT X = ZEP < N ) 

1140 FOR J=1 TO N 
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•r i. 



1 150- X< J;'=MOJ, 1 ) 

1160 Y(J)=M<J>2> 

1170 NEXT J 
1180 INPUT Q 

1130 REM "CALCULATE LOGARITHMS 
F X- HND Y- VALUES IF MECES 
fl R Y 11 

1200 IF Q> =2 THEN 1280 



1210 FOR J=1 TO H 

1220 LET Y< J)=LOG<M< J, 2> > 

127.0 NEXT J 

1240 IF Q= 1 THEN 1280 

1250 FOP. J=1 TO N 

1260 LET X < J > =L OG C M < J , 1 > > 

1270 NEXT J 

1280 REM "CALCULATE ELEMENTS OF 
A-MATRIX AND D-VECTOR" 

1230 LET N 1 =Q 

1300 IF Nl; =2 THEN 1326 

1310 LET N 1 =2 

1320 MAT A=2ER<N 1 , Nl ) 

1330 MAT D = 2ER'sNl > 

1340 FOR J=1 TO HI 
1350 FOR K=1 TO Nl 
1360 IF J+K> 2 THEN 1330 
1370 LET A<J,K>*N 
1380 GOTO 1420 
133© FOR L=1 TO N 

14 00 LET A (. J j IO = A (. J / K ^ + X ( L ) ' s C J ■+ 1- 
-2 ) 



1410 NEXT L 

1420 NEXT K 

1430 FOR L= 1 TO N 

1440 IF J>1 THEN 14?0 

1450 LET D< J>=0< J)+Y<L> 

1460 GOTO 1480 

14 70 LET D*: J>=D< J > +Y< L'i .4:X (L > J- 
1 > 



1480 NEXT L 
1430 NEXT J 

1500 REM " SOLVE SIMULTANEOUS LI 
NEAR EQUATIONS" 

1510 MAT B = I N V < A > 

1520 MAT C=B*D 

1530 REM "PRINT EQUATION FOR C'JR 





VE FIT" 








140 


PRINT U 


i“- 


ING 


1.550 


1 5 0 


IMAGE / 


t 


li 


_ ii 


> f. 0 


IF Q>1 


7 


H E N 


1679 


70 


LET Cl* 


£ 


xpa 


XD) 


30 


IF 0=1 


T 


HEN 


1 630 


3 0 


PRINT U 




I H G 


1600 



1680 IMAGE 5 X , " X X X p 0 HER FI I N f. T I fi N 
*** " , . , X 
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Cm O 



1610 PRINT "V=" :Cli "tX^" ;C<2> 
1620 GOTO I860 
1630 PRINT USING 1640 
1640 IMAGE 3X, " ***EXPONENT I AL FU 
NOTION***" , ,X 

1650 PRINT " Y = " ; Cl ; "*EXP< " ; C<2> > 
"*X> " 

1660 GOTO I860 
1676. IF C 1 2 ) > =0 THEN 1740 
1630 PRINT USING 1630 
1636 IMAGE 3X , " ***POL YNOM I AL FUN 
CT I ON*** " 

1780 PRINT USING 1710 ; Q-l 
1710 IMAGE 10X, "< ORDER = "..0>" )" > 
/ , X 



” v =" ,C(1); 0(2 






1720 PRINT v 
1730 GOTO 1730 
1740 PRINT USING 1750 
1750 IMAGE 3 X .• " * * * POL V H 0 M I A L FUN 
CT I ON*** " 

1760 PRINT USING 1778 ; Q-l 
1770 IMAGE 1 OX , " < ORDER = u ,0," , 

/ , X 

1730 PRINT " Y=" 0(1 );"+"; 0(2) ; "* 
X u 



1730 


IF 0 = 2 


THEN I860 




1 3 0 0 


FOR J = 3 


TO 0 




1 3 1 0 


IF C(J> 


> = 0 THEN 1 


o 


1320 


PRINT 0 


< J) ; " *>:'-" ; 


J 


1330 


LOT 0 13 


50 




1340 


PRINT " 


+ " 0 < J ) ;• " * 


X 


1 85=10 


NEXT J 






1 360 


REM •'PRINT INPUT 


ij 




X AND V 


AND CALCU 


L 




ES OFY " 






1370 


IF Q> =2 


THEN 1350 




1330 


F O.R J = 1 


TO N 




1 330 


LET Y< J)=EXP(Y( J) 


> 


1 3 0 0 


NEXT J 






1310 


IF 0=1 


THEN 1350 




1320 


FOR J= 1 


TO N 




1330 


LET X<J 


)=EXP(X( J) 


) 


1 3 4 0 


NEXT J 






1350 


PRINT 






136 O 


PRINT !J 


S I N ij 1370 




13 70 


IMAGE 


, 4 X , "X" , 6X 


i 




• 5X , " Y 


(CALC) " 




1330 


PRINT lj: 


SING 133Q 




1 3 3 0 


IMAGE 3! 


X, " * * * " . 5X 


i 




, 5X, "**: 


* ***** ° 




2 0 0 0 


LET 3=0 






2 O 1 0 


FOR J=1 


TO N 




2020 


IF 0> =2 


THEN 2080 




2 O 3 0 


IF 0=1 ' 


THEN 2-0 6 0 




204 O 


LET Y 1 = i 


: i *x < j ) -c ( ; 


r, 



J-l ; 



.'ALU 



V (ACT)" 
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2050 GOTO 2120 

2O60 LET Y1=C1*EXP<C<2>*X< J> > 
2070 GOTO 2120 
2080 LET Y 1 = C ( 1 ) 

2020 FOR K=2 TO Q 

2100 LET Y 1 = Y 1 + C C K> t X < J > A < K - 1 > 

2110 NEXT K 

2120 LET S=S+CYCJ)-Y1 > •'' 2 
2130 PRINT USING 2140 ; XCJ>,YCJ 
) , Y 1 

2140 I MAGE X ■ DD .ODD, 5X .• DO . 00 , 6X , 
DO ODOn 
2150 NEXT J 
2160 PRINT 

2170 PRINT "SUM OF SQUARE ERRORS 
= “ ; S 

2180 CLEAR G DISP "DO YOU DESIRE 
A PLOT?" 

2180 DISP "CENTER Y OR N>° 

220O INPUT PS 

2210 IF P $ = " N " THEN 3680 

2220 GOSUB 2258 

2238 PAUSE 

2240 GOTO 3620 

2250 DIM Y1 C4> , A1 C5> , Ml <5> 

2260 Y K1),Y1<3) = INF G Y 1 C 2 > , Y 1 C 
4 > = - I N F 6 A1 < 1 ) , fll <2> / hi (3> 
, A1 C4X, A1 < 5 > = 0 



2 2 i** 0 


F L> R 


J 


= 1 


TO 


N 










2280 


IF 

i \ 


Y 1 


c 2 ; 


> <x c 


■J> 


THEN 


Y 1 


< 2 > 


X 

it 


2220 


.J } 

IF 


Y 1 


< i : 


> y x c 


J> 


THEN 


Y 1 


( 1 > 


= X C 


23O0 


• J ? 

IF 


Y 1 


C4: 


> <YC 


J) 


THEN 


VI 


S 4 > 


= V < 




J> 


















2310 


IF 


Y 1 


< 3 ; 


> > Y C 


J> 


THEN 


Y 1 


r* 3 '* 


it 

h: 

/-N 




J> 


















2320 


A 1 C 


1 > ■ 


= A 1 


C 1 ) 


+ XC 


J) G 


A 1 


i*’ 2 


= h 1 




2 > 


+ X 1 


C J : :■ 


t X c 




G A 1 C 


3 *•' : 


=fll 


c 3 :■ 



+ Y < J > 

2330 A 1 C 4 > = A 1 C 4 > +Y C J > 1: Y C J G A1C 
5 > = A 1 <5>+XC J>*Y< J>. 

2340 NEXT J 

2350 M 1 < 1 > = A 1 < 1 ) / N G M1C2) = CA1C2 
: , -Hl(l)|:ftl<l)/H)/(fi-l> 6 Ml 
C 3 > = A 1 C 3 > ■••'N 

2360 Ml C 4 > = C A 1 C 4 >-A 1 C3>*A1 <3>/N) 
/ C N - 1 > G Ml < 5 > = C A 1 < 5 > - A 1 C 1 > 
*• A 1 C 3 > N > / Hi C 2 > - A 1 C 1 ,■ t A 1 < 1 
> s N ) 

2370 R 5 = M 1 C 3 > -M 1 C 1 > % M 1 < 5 > G 05 = 0 

2380 CLEAR G D3= 1 G DISP "AUTO X 
-SCALING? = Y / N " 
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ro to ror o ro ro rororo ro ro roro ro roro ro r-« ro roro ro roro roro ro rororo 

CT*» CTi »7i Ch CT» Cf> Ch Ch i7n CO CO <71 CO cH lH CH '01 • >1 4* 4 * 4* 4- 4^ 4^ 4* OJ 

ro ch lh 4^ i\j ro ♦-* ci 'vO •:•:• - j ch cn -u »\i ro ^ •:*:• «x» - j ch c/i 4^ oj ro ►— c» ko 

d» d> d 1 d* d* d d' d' d* d» d‘ d* d d' d d' © d' ‘*0 O d 1 d> © Q O O © d d> d 



INPUT F:J 

ON FHR GOTO 23S0 , 24 1 0 , 2420 

x i =v i < i > o ;: 2 =yi < 2 ; 0 goto 

2470 

CLEAR & DISP "ENTER MINI NUN 
VALUE FOR X-AXiS" 

INPUT XI 

CLEAR 0 DISP "ENTER MAXIMUM 
VALUE FOR' X-AXIS" . 

INPUT X2 

IF X 1 > = X 2 THEN BEEP 0 GOTO 
2 420 

CLEAR 0 DISP "DO YOU NISH V 
ERTICAL OR HORIZONTAL LABEL 
S ON THE X-AXIS? = V-H" 

INPUT RT 

IF NOT LENCRT> THE!- BEEP 0 
G 0 TO 2 4 7 0 

S.- = ! S IF u z C f ’■ R $ C 1 > 1 3 > = " V " 
T H z. N ’i< 0—0 o 1 j 0 T U 3 Zt 2! 0 
IF UPcitipici, 13>#"H" THEN B 
EEP 0 GOTO 2470 
CLEAR 0 DISP "DESIRED NUMBE 
R OF INTERVALS ON X-AXIS? < 
<=16> " 

INPUT LI 

I F *■ L 1 < 1 OR L 1 > 1 6 OR LI # I NT C 
LI I- THEN BEEP 0 GOTO 2520 
CLEAR 0 DISP "DESIRED NUMBE 
R OF SUBDIVISIONS BETWEEN L 
ABELED INTERVALS?" 

INPUT L3 

SS=0 6 IF L3=0 THEN SB=1 0 
L3= 1 

IF L3 < 1 OR L3> L 1 THEN BEEP 
0 GOTO 2550 

CLEAR 0 DISP "DO YOU DESIRE 
AUTOMATIC SCALING OF THE Y 
-AXIS? = Y/N" 

INPUT R t 

0 N F N R G 0 T 0 2 5 ? O ,2678, 2 6 2 O 
CLEAR 0 DISP “ENTER MINIMUM 
DESIRED VALUE FOR Y-AXIS" 
TNPUT Y 1 

CLEAR 0 DISP "ENTER MAXIMUM 
DESIRED VALUE FOR Y-AXIS" 
INPUT Y2 

IF Y 1 > Y 2 THEN BEEP 0 GOTO 2 

620 ELSE 26S0 

Y-l = Y 1 3 >. 0 Y2 = Y 1 4 > 

CLEAR 0 DISP, "ENTER DESIRED 
NUMBER* OF INTERVALS ON THE 
Y-AXIS <<=1 2 1'" - 
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•>j oi ojco rorororo to r\>f\> ro roro roro rv roruro roro roro roro ro r o 10 ro to to ro to to 

cd cd cd cd 'iTt 'iii •£» •.£* *D •‘O 'ji* '£• *-o *o co co cococo coco coco ~j ~%j -*.j -sj -j-.j -nj *xjc^ 

•>j ro *-*■ cd ».£• co "sj Ot fj] -t± gj ro »— * cd *x* co ~ j Or cn x* oi ro *-•> cd *o co ~j *.n -u oj ro ►-* cd •£• 

CD 'S< CDO »2D »ID »ID *3D CD CD 'ID CD CD CD 'ID 'HZ' O *ID CD CD CD CD 'ID CD 'X' O CD 'ID 'ID 'ID CD 'ID CD CD 



INPUT L2 

IF L2 < 1 OR L2> 1 2 OP L2#INT< 
L2> THEN BEEP G GOTO 268Q 
CLEAR G DISP "ENTER DESIRED 
NUMBER OF SUBDIVISIONS BET 
WEEN LABELED INTERVALS" . 
INPUT L4 

S7=0 6 IF L4=0 THEN S7=l G 
L 4= 1 

IF L4 < 1 OR L4> 1 2 THEN BEEP 
G GOTO 2710 

Z1=X2-X1 G Z3»IHT<£00/L1>*L 
1 G Z2=Y2~Y 1 G Z4=INT(144/L 
2>*L2 

D1=ZKZ3 G D 2 = Z 2 y Z 4 G Z5 = K2 
+ < 2 0 7 - Z 3 > * 0 1 G X0=X 1-43*01 
S Y0=Y 1 -40*02 

GCL-EAR G SCALE NO , Z5, Y0, Y2+ 
< 151-Z4>*02 

X H X IS Y 1 > Z 1 / L 1 > X 1 .• X2 + D 1 

YAK IS X 1 .• Z2/L2* Y1 , Y2 + D2 

IF S3 THEN 2968 

W 1 =LGT < ABS < Z 1 /L 1 *L3 ) > 

W=LGT(ABS(X1+<X1=0))>+1 

IF U > 5 - <. S G N C X 1 ) = - 1 > OR Wl<- 

3 THEN 2920 

IF S9 THEN 3520 

LOIR 90 

FOR J = X 1 TO X2+D1 STEP INI 
1*L3 

MOVE j-4*01 , Y O 
GO SUB 3320 • 

LABEL V £ C 1 > V 3 
NEXT J 
GOTO 2960 
MOVE NO, Y0 
LDIR O 

I1=X1 G GOSUB 3450 G I 1=ZK 
LI G Z* = Vf G GOSUB 3450 
LABEL "KMIN="&ZTC1, J03^" =TI 
CS=“&V*C1 , 103 
IF S7 THEN 3140 
Wl=LGT<fiBSvZ2/L2*L'4>> 

H = LGT CABS < Y 1 + < Y 1 =0 > ) ) + 1 
IF U> 5- < SGN < Y 1 > =- 1 > OR W 1 <- 
3 THEN 3O70 
LDIR O 

FOR J = Y 1 TO Y2 + D2 STEP Z2/L 
2*L4 

MOVE XO > J-4 *02 
GOSUB 3320 
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304;> LABEL VTC1,V3 

3 0 5 C NEXT J 

306C GOTO 3148 

3070 LOIR 90 

38S0 I 1 = Y 1 6 GOSUB 3450 

3 0 9 0 M 0 V E X 0+12 1 0 1 , Y 1 

3 1 ©0 LABEL " YM I N = " &V* C 1 ,103 

3119 MOVE X G +24*01, Y1 

7120 1 1 = 2 2 , ' L 2 6 GOSUB 3450 

3130 LABEL " T I CS= " & V* C 1 ,103 

3140 PENUP G LDIR 0 

3150 MOVE X 1 ) , Y ( 1 > 

3160 FOR J=1 TO N G MOVE X < J ) 
D 1 , Y < J ) -4*02 (2 LABEL " * " 
NEXT J G BEEF' 

3170 FOR L=0 TO 100 
3 1 SO X 4 = Y 1(1) + L* < < Y 1 < £ > -Y 1(1) 
00 ) 

3190 IF 0 > = 2 THEN 3250 
320O IF Q= 1 THEN 3230 
3210 LET Y4=C1*X4-C<2> 

3220 GOTO 3290 

3230 LET Y4=C1*EXP<C<2)*X4> 

3240 GOTO 3290 

3250 LET Y4=C < 1 ) 

3260 FOR K=2 TO Q 

3270 LET V 4 = Y 4 + C < K > * X 4 < K - 1 > 

3280 NEXT K 

3290 PLOT X 4 .. Y 4 

330O NEXT L 

3310 BEEP 6 G 0 T 0 339 O 

3320 VT = 11 " 6 X = J 

3330 V f = V A L T < X ) 

3340 tF POS < VS , " E " ) THEN 3408 
3350 G9 =L GT < ABS < X+ < X = 8 ) ) > 

3360 IF LEN<V$)>5 AND ABS<G9> 
<SGN<X>«-1> THEN V=5 G V 

» Q RETURN 

3370 IF LEN < Vf ) <5 THEN V=LEN< 
0 RETURN 

3386 Vf Cl , 53=VRL$ <X) G V = 5 
3390 GRAPH G RETURN 
34O0 E6=P0S<Vf , "E" X 
3410 IF V f C 1 , 10 = ”- " THEN V4C3 
$■ C E 0 3 G G 0 T 0 34 3 O 
3420 V T C 2 3 = V 4 L E 0 5 
34 36 V = L E N < V T > G IF V>5 THEN 
NT USING 3560 ; J G V = 5 

_ H M 

3440 RETURN 
3450 Vf="" 

3460 Vf = VRLT-< II) 



-24 

G 

) ✓ 1 



4 - 

.j ii 

VT) 

3=V 

F'RI 
G V 
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3478 IF POSCVS, "E" > THEN 3500 
3480 v$ C 1 , 1 0T = VAL$ < I 1 > 

3490 -RETURN 

33O0 Vf C6, 1 03=V$CPOS < 0# .. M E M >3 
3510 RETURN 
3520 LOIR 0 <2 L 9 = - I H F 
3530 FOR J = X 1 TO X 2 + D i STEP 21-"L 
14L3 

3540' ij 0 S U 0 3320 

3550 IF L 9 > J - < V * 4 + 6 > 4 0 1 OR L9>Z5 
+ 0-04:3 >4:01 THEN PRINT US IN 
G 3560 i? GOTO 3610 
3568 IMAGE "LABEL DELETED AT “ ,7 
D . 4D 

3570 MOVE J+ C 2 - V 4 4 > 40 1 , V 1 - 1 24D2 
3 5 3 0 L 9 = J + < V44+2 > 4 D 1 
3590 IF L9‘> Z5 THEN MOVE Z5+C2-V* 
8 > 4 D 1 j Y1-124D2 
3600 LABEL V4C1.V3 
3610 NEXT J 
3620 GOTO 2960 
3630 DEF FNR 

364G IF NOT LEN < R$ > THEN 11 = 1 6 
GOTO 3660 

3650 I 1 = P 0 S < "YN" , IJPC$<R*C1 . 13> > + 
[ 

3660 IF 11 = 1 THEN BEEF- 
3670 F N R = I 1 
3630 FN END 

3690 CLEAR 3 DISP "DO YOU WISH T 
0 FIT A DIFFERENT TYPE CURV 
E ? < Y / N > " 

370O INPUT US 

3710 IF W$="Y M THEN 3730 

3720 GOTO 3740 

3730 GOTO 1020 

3740 'NEXT I 

3750 ASSIGN# 1 TO 4 

3760 BEEP 

3770 CLEAR 3 DISP " CAL I BRAT 

ION COMPLETE” 
w' t 8 0 END 
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APPENDIX J 



HEAT TRANSFER DATA ACQUISITION COMPUTER PROGRAM LISTING 



10 MASS STORAGE IS M : D 7 C>1 ” 

SHORT VC5G> 

TO DIM C<3> ,7(80) ,R(S) .L(3>.Q!*C1003 
40 CIS*" “ 

CO CLEAR * DISP " QUICK" 

DISP " " 

70 DISP "THIS PROGRAM IS FOR USE IN SAMPLING FOR FUTURE PROCESSING " 

SO DISP "OPERATING PARAMETERS AND TC OUTPUTS ARE STORED ON DISK FOR FUTURE REDUC 
T I ON" 

PO PAUSE 

100 CLEAR 0 DISP "WHAT IS THE DATE? (MO/DA/YR) " 

110 INPUT D1S 

ICO CLEAR 0 DISP "WHAT IS THE TIME? (MILITARY) “ 

130 INPUT Tl* 

140 CLEAR 0 DISP "WHAT IS THE BED WIDTH? (INCHES)" 

130 INPUT X 

160 CLEAR 0 DISP "WHAT IS THE STATIC BED HIEGHT? (INCHES)" 

170 INPUT SIS 

ISO CLEAR 0 DISP "WHAT WAS THE HIGHEST NUMBERED PROBE UTILIZED? (TWO DIGITS)" 

190 INPUT LS 
200 L-YAL(LS) 



210 

220 


CLEAR 

INPUT 


0 DISP 
VI 


"WHAT 


IS 


THE 


VOLTAGE 


TO THE 


RIGHT-HAND HEATER? 


(VOLTS) 


230 

240 


CLEAR 

INPUT 


0 DISP 
11 


"WHAT 


IS 


THE 


CURRENT 


TO THE 


RIGHT-HAND HEATER? 


(AMPS) " 


250 

260 


CLEAR 

INPUT 


0 DISP 
V2 


"WHAT 


IS 


THE 


VOLTAGE 


TO THE 


LEFT-HAND HEATER? 


(VOLTS) " 


270 

280 


CLEAR 

INPUT 


0 DISP 
12 


"WHAT 


IS 


THE 


CURRENT 


TO THE 


LEFT-HAND HEATER? 


(AMPS)" 


290 

300 


CLEAR 

INPUT 


0 DISP 
A 1 


" WHAT 


IS 


THE 


AMBIENT 


TEMPERATURE? (DEG F) " 




310 

320 


CLEAR 

INPUT 


0 DISP 
FI 


"WHAT 


IS 


THE 


AIR FLOW-RATE? 


(7.) " 





330 F2-F1 1.373 

340 CLEAR 0 DISP "WHAT IS THE BED EXPANSION? (INCHES)" 

330 INPUT El S 

360 CLEAR 0 DISP "ENTER ANY COMMENTS YOU WISH RECORDED (MAX 100 CHARACTERS)" 

370 INPUT 01 SC 13 

380 CLEAR 0 DISP "WHAT DO YOU WISH TO CALL YOUR DATA FILE? (SIX CHARACTERS MAX) 1 
390 INPUT FS 
400 HS» " AF01 AL " 

410 HSX73- M 79" 

420 HSC93-"VT1 VR3TE0" 

430 CLEAR 709 

440 OUTPUT 709 ; HS 

450 CREATE FS, 1,1300 

460 ASSIGN# 1 TO FS 

470 ASSIGN# 2 TO " C0EFF2 : D700 " 

480 PRINT# 1 } D1S,T1*, X,S1S,LS,L, VI, II, V2, 12, A1 , Fl , F2, El S , 01 S 
490 PRINT USING 500 

500 IMAGE H tmtttttxttttmttttttttnttt n 

310 PRINT USING 320 

320 IMAGE 3/ , 3X, " **EX PER I MENTAL DATA**" 

330 PRINT USING 340 ; FS 
540 IMAGE 8X , "RUN #",X,6A 
330 PRINT USING 360 j D1S,T1S 
560 IMAGE 6X,BA, X, X,4A 
370 PRINT USING 380 

380 IMAGE 2/,2X, "TC#",3X, " OUTPUT ",3X, "TEMP" 
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z~0 PP I NT US I N3 600 

rv_- IMAGE ?X , " < mV; " . 5X , " . DEG F^'lISTaoO 
c 1 < PRINT USING ©20 

©20 IMAGE 2X « " M , 3X . " " . 5> . " " 

©30 CLEAR »•' DISP “PRESS CONTINUE WHEN READY TO BEGIN SAMPLING" 

©40 PAUSE 

©50 OUTPUT 70° : ” TEC " 

©oO FOP 1=1 TO 7? 

©70 OUTPUT 70? ; " AS" 

680 ENTER 70? ; R 
690 V(I)-R«1000 
700 NEXT I 

710 OUTPUT 709 ; " TE1 TE" 

720 ENTER 709 ; S 
730 BEEP 

740 CLEAR <? DISP "SAMPLING COMPLETE" 

750 DISP "RESULTS WILL BE PRINTED MOMENTARILY" 

760 READ# 2, 1 
770 FOR 1*1 TO 79 

780 READ# 2 ; C < 1 > , C <2> , C (3) , C <4> , C <5> 

790 IF I >L THEN 810 

800 GOTO 320 

810 IF I <*67 THEN 940 

820 T (I ) -C < 1 ) +C<2) *V< I >+C (3) *V ( I ) *V(I> +C (4) *V < I ) *V ( I ) *V < I ) +C <5> *V ( I ) *V (I ) *V < I ) 
(I) 

830 T ( I ) — T < I ) f 9/5+32 
840 PRINT# 1 ; V(I> , T(I> 

850 IF I #41 THEN 920 
860 PRINT USING 870 

870 IMAGE 3/,2X, "TC#'\3X, "OUTPUT", 5X, "TEMP" 

880 PRINT USING 890 

890 IMAGE 9X , " (mV) " , 5X , " <DEG F) " 

900 PRINT USING 910 

910 IMAGE 2X," " , 3X, " ",5X," " 

920 PRINT USING 930 ; I , V < I ) , T < I > 

930 IMAGE 2X , DD , 3X , DD. DDD ,4X, DDD , DDDD 

940 NEXT I 

950 ASSIGN# 1 TO * 

960 ASSIGN# 2 TO * 

970 R<1)*<T<1)+T(2>+T<3) )/3 
980 R < 2 ) * < T ( 5 ) +T ( 6 ) +T ( 7 ) +T ( 8 ) ) / 4 
990 R(3)*(T<9)+T(10)+T<11)+T(12) >/4 
1000 R ( 4 ) * ( T (13) +T (14) +T (15) +T (16) > /4 
1010 R(5>-<T<17)+T<18)+T<19)+T<20> >/4 
1020 L < 1 ) — <T <21 > +T (22) +T (23) +T (24 ) ) /4 
1030 L ( 2 ) = < T ( 25 ) +T ( 26 ) +T ( 27 ) +T ( 28 ) ) /4 
1040 L (3) * (T (29) +T (30) +T (31 ) +T (32) ) /4 
1050 L (4) - <T (33) +T (34) +T (35) +T (36 ) ) /A 
1060 L ( 5 ) * ( T ( 37 > +T ( 38 ) +T ( 39 ) +T ( 40 ) ) / 4 
1070 G1-(T(41)+T(42)+T(43) )/3 
1080 G2-(T(44)+T(45)+T(46 ) ) /Z 
1090 PI* (T (47) +T (48) +T (49) ) /3 
1100 Bl-0 
1110 Xl-0 

1120 FOR 1-30 TO L 

1130 IF 1-55 THEN 1130 

1140 GOTO 1170 

1150 Xl-1 

1160 GOTO 1260 

1170 IF 1-65 THEN 1190 
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; 1 so 

: : =>o 

z 1 " 

1 zee 

1240 

i 250 

1260 

1270 

1280 

1290 

1 300 

1310 

1 320 

1 330 

1 340 

1350 

1 3©0 

1370 

1380 

1390 

1 400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

6160 



GOTO 1210 
GOTO 1 2©0 

I” i=oo then izgo 
GOTO 1250 
x 1=3- 

GGTC 1260 
B1=B1*T ( I > 

NEXT I 

N1-L-50+1-X1 

Hl = <R ( 1 ) +R (2) +R <3) +R (4) +R (5) ) /5 

H2= <L ( 1 ) +L <2) +L (3) +L (4) +L <5) ) /5 

BI=B1/N1 

E 1 = V 1*11*3.4121 

E2*V2*I2*3. 4121 

IF XI >1 THEN 1360 

P2=T (55) 

GOTO 1400 
IF XI >2 THEN 1390 
P2= <T (55) + T (65) ) /2 
GOTO 1400 

P2= (T (55) +T < 65) +T (66) ) /3 
PRINT USING 1410 : X 

IMAGE 3/, "BED WIDTH = ",DD.DD, X. "IN" 

PRINT USING 1430 ; SI* 

IMAGE "STATIC BED HEIGHT = " , 5A, " IN" 

PRINT USING 1450 ; F2 

IMAGE "AIR FLOW RATE * " , DDD. DD, " CFM" 

PRINT USING 1470 ; A1 

IMAGE "AMBIENT TEMPERATURE- " , DD . DD , " DEG F" 
PRINT USING 1490 ; VI 

IMAGE "RIGHT HEATER VOLTAGE* " , DD. D, " VOLTS" 
PRINT USING 1510 ; II 

IMAGE "RIGHT HEATER CURRENT = " , D. DD, " AMPS" 
PRINT USING 1530 ; V2 

IMAGE "LEFT HEATER VOLTAGE * " . DD . D , " VOLTS" 
PRINT USING 1550 ; 12 

IMAGE "LEFT HEATER CURRENT = " , D. DD, " AMPS" 
PRINT USING 1570 ; El* 

IMAGE "BED EXPANSION « ",4A, " INCHES" 

PRINT USING 1590 ; 01* 

IMAGE "COMMENTS: ", 100A 

END 
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appendix k 



HEAT TRANSFER DATA REDUCTION COMPUTER PROGRAM LISTING 



10 SHORT V(30> 

20 PRINTER IS 701,100 

21 PRINT CHR* (2*7) 1L" 

30 p R I NT CHR* (27) VU 6D" 

40 Din C<5) ,T(80) ,R(5> ,LC5> .OISCIOCO 
50 01 ” 



60 E p «0 



70 P9»0 

80 CLEAR t? DI5P " PROCESS" 

90 DISP " " 

100 DISP "THIS PROGRAM IS FOR USE IN PROCESSING DATA WHICH HAS BEEN PREVIOUSLY R 
ECORDED. " 

110 DISP " " 

120 DISP "HEAT TRANSFER COEFFICIENTS ARE CALCULATED, DISPLAYED, AND PRINTED." 

130 PAUSE 

140 CLEAR 0 DISP "WHAT IS THE NAME OF THE DATA FILE? (SIX CHARACTERS MAX)" 

150 INPUT F* 



160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

290 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 



CLEAR 0 DISP "PRESS CONTINUE TO BEGIN PROCESSING" 

PAUSE 

ASSIGN# 1 TO F* 

READ# 1 ; D1*,T1*, X.S1*.L*,L, VI, II, V2. 12, Al , FI , F2, El *, 01* 

F2-F 1 * . 324 54+ .795155 

CLEAR 0 DISP "DO YOU WISH A PRINT OUT OF THE EXPERIMENTAL DATA? (Y/N) " 
INPUT RB* 

IF R8t«"Y" THEN 260 
B9-1 

GOTO 1110 

CLEAR 0 DISP " COMMENT STATEMENT: " 

DISP " " 

DISP 01* 

DISP " " 

DISP "DO YOU WISH TO MAKE ANY CORRECTIONS? (Y/N)" 

INPUT R9* 

IF R9*«"N" THEN 360 
DISP " " 

DISP "ENTER CORRECTED COMMENT STATEMENT- 
INPUT 01* 

i 

i 

! x* f PRINT EXPERIMENTAL CONDITIONS*** 

i 

PRINT CHR* (12) 

PRINT USING 420 

IMAGE 3/,30X, " **EX PER I MENTAL DATA**" 

PRINT USING 440 ; F* 

IMAGE 35 X , "RUN #",X,6A 
PRINT USING 460 j D1*,T1* 

IMAGE 33X , 8A, X, , X , 4A 
IF E9*l THEN 1030 
PRINT USING 490 ? F* 

IMAGE 5/ , 22X , "RUN NUMBER ~",20X,6A 
PRINT USING 510 ; X 

IMAGE / , 22X , "BED WIDTH 14X, DD. DD, X ," INCHES" 

PRINT USING 530 ; SI* 

IMAGE /,22X, "STATIC BED HEIGHT 6X , 5A, X ," INCHES" 

PRINT USING 550 ; F2 

IMAGE /,22X,"AIR FLOW RATE , 12X, DDD.DD, X, "CFM" 

PRINT USING 570 > Al * 
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IMAGE 


. 22>. "AMBIEN- "Emp£ 


RAT JR- = '■ , 5X , D2 . DD . X . 


h 

irj 

id 

n 


59 v 


=•= : n7 


USING 5®0 : VI 








I MAGE 


. 22X . "RIGHT HEATER 


V0_TA3E = " . 5X . DD. D . X 


. " V0_ S 


=•00 


=*r:n~ 


US I NG slO : 11 






s 1 C' 


I MAGE 


. 22X . " R I GHT HEATER 


CURRENT = 1 , © X . D . DD . X 


. "AMPS" 


-20 


print 


-SING ©30 : V2 






r 30 


: mage 


/ . 22X . "LEF T HEADER 


VOLTAGE . ©X . CD. D« X . 


"VO^TS" 


-40 


PRINT 


USING 650 : 12 






©50 


IMAGE 


/ . 22 X . "LEFT HEATER 


CURRENT =".7X,D.DD. X. 


■AMPS" 


©©0 


PRINT 


USING ©70 : El* 






670 


IMAGE 


/ . 22 X . " BED EXPANSION 1 1 X , 4A. X. " INCHES 


" 


680 


29=LEN (01* ) 






690 


c:scD9^i . ioo:=" ■■ 






700 


IF D 9 


30 THEN 740 






710 


PRINT 


USING 720 s 01*11.293 




720 


IMAGE 


/. 22X . "COMMENTS: " , X 


, COA 





730 3GTC 970 
740 K9*25 

750 IF C2*:K9,K'9D* m " THEN 790 
7»0 k'?«K9-i 
770 GOTO 750 

730 PRINT USING 790 : Q1SC1.K93 

790 IMAGE /.22X. "COMMENTS: X.40A 

900 D2-K9-35 

910 IF D9>D2 THEN B70 

920 KB-K9*1 

SCO IF K8>D9 THEN 970 

940 PRINT USING 850 ; 01 *CK9-*-l . D93 

950 IMAGE 22X.50A 

960 GOTO 970 

970 L9=C2 

SSO IF 01«CL9,L93«" " THEN 910 
890 L9-L9-1 
900 GOTO 880 

910 PRIN~ USING 920 : Q1SCK9«-1 . L93 
920 IMAGE 22X.50A 
930 LB-L9-1 

940 IF LS>D9 THEN D9»L8 

950 PRINT USING 960 : Gl*CL9+l.D93 

960 IMAGE 22 X , 50 A 

970 PRINT CHRS < 12) 

980 E9«l 
°90 GOTO 360 
1000 ! 

1010 ! * PRINT THERMOCOUPLE READINGS**! 

1020 1 

1030 PRINT USING 1040 

1040 IMAGE /. COX, "THERMOCOUPLE READINGS" 

1050 PRINT USING 1060 

1060 IMAGE /, 13X, "TC#" . 3X . "OUTPUT" . 4X, “TEMP", 14X, "TC#’\3X, "OUTPUT ", 4X , "TEMP" 
1070 PRINT USING 1080 

1080 IMAGE 20X, M (mV) " , 4X. " (DEG F ) ", 1 9 X , " (mV) " , 4 X , " C DEG F) " 

1090 PRINT USING 1100 

1 lOO IMAGE 13X. " ** . 3X , " ",4X, ” ", 14X. " M ,3X. " ".4X, " " 

1110 FOR 1-1 TO 79 
1120 IF I >L THEN 1140 
1130 GOTO 1150 
1140 IF K-67 THEN 1160 
1150 READ# 1 : V<I> ,T(I) 

1160 NEXT I 
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12 ”0 IF B*=l THEN 13^0 

: i so w i *0 
: i qcj z i =o 

: 300 “Or. f =1 TO 40 

1010 IF W 1 = 1 THEN J*J + 1 6' GOTO 1030 

1020 J*K+40 

1030 IF J >L THEN 1050 

1040 GOTO 1 060 

1050 IF J =67 THEN 1340 

1260 IF J >79 THEN Zl*l 

1070 IF Z 1 =1 THEN 1310 

1080 PRINT USING 1290 s K. V <K) . T <K) , J . V < J ) , T ( J) 

1090 IMAGE 1 3X . DD . 3X , DD . DDD, 4X , DDD . DD, 1 3X . DD. 3X. DD. DDD, 4X, DDD. DD 
1300 GOTO 1330 

1310 PRINT USING 1320 ; K.VOO.TOO 
1320 IMAGE 13X.DD.3X, DD. DDD , 4X , DDD. DD 
1330 GOTO 1360 



1340 J-J+l & Wl»l 



1350 GOTO 1250 

1360 NEXT K 

1370 ASSIGN# 1 TO * 

1380 PRINT CHRS (12) 

1390 R ( 1 ) * (T ( 1 ) ***T (2) +T (3) ) /3 
1400 R ( 2) = ( T ( 5) +T ( 6 ) +T ( 7) +T (8))/4 
1410 R(3)«(T<9) +T< 10>+T (1 1 )+T (12) ) /4 
1420 R<4)«(T(13)+T(14)+T(lS)+T(16) ) /4 
1430 R ( 5) = ( T ( 1 7 ) +T ( 1 8 ) +T ( 1 9) +T ( 20 ) ) / 4 
1440 L ( 1 ) c < T (21) +T ( 22 ) +T (23) +T (24) ) /4 
1450 L (2) =* (T (25) +T (26) +T ( 27) +T (28 ) ) /4 
1460 L ( 3) = (T ( 29 ) +T (30) +T (31 ) +-T (32 ) ) /4 
1470 L ( 4 ) * (T ( 33) +T (34) +T ( 35) +T (36) ) /4 
1480 L (5) * (T (37) +T (38) +T (39) +T (40) ) /4 
1490 Hl«(R(l)+R(2)-*-R(3)+R(4)+R(5) ) /5 
1500 H2*(L(1)+L(2)+L(3)+L(4)+L(5) ) /5 
1510 IF B9-1 THEN 3500 
1520 ! 

1530 ! ***DISPLAY HEATER TEMP SNAPSHOT (R-H) *** 

1540 ! 



1550 


PRINT 


1560 


IMAGE 


1570 


PRINT 


1580 


IMAGE 


1590 


PRINT 


1600 


IMAGE 


1610 


PRINT 


1620 


IMAGE 


1630 


PRINT 


1640 


IMAGE 


1650 


PRINT 


1660 


IMAGE 


1670 


PRINT 


1680 


IMAGE 


1690 


PRINT 


1700 


IMAGE 


DDD. DD 


1710 


PRINT 


1720 


IMAGE 


1730 


PRINT 


1740 


IMAGE 


D. DD 





USING 1560 

3/ , 30 X , " ** EXPERIMENTAL DATA***' 

USING 1580 ; F* 

35X , “RUN#" . X , 6 A 
USING 1600 ; D1S.T1* 

33X , 8A, X , , X , 4A 
USING 1620 

3/, 32X, “RIGHT-HAND HEATER" 

USING 1640 

31 X, "TEMPERATURE PROFILE" 

USING 1660 

/ , 25X , "COL 1 , 3X , "COL 2",3X,"C0L 3",3X,"C0L 4",7X,"AVE" 

USING 1680 

23X, " " 

USING 1700 5 T<1) ,T(2) ,T<3) ,R(1) 

17X, "ROW 1",X, " ! ", X,DDD.DD,2X,DDD.DD,2X,DDD,DD,2X, •'******% 2X, M 

USING 1720 
23X , " J " , 33X . " ! " 

USING 1740 ; T (5) , T (6) , T (7) , T ( 8) , R (2) 

17X, "ROW 2". X, " I " , X , DDD . DD , 2X , DDD . DD, 2X , DDD. DD, 2X , DDD. DD, 2X , " ! " 



" 2X 



2X. DD 
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: "50 

1 7 

: "30 

Z . DD 

i 7°0 

1 BOO 

1810 

1220 

D. DD 

1 820 

1840 

1850 

1860 

D. DD 

1870 

1880 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

1960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2130 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

% 



PRINT USING 1 ~©0 
IMAGE 23 X , " : " . 33X . " I " 

=, R I NT USING 1 730 : T ( 9 > . T < 1 0 ) . T ( 1 1 ) . T ( 1 2 ) . R ; 3 ) 

IMAGE 1 7X . " ROW 3” . X. " i " . X. DDD. DD . 2X . DDD. DD. 2X . DDD . DD . 2X , DDD. DD. 2X . 1 ' . 2* , 2D 

PRINT USING 1800 
IMAGE 23X . " I ", 33X , " ! " 

PRINT USING 1820 : T * 1 3 ) . T ( 1 4 ) . T ( 1 5 ) . T < 1 6 > , R < 4 ) 

IMAGE 17X. "ROW 4" . X. " i " , X. DDD. DD, 2X . DDD. DD, 2X . DDD. DD, 2X , DDD. DD. 2X. " ! ” . 2> . DC 

PRINT USING 1840 
IMAGE 23X , ” ! " , 33X , " \ " 

PRINT USING 1860 ; T ( 1 7) , T ( 1 8) , T < 19) , T (20) . R (5) 

IMAGE 17X, "ROW 5" . X, " ! ", X, DDD. DD , 2X , DDD . DD. 2 X, DDD. DD, 2X, DDD. DD , 2X , " ! ” . 2X . DD 
PRINT USING 1880 
PRINT USING 1900 

IMAGE /, 16X. "NOTE: TEMPERATURES SHOWN ARE IN DEGREES FAHRENHEIT" 

PRINT USING 1920 

IMAGE 22 X , " AND ARE ARRANGED AS SEEN FROM WITHIN THE BED" 

PRINT USING 1940 

IMAGE 21/-. 17X, "NOTE: PLOTTED TEMPERATURES ARE AVERAGES OF EACH" 

PRINT USING 1960 

IMAGE 23X, "HORIZONTAL ROW" 

PRINT CHRS (12) 

i 

! ***PLOT TEMP GRAPH (R— H ) tt 

i 

PRINTER IS 2 

PRINT USING 2030 ; F* 

IMAGE 4 / , 1 OX , " RUN# " , X , 6A 
PRINT USING 2050 
IMAGE 3/ 

PRINTER IS 701,100 
Z9=IP(R(3) ) 

Z8-RMDCZ9. 5) 

IF Z8*0 THEN 2120 
Z9=Z9-1 
GOTO 2080 
Z6-Z9-15 
Z7*Z9+1 5 
G CLEAR 8 SETGU 
LOCATE 25. 118.20.80 
! FRAME 

SCALE Z6.Z7.0, 10 
FXD 0 

LAXES 1 f 1 . Z6 , 0, 5, 1 
Z5*Z7— 1 6 
MOVE Z5 f —3 
DEG 

LDIR O 3 LORG 4 
LABEL "TEMP (DEG F) • 

Z4«Z6-3 

MOVE Z4,5 

LDIR 90 9 LORG 4 

LABEL "VERT POSITION (IN)" 

Z3-Z6-+-1 4 
MOVE Z3, 12 
LDIR 0 0 LORG 4 
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2320 LhBEl. " R I GHT — t-iAND HEATER" 

zzzo 22= lo- : 4 



2240 MOVE 22,11 
2250 LOIR 0 G’ L2RG 4 
Z2o0 i-ADEL " TEM C vs DEPTH 
2270 3ETUU 
2280 V=9 

2290 FOR 1=1 TO 5 
2400 MOVE R ( I > . Y 
2410 Y«Y-2 



2420 LORG 5 
2430 LABEL 'T 1 
2440 NEXT I 
2450 COPY 

2460 ! ' t x *D I SPLAY HEATER TEMP SNAPSHOT (L-H>*»* 



2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

2600 

2610 

2620 

2630 

D.DD 

2640 

2650 

2660 

2670 

D.DD 

2680 

2690 

2700 

2710 

D.DD 

2720 

2730 

2740 

2750 

D.DD 

2760 

2770 

2780 

2790 

D.DD 

2800 

2810 

2B20 

2B30 

2840 

2B50 

2S60 



PRINT USING 2490 

IMAGE 3/. 30X, " < ^EXPERIMENTAL DATA**' 1 
PRINT USING 2510 ; F$ 

IMAGE 35X, "RUN#" , X, 6A 
PRINT USING 2530 ; D1*,T1* 

IMAGE 33X , BA X , 4 A 

PRINT USING 2550 

IMAGE 3/,33X, "LEFT-HAND HEATER" 

PRINT USING 2570 

IMAGE 31 X, "TEMPERATURE PROFILE" 

PRINT USING 2590 

IMAGE / , 25X , "COL 1",3X,"C0L 2",3X,"C0L 3",3X,"C0L 4",7X, " AVE" 

PRINT USING 2610 

IMAGE 23X, " " 

PRINT USING 2630 ; T (21 ) . T (22) , T (23) „ T (24 > , L ( 1 ) 

IMAGE 1 7X , "ROW 1 " , X , " i " , X , DDD . DD . 2X , DDD . DD , 2X , DDD. DD. 2X , DDD. DD, 2X , " ! " , 2X , DD 

PRINT USING 2650 
IMAGE 23 X, " I " , 33X , " ! " 

PRINT USING 2670 ; T (25) . T (25) , T (27) , T (28) , L (2) 

IMAGE 1 7X , "ROW 2",X, " ! " . X . DDD , DD , 2X , DDD . DD , 2X , DDD . DD , 2X , DDD . DD, 2X , " J ",2X, DD 

PRINT USING 2690 
IMAGE 23X, " ! " , 33X, " ! " 

PRINT USING 2710 ; T (29) , T (30) , T (31 > . T (32) , L (3) 

IMAGE 1 7X • "ROW 3",X, " I " , X , DDD. DD, 2X , DDD . DD* 2X , DDD. DD, 2X. DDD. DD, 2X , " : " , 2X . DD 

PRINT USING 2730 
IMAGE 23X, " i ", 33X, " J " 

PRINT USING 2750 s T (33) , T (34) , T (35) , T (36) , L (4 ) 

IMAGE 17X, "ROW 4". X, " ! " , X , DDD . DD , 2X , DDD . DD , 2X , DDD . DD , 2X , DDD . DD , 2X , " ! " ,2X, DD 

PRINT USING 2770 
IMAGE 23X, " ! " ,33X, " S M 

PRINT USING 2790 ; T (37) , T (38) , T (39) , T (40) , L (5) 

IMAGE 17X, "ROW 5", X, " ! " , X , DDD . DD , 2X , DDD . DD , 2X , DDD . DD , 2X , DDD . DD , 2X , " ! " f 2X,DD 
PRINT USING 2810 

IMAGE 23X," " 

PRINT USING 2B30 

IMAGE / , 1 6X , "NOTE : TEMPERATURES SHOWN ARE IN DEGREES FAHRENHEIT" 

PRINT USING 2850 

IMAGE 22 X, "AND ARE ARRANGED AS SEEN FROM WITHIN THE BED" 

PRINT USING 2870 
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2970 in AGE 21 ' , 1 7X . "NOTE: -‘LOTTED TEMPERATURES ARE AVERAGES QF EACH*' 

23 SO PRINT USING 29^0 

23P0 IMAGE 23 X. "HORIZONTAL ROW*’ 

2-00 PRINT 2 HR 4 (12) 



2 C 10 1 

2-20 1 u='LGT TEMP GRAPH (L-H) ** " 

2-30 ' 



2°Au 

2950 

2 Q 6U 

2-70 
2^30 
2990 
3000 
30 1 0 
3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3 1 00 
3110 
3120 
3130 
3140 
3150 
3160 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
3280 
3290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
3370 
3380 
3390 
3400 
3410 
3420 
3430 
3440 
3450 
3460 



PRINTER IS 2 

PRINT USING 2960 : F* 

IMAGE 4/, 1 OX, "RUN# ",X,6A 
PRINT USING 2980 
IMAGE 3/ 

PRINTER IS 701,100 
Z9*IP (L (3) ) 

Z8=RMD CZ9, 5) 

IF Z8=*0 THEN 3050 

Z9=Z9-1 

GOTO 3010 

Z6=Z9-15 

Z7=Z9+ 1 5 

GCLEAR 6* SETGU 

LOCATE 25, 118, 20, 80 

! FRAME 

SCALE Z6, Z7, 0, 10 
FXD 0 

LAXES 1, 1, Z6, 0,5, 1 

Z5=Z7-16 

MOVE Z5,-3 

DEG 

LDIR 0 8 LORG 4 

LABEL M TEMP (DEG F) - 

Z4=Z6— 3 

MOVE Z4 , 5 

LDIR 90 8 LORG 4 

LABEL M VERT POSITION <IN>" 

Z3=Z6+1 4 

MOVE Z3, 12 

LDIR 0 8 LORG 4 

LABEL ” LEFT-HAND HEATER” 

Z2-Z6+14 
MOVE Z2,ll 
LDIR 0 8 LORG 4 
LABEL ’’TEMP vs DEPTH” 

SETUU 

Y=9 

FOR 1=1 TO 5 

MOVE L < I ) , Y 

Y=Y-2 

LORG 5 

LABEL 

NEXT I 

COPY 

i PRINT PAGE HEADING 
PRINT USING 3410 

IMAGE 3/, 30X, ” XtEXPER I MENTAL DATA* * " 
PRINT USING 3430 ; F % 

IMAGE 35X, ’’RUN#" , X, 6A 
PRINT USING 3450 ; D 1S,T1S 
IMAGE 33 X , 8A, X , " — ” , X , 4 A 
PRINT USING 3470 
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34^0 IMAGE 2 * . 34 X . " FLU I D I ZED BED " 

3 a 30 PRINT USING 3490 

3490 IMAGE 31 X, "TEMPERATURE PROFILE" 

7500 1 READ BED SIZE 
3510 01 *=FS C 1 , 2 3 
7320 IF Cl *="12" THEN 4590 

3570 IF Gl*=" 10" THEN 4220 

3540 IF G1 *="08" THEN 3090 

3350 IF B9=l THEN 3850 

3560 • 

7570 ! * *D I SPLAY 6 IN BED TEMP SNAPSHOT** 

3580 | 



3590 


PRINT 


USING 


3600 


3600 


IMAGE 


2/ , 27X 


II m T r - • • 


36 1 0 


PRINT 


USING 


3620 


3620 


IMAGE 


27X , " { 


" . 25X , " : " 


3630 


PRINT 


USING 


3640 


3640 


IMAGE 


27X , •• ! 


" . 23 X , " : " 


3650 


PRINT 


USING 


3660 ; T (56) . T (53) . T (50) 


3660 


IMAGE 


27 X . ■■ J 


", X , DDD. DD, 2X , DDD. DD, 2X , DDD. DD , 2X . " J " 


3670 


PRINT 


USING 


3680 


3680 


IMAGE 


27X , " : 


" . 25X , " ! " 


3690 


PRINT 


USING 


3700 


3700 


IMAGE 


27 X , " : 


" . 25X , " : " 


3710 


PRINT 


USING 


3720 s T (51 ) 


3720 


IMAGE 


27X , " ! 


", 17.X, DDD. DD,2X, " ! " 


3730 


PRINT 


USING 


3740 


3740 


IMAGE 


27 X , " J 


" , 25X , " { " 


3730 


PRINT 


USING 


3760 


3760 


IMAGE 


27 X , ■■ : 


" , 25X , " ! " 


3770 


PRINT 


USING 


3780 : T (57) , T (54) , T (52) 


3780 


IMAGE 


27 X , " ! 


" , X , DDD. DD , 2X . DDD . DD, 2X , DDD. DD, 2X , " ! " 


3790 


PRINT 


USING 


3800 


3800 


IMAGE 


27X , " ! 


" , 23 X , " } " 


3810 


PRINT 


USING 


3830 


3820 


PRINT 


USING 


3840 


3830 


IMAGE 


27 X , " ! 


" , 25X • " ; " 


3840 


IMAGE 


27 X , 


— " 


3850 


! CALCULATE 


AVERAGE BED TEMP 


3860 


B1=(T(50)+T( 


51 )+T (S2)+T(53)+T(54)+T (56)+T (57) ) /7 


3870 


GOTO 4900 




3880 


! 






3890 


i 


* *DISPLAY 0 IN BED TEMP SNAPSHOT** , 


3900 


IF B9*l THEN 


l 4180 


3910 


i 






3920 


PRINT 


USING 


3930 










O 7 ■J'KJ 


1 IrHut 


Ui/ • aijA 


* 


3940 


PRINT 


USING 


3930 


3930 


IMAGE 


25X, " ! 


" , 29 X , " ! " 


3960 


PRINT 


USING 


3970 


3970 


IMAGE 


23 X , " : 


" . 29 X , " : - 


3980 


PRINT 


USING 


3990 ; T ( 56 ) , T ( 53 ) , T ( 50 ) 


3^990 


IMAGE 


23X, " J 


",3X,DDD.DD,2X, DDD. DD,2X . DDD. DD, 2X , " ! 


4000 


PRINT 


USING 


4010 


4010 


IMAGE 


23X , " : 


" , 29X , " ; " 


4020 


PRINT 


USING 


4030 


4030 


IMAGE 


25X, - ! 


" . 29 X , " ! " 


4040 


PRINT 


USING 


4050 ; T (51 ) 


4050 


IMAGE 


23X, " ! 


" . 21 X , DDD. DD, 2X , " } " 


4060 


PRINT 


USING 


4070 
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-* v.; i.J 


IMAGE 


25 X , " : 


" , 2°X , " ; 


-030 


PRINT 


USING 


4090 


-o=?o 


I MAGE 


25 X . '• : 


” « 2 C ’X . " : 


- : c»*:> 


PRINT 


USING 


411 0 : t 


4 1 1 0 


IMAGE 


25 X , •’ ; 


" , 5X , DDD 


4 : 20 


PR I NT 


USING 


4130 


4 1 30 


IMAGE 


25X , ” ! 


".2PX, " : 


4 1 40 


PRINT 


USING 


4 1 50 


4150 


IMAGE 


25X , •• ; 


" , 29X , " ; 


4160 


PRINT 


USING 


4170 


4 170 


IMAGE 


25 X , 




4 1 BO 


! CALCULATE 


AVERAGE ] 


4190 


B 1 = < T ( 50 ) +T ( 


51 ) +T <52 



‘ C !7 ■ , T ( 54 ) , 7 ( 52 



:x, n : " 



4 200 
4210 
4220 
4230 
4240 



GOTO 4900 

i 

* * *D I SPLAY 10 IN BED TEMP SNAPSHOT* * * 
IF 99=1 THEN 4530 



4250 


PRINT 


USING 


4260 


4260 


IMAGE 


2/ , 23 


X, " 


4270 


PRINT 


USING 


4280 


4280 


IMAGE 


23X, " 


: •• , 33: 


4290 


PRINT 


USING 


4300 


4300 


IMAGE 


23X « " 


: •• . 33: 


4310 


PRINT 


USING 


4320 


4320 


IMAGE 


23X, " 


! M ,X,] 


4330 


PRINT 


USING 


4340 


4340 


IMAGE 


23 X . " 


! " 4 33! 


4350 


PRINT 


USING 


4360 


4360 


IMAGE 


23X , ” 


! " , 33! 


4370 


PRINT 


USING 


4380 


4300 


IMAGE 


23X , " 


: ",x,i 


4390 


PRINT 


USING 


4400 


4400 


IMAGE 


23X , ” 


1 ", 33! 


4410 


PRINT 


USING 


4420 


4420 


IMAGE 


23X, •* 


! “ 4 33! 


4430 


PRINT 


USING 


4440 


4440 


IMAGE 


23X , " 


: " • x , ) 


4450 


PRINT 


USING 


4460 


4460 


IMAGE 


23 X , " 


! 11 4 33! 


4470 


PRINT 


USING 


4480 



44S0 

4490 

4 500 

4510 

4 520 

4530 

4540 

4550 

4560 

4370 

4580 

4 590 

4600 

4610 

+T ( 63 

4 620 



'X , " ! M 



33X , 1 



IMAGE 23 X , " 

PRINT USING 4500 

IMAGE 23 X, " H 

PRINT USING 4520 

IMAGE / , 15X , “NOTE; TEMPERATURES SHOWN ARE IN DEGREES FAHRENHEIT" 

PRINT USING 4540 

IMAGE 21 X , "AND ARE ARRANGED AS SEEN FROM ABOVE THE BED" 

! CALCULATE AVERAGE BED TEMP 

Bl*<T<50)+T<5l>+T<52)+T<53)+T(54)+T<56)+T<57)+T<58)+T(59)+T<60> >/10 
GOTO 4950 

i 

! * * *DI SPLAY 12 IN BED TEMP SNAPSHOT*** 

! CALCULATE AVERAGE BED TEMP 

B1 = <T<50)+T(31)+T(52>+T(33)-»*T(54>+T<56)+T(37)+T(58)+T(59)+T<60>+T(61)^T<62) 

) ) /13 

IF B9 = 1 THEN 4950 



4630 ! 

4640 PRINT USING 4650 
4650 IMAGE 2/,19X," 
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4©©0 


PRINT 


USING 


4 e 7 o 




4© 70 


I MAGE 


1 9 X . " : 


" . 4 1 X . " 


: " 


4oS0 


PRINT 


USING 


4o 90 




4 r g 0 


IMAGE 


i*y. " : 


" . 4 1 >. , " 


: ’* 


4700 


pc- 1 NT 


USING 


4^10 : 


1 i 61) , T (58) . T (5o) . 7 (53) . T (50) 


4710 


I MAGE 


ipx, •• : 


M . X, DDD 


. DD, 2X. DDD. DD. 2X. DDD. DD. 2X . DDD. DD. 2X . DDD. DD. 2X 


4^20 


PRINT 


USING 


4730 




4730 


IMAGE 


IPX, " : 


" , 4 1 X , M 


: " 


4740 


PRINT 


USING 


4750 




4 750 


IMAGE 


1 9 x , •• : 


" . 4 1 X . M 


! " 


4760 


PRINT 


USING 


4770 ; 


T (62) , T (59) . T (51 ) 


4770 


IMAGE 


19X, " ! 


M , X , DDD 


. DD, 2X, DDD. DD, 18X , DDD. DD, 2X, " ! " 


4780 


PRINT 


USING 


4790 




4790 


IMAGE 


1 9 X . " ! 


" . 4 1 X , " 


; ** 


4000 


PRINT 


USING 


4810 




4810 


IMAGE 


i9x, •• : 


M , 41 X, M 


! " 


4820 


PRINT 


USING 


4830 ; 


T (63) , T (60) , T (57) , T (54) , T (52) 


4830 


IMAGE 


i9x, •• : 


" , X , DDD 


. DD, 2X , DDD. DD. 2X, DDD. DD, 2X . DDD. DD, 2X , DDD. DD, 2X 


4840 


PRINT 


USING 


4850 




4850 


IMAGE 


1 9X . " ! 


" , 4 1 X , " 


• " 


4860 


PRINT 


USING 


4870 




4870 


IMAGE 


19X. M : 


" . 41 X, " 


• II 


4880 


PRINT 

IMAGE 


USING 

1 O Y •• _ 


4890 




V \,1 


A ' A • 






4900 


IF B9«l THEN 


1 4960 




4910 


PRINT 


USING 


4920 




4920 


IMAGE 


/. 15X , 


'•NOTE: 


TEMPERATURES SHOWN ARE IN DEGREES FAHRENHEIT " 


4930 


PRINT 


USING 


4940 




4 940 


IMAGE 


21 X," AND ARE 


ARRANGED AS SEEN FROM ABOVE THE BED" 



4950 J 

4960 ! * * ^CALCULATE RESULTS** 

4970 ! 

4980 • CALCULATE AVERAGE AIR INLET TEMP 
4990 P1»(T (47) +T(48) +T<49) ) /3 
5000 ! CALCULATE AIR OUTLET TEMP 
5010 P2=T (55) 

5020 ! CALCULATE RH HEATER GUARD TEMP 
5030 G1*(T(41)+T( 42) **-T (43) ) /3 
5040 ! PRINT PAGE HEADING 
5050 PRINT CHR* (12) 

5060 PRINT USING 5070 

5070 IMAGE 3/ • 30X , " * *CALCULATED RESULTS** M 

5080 PRINT USING 5090 ; FS 

5090 IMAGE 35X, "RUN #",X,6A 

5100 PRINT USING 5110 ; D1*,T1S 

5110 IMAGE 33X, 8A. X, "- M , X,4A 

5120 PRINT USING 5130 

5130 IMAGE 2/.32X, "RIGHT-HAND HEATER “ 

5140 ! PRINT AVERAGE HEATER TEMPERATURE 
5150 PRINT USING 5160 : HI 

5160 IMAGE 2/, 15X, M AVERAGE HEATER TEMPERATURE - " , 12X , DDD . DD , X , M DEG F" 

5170 • CALCULATE AND DISPLAY ELECTRICAL ENERGY INTO HEATER 
5180 E1«V1 *11*3.4121 

5190 ! CALCULATE AND DISPLAY LOSS OUT OF BACK OF HEATER 
5200 G3* ( HI -G1 )*. 1703163 

5210 ! CALCULATE CONVECTION HEAT TRANSFER COEFFICIENT FROM GUARD TO ATMOSPHERE 
5220 Cl =G3/ ( . 347222* (Gl-Al ) ) 

5230 ! CALCULATE LOSS FROM SIDES OF HEATER TO ATMOSHPERE 
5240 S3-C1 *. 069444* <T (72) -Al) 

5250 ! CALCULATE LOSS FROM TOP STRIP TO ATMOSPHERE 
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52 rO J2=C 1 * . 0205333* C T oo)-Al; 

i— - Al.CJwA"E LOSS ” R 0 M BOTTOM STP I r TG ATMOSPHERE 

7230 37 «C! *.0205773* “ • ” 1 ’ -41 ) 

7 C' 1 *"- CALCULATE TG7Al_ _OSSES 37 ATMOSPHERE 

3300 r-4=G3— 53**-U2* , S7 

2710 1 CALCULATE HEAT FLUX THRU RH PLATE 
5320 P3=E1-H4 
5330 V**0 

3340 1 CALCULATE SURFACE TEMPERATURE OF PLATE 
3350 P5=H1-P3*. 0208333/ ( (232-. 032* (Hi -70) ) * . 347222) 

5360 1 CALCULATE HEAT TRANSFER COEFFICIENT FOR PLATE TO BED 
5330 P7=P3/ < . 347222* (P5-B1 ) ) 



5330 ! CALCULATE LOSS FROM SIDE STRIPS TO BED 
5390 C3=P7*. 069444* (T (73) -B1 ) 

5400 * CALCULATE LOSS FROM TOP STRIP TO BED 
5410 S9»F7t. 0208333* (T (69) -Bl) 

5420 ! CALCULATE LOSS FROM BOTTOM STRIP TO BED 
5430 S5»P7*. 0208333* (T (70) -Bl ) 

5440 1 CALCULATE TOTAL LOSSES INTO BED 
5450 H3=C3+S9+S5 

5460 ! CALCULATE NEW HEAT FLUX THRU PLATE 
5470 X9*E1-H4-H3 



5480 V9*V9+1 



5490 ! COMPARE NEW PLATE HEAT FLUX WITH OLD AND REPEAT CALCULATION IF DIFFERENCE 
EXCEEDS 0.01 
5500 X8*ABS (P3-X9) 

5510 IF X8>=. 01 THEN P3-X9 0 GOTO 5340 

5520 J PRINT RESULTS 

5530 PRINT USING 5540 ; P5 

5540 IMAGE 15X, "AVERAGE PLATE SURFACE TEMPERATURE *= " « 5X , DDD. DD . X . " DEG F” 

5550 PRINT USING 5560 : Bl 

5560 IMAGE 1 5X .’* AVERAGE BED TEMPERATURE - ",1 5X , DDD . DD , X , "DEG F M 
5570 PRINT USING 5580 : El 

5580 IMAGE 15X, "ELECTRICAL ENERGY INTO HEATER *" , 8X , DDD. DD, X , "BTU/HR" 

5590 PRINT USING 5600 : H4 

5600 IMAGE 15X,**L0SS FROM HTR BACKING TO ATM «*\ 9X , DDD. DD , X , "BTU/HR" 

5610 PRINT USING 5620 ; Cl 

5620 IMAGE 15X."h(HTR BACKING-ATM) , 8X , DDD . DD , X , "BTU/HR— FT~2-DEG F" 

5630 FRINT USING 5640 ; H3 

5640 IMAGE 15X, f, L0SS FROM HTR PERIMETER TO BED *" , 7X , DDD. DD, X , "BTU/HR” 

5650 PRINT USING 5660 ; F3 

5660 IMAGE 15X,"HEAT FLUX THRU PLATE 1 7X , DDD. DD , X , " BTU/HR " 

5670 PRINT USING 5680 ; P7 
5680 IMAGE 15X , "h (PLATE-BED) 

5690 ! 



14X,DDD.DD, X, "BTU/HR-FT^C-DES F" 



5700 ! CALCULATE AND DISPLAY RESULTS FOR LEFT-HAND HEATER 

5710 ! 

5720 PRINT USING 5730 

5730 IMAGE 2/ . 33X , "LEFT-HAND HEATER** 

5740 ! CALCULATE LH HEATER GUARD TEMPERATURE AVE 
5750 G2«(T(44)+T(45>+T(46) > /3 

5760 ! CALCULATE ELECTRICAL ENERGY INTO HEATER 
5770 E2-V2*I2X3.4121 

5780 ! CALCULATE LOSS FROM BACK OF HEATER TO ATMOSPHERE 
5790 G4«(H2-G2) *. 1703163 

5800 ! CALCULATE CONVECTION HEAT TRANSFER COEFFICIENT FROM GUARD TO ATMOSHPERE 
5810 C2*G4/ (.347222* (G2-A1) ) 

5820 ! CALCULATE LOSS FROM SIDES OF HEATER TO ATMOSPHERE 
5830 S4-C2*. 069444* (T (78) -Al) 

5840 ! CALCULATE LOSS FROM TOP STRIP TO ATMOSPHERE 
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5850 
5 a 90 
3570 
3850 
55 Q 0 
5900 
5 Q 1 0 
5920 
5930 
5 Q 40 
5950 
5960 
5970 
5980 
5990 
6000 
6010 
6020 
6030 
6040 
6050 
6060 
0 . 0 1 
6070 
6080 
6090 
6100 
6110 
6120 
6130 
6140 
6150 
6160 
6170 
6180 
6190 
6200 
6210 
6220 
6230 
6240 
6250 
6260 
6270 
6280 
6290 
6300 
63 1 0 
6320 
6330 
6340 
6350 
6360 
6370 
6380 
6390 
6400 
6410 
6420 
6430 



U3— C2* . 0206333* ( T ( ~o > —A 1 ) 

1 CALCULATE i_3S5 FROM BOTTOM STRIP TO ATMOSPHERE 
38=U3 

1 CALCULATE TOTAL LOSSES TO ATMOSPHERE 

i-i c =04 —'Z.a, 

* CALCULATE HEAT FLUX THRU LH PLATE 
P4=E2-H6 

1 CALCULATE SURFACE TEMPERATURE OF PLATE 
F‘6»H2-F’4* .0208333/ ( (232-. 032* (H2-70) ) * . 347222) 

! CALCULATE HEAT TRANSFER COEFFICIENT FOR PLATE TO BED 
P0=P4/ ( . 347222* (P6-B1 ) ) 

? CALCULATE LOSS FROM SIDE STRIPS TO BED 
C4-P8* . 069444* (T (79) -Bl ) 

' CALCULATE LOSS FROM TOP STRIP TO BED 
U 1 *P8 * . 0208333* ( T ( 77 ) -B 1 ) 

* CALCULATE LOSS FROM BOTTOM STRIP TO BED 
S6-U 1 

! CALCULATE TOTAL LOSSES TO BED 
H5-C4-HJ1+S6 

‘ CALCULATE NEW HEAT FLUX THRU PLATE 
X7-E2-H6-H5 

! COMPARE NEW HEAT FLUX TO OLD AND REPEAT CALCULATION IF DIFFERENCE EXCEEDS 
X6-ABS (P4-X7) 

IF X 6 >= . 01 THEN P4=X7 £ GOTO 5920 
« PRINT RESULTS 
PRINT USING 6110 : H2 

IMAGE 2/, 15X, “AVERAGE HEATER TEMPERATURE , 12X , DDD, DD, X , "DEG F” 

PRINT USING 6130 ; P6 

IMAGE 15X, "AVERAGE PLATE SURFACE TEMPERATURE , 5X , DDD. DD, X , "DEG F" 

PRINT USING 6150 : Bl 

IMAGE 1 5X , *' AVERAGE BED TEMPERATURE * " , 1 5X , DDD. DD, X , “DEG G" 

PRINT USING 6170 : E2 

IMAGE 15X, "ELECTRICAL ENERGY INTO HEATER , 8X , DDD. DD, X , "BTU/HR" 

PRINT USING 6190 : H6 

IMAGE 1 5X , "LOSS FROM HTR BACKING TO ATM *" , 9X , DDD. DD, X , "BTU/HR'* 

PRINT USING 6210 ; C2 

IMAGE 15X, **h (HTR BACKING-ATM) »",8X, DDD. DD, X, " BTU/HR— FT^2— DEG F" 

PRINT USING 6230 s H5 

IMAGE 15X , "LOSS FROM HTR PERIMETER TO BED **' « 7X , DDD. DD, X , "BTU/HR** 

PRINT USING 6250 ; P4 

IMAGE 15X,"HEAT FLUX THRU PLATE - “ , 1 7X , DDD. DD, X , "BTU/HR M 
PRINT USING 6270 ; P8 

IMAGE 15X, "h (PLATE-BED) * M , 1 4X , DDD . DD , X , " BTU/HR— FT~2— DEG F M 

i 

! CALCULATE ENERGY BALANCE VALUES 

i 

! CALCULATE AIR DENSITY 
X5» 1.325*29, 92/ (P2+459. 69) 

! CALCULATE ENERGY CARRIED OUT OF BED BY AIR STREAM 
A5-F2*X5*60*.241* (P2-P1) 

! CALCULATE LOSS THRU FRONT AND BACK FACES 



L6-.448*X* (T (75) — 


T (74) ) - 


* 










! CALCULATE 
A6-A5+L6 


TOTAL 


ENERGY 


LEAVING ! 


BED 






! CALCULATE 
B3-E1-H4 


TOTAL 


ENERGY 


INTO 


BED 


FROM 


RH 


HEATER 


f CALCULATE 
B4-E2-H6 


TOTAL 


ENERGY 


INTO 


BED 


FROM 


LH 


HEATER 


) CALCULATE 


TOTAL 


ENERGY 


INTO 


BED 


FROM 


ALL 


. SOURCES 
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c*440 A7sE.“-**p4 

r'-.'O ! CA-CJLATE CONVECTION HEAT TRANSFER COEFFICIENT F/p FACES-ATMQSPHEF 
3 A 0 O jo=Lo/ ;C4*>;» <T (74 ) -A1 , > 

-4-0 ' CALCULATE CONVECTION HEAT TRANSFER COEFFICIENT F/P FACES-BED 
o430 D?=L6/ (C4*X* (Bl-T (75) ; ) 

©4*0 * CALCULATE SUPERFICIAL VELOCITY 
to 500 U=F2/ Co* X/ 144) /60 

©510 ' CALCULATE SUPERFICIAL MASS VELOCITY 
6500 G=U*X5*3600 

©530 ! CALCULATE PARTICLE REYNOLDS NUMBER 

6540 R=G/3600*. 0009166667/. 000000396/32. 174 

o550 ! PRINT RESULTS 

6560 PRINT USING 6570 

6570 IMAGE 2/ . 34 X , "ENERGY BALANCE" 

6580 PRINT USING 65<?0 ; Pi 

65P0 IMAGE 2/, 15X, "AVERAGE AIR INLET TEMPERATURE **" , 9X , DDD. DD, X , "DEG F" 
6600 PRINT USING 6610 ; P2 

6610 IMAGE 1 5X , "AIR OUTLET TEMPERATURE = " , 1 6X , DDD . DD , X , "DEG F" 

6620 PRINT USING 6630 ; Bl 

6630 IMAGE 15X, "AVERAGE BED TEMPERATURE « " , 15X , DDD . DD , X , " DEG F" 

6640 PRINT USING 6650 : A5 

6650 IMAGE l5X,"q(AIR OUT) =" , 27 X , DDD. DD, X , "BTU/HR" 

6660 PRINT USING 6670 ; L 6 

6670 IMAGE 15X."q(L0SS F/R WALL) 21 X , DDD. DD, X , "BTU/HR" 

6680 PRINT USING 6690 ; A 6 

6690 IMAGE l5X,"q (TOTAL OUT OF BED) *=" , 18X , DDD. DD . X , " BTU/HR" 

6700 PRINT USING 6710 ; A7 

6710 IMAGE 15X,"q (TOTAL INTO BED) »•' , 20X , DDD. DD, X , "BTU/HR" 

6720 PRINT USING 6730 ; U 

6730 IMAGE 15X, "SUPERFICIAL VELOCITY 17X, DDD. DD, X, "FT/SEC" 

6740 PRINT USING 6750 ; G 

6750 IMAGE 15X , " SUPERFIC I AL MASS VELOCITY »- , 7X, DDD. DD, X , "LeM/HR-FT^" 
6760 PRINT USING 6770 ; R 

6770 IMAGE 1 5X , "PARTICLE REYNOLDS NUMBER * M , 20X , DDD . DD 
6780 END 
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